Praj-Paramit

hr2dy.f90

integer, parameter :: nx=36, ny=60, nv=27, ndt=8  ! ndt=24hr/3hr=8
real:: data(nx,ny),dvar(nx,ny,nv)
integer:: iyr,iyrs,nmo(12),irec,orec
character*80 infil,oufil
character*8 var(27)
data nmo /31,28,31,30,31,30,31,31,30,31,30,31/
!data var /u10m,v10m,uvdrag,tg,tlef,t2m,q2m,ssw,rsw,tpr,evp,runoff,scv,sena,
!         + flw,fsw,flwd,sina,prcv,psb,zpbl,tgmax,tgmin,t2max,t2min,w10max,ps_min /

   DO iyr=1971,2000
   DO imo=1,12
      iyrs=int(iyr/10)*10
      write(*,*) 'Year = ',iyr, 'YearS = ',iyrs
      write(infil,'(a,i4.4,a,i4.4,i2.2,a)') '/nnas9/jshin/water21/dynds/obs',iyrs,'s_eman_kor/output/SRF.',iyr,imo,'0100'
      write(oufil,'(a,i4.4,i2.2,a)') 'dy_SRF.',iyr,imo,'0100'
      write(*,*) infil
      open(1,file=infil,form='unformatted',access='direct',recl=nx*ny,status='old')
      open(2,file=oufil,form='unformatted',access='direct',recl=nx*ny,status='unknown')

      irec=0
      orec=0
      nt=nmo(imo)
      if(mod(iyr,4).eq.0.and.imo.eq.2)  nt=29

   DO itm=1,nt
       dvar=0.
!       pause

     DO ihr=1,ndt
       DO ivr=1,27
          irec=irec+1
          read(1,rec=irec) data
          dvar(:,:,ivr)=dvar(:,:,ivr)+data(:,:)
!         if(ivr.eq.1) write(*,*) dvar(18,30,ivr),data(18,30)
       ENDDO ! ivr (Variable)
     ENDDO ! ihr

!      write(*,*) dvar(18,30,1),dvar(18,30,1)/ndt
!      write(*,*) dvar(18,30,1)

       DO ivr=1,27
         orec=orec+1
         write(*,*) iyr,imo,nt,irec,orec
         write(2,rec=orec) dvar(:,:,ivr)/ndt
       ENDDO ! ivr (Variable)
    ENDDO ! itm
   ENDDO  ! imo
         close(1)
         close(2)
   ENDDO  ! iyr

end

i=1
ii=1
jj=1
itot=23
*---------------------------
* ii = color
* jj = style
*---------------------------
iitot=7
jjtot=3
*
while (i<=itot)
ii=math_mod(i,iitot)
if(ii=0) ; ii=iitot ;endif
jj=math_nint(i/iitot)+1
say i' 'ii' 'jj
pull dummy
*
'q file 'i
line=sublin(result,2)
ctl=subwrd(line,2)
model=substr(ctl,11,10)
say model

is=(i-1)*12*5+1
ie=i*12*5

*
  x0=8.5
  y0=7.5
  x1 = x0
  x2 = x0+0.5
  x3 = x0+0.8
  y1=y0-i*0.2+0.05
  y2=y0-i*0.2
*
'set ccolor 'ii
'set cstyle 'jj
'set cmark 0'

'set line ' ii ' 'jj
'draw line 'x1' ' y1' ' x2' ' y1
'draw string 'x3 ' 'y2' 'model
*
'set grads off'

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
* script_math_demo.gs
*
* Illustrates the use of the scripting language math funcitons
* Written by M. Fiorino
* Modified by J.M. Adams
* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *

doformat = 1
dotrig = 1
doint = 1
domisc = 1
dovalnum = 1

* FORMATTING NUMBERS
if (doformat)
  print ' '
  print 'Formatting numbers:'
  v = 3.1456
  fmt = '%-6.1f'
  rc = math_format(fmt,v)
  print fmt' of 'v' = 'rc

  fmt = '%-6.2f'
  rc = math_format(fmt,v)
  print fmt' of 'v' = 'rc

  fmt = '%6.2f'
  rc = math_format(fmt,v)
  print fmt'  of 'v' = 'rc

  fmt = '%06.3f'
  rc = math_format(fmt,v)
  print fmt' of 'v' = 'rc

  fmt = '%6.0g'
  rc = math_format(fmt,v)
  print fmt'  of 'v' = 'rc

  fmt = '%12.3e'
  rc = math_format(fmt,v)
  print fmt' of 'v' = 'rc
  print ' '
endif

* TRIG FUNCTIONS
if (dotrig)
  print 'Trigonometric functions:'
  pi  = 3.1415926
  d2r = pi/180
  r2d = 1/d2r
  dang = 45
  ang = dang*d2r
  fmt = '%-6.4f'

  cos  = math_cos(ang)
  cosf = math_format(fmt,cos)
  print 'cos of 'dang' = 'cosf

  sin  = math_sin(ang)
  sinf = math_format(fmt,sin)
  print 'sin of 'dang' = 'sinf

  tan  = math_tan(ang)
  tanf = math_format(fmt,tan)
  print 'tan of 'dang' = 'tanf

  ang = math_acos(cos)
  dang = ang * r2d
  dangf = math_format(fmt,dang)
  print 'acos of 'cosf' = 'dangf

  ang = math_asin(sin)
  dang = ang * r2d
  dangf = math_format(fmt,dang)
  print 'asin of 'sinf' = 'dangf

  ang = math_atan(tan)
  dang = ang * r2d
  dangf = math_format(fmt,dang)
  print 'atan of 'tanf' = 'dangf

  u =  1
  v = -1
  rc = math_atan2(u,v)
  rc2 = rc*r2d
  rc2f = math_format(fmt,rc2)
  print 'atan2 of ('u','v') = 'rc2f

  dang = 30
  ang = dang*d2r

  cosh = math_cosh(ang)
  coshf = math_format(fmt,cosh)
  print 'cosh of 'dang' = 'coshf

  sinh = math_sinh(ang)
  sinhf = math_format(fmt,sinh)
  print 'sinh of 'dang' = 'sinhf

  tanh = math_tanh(ang)
  tanhf = math_format(fmt,tanh)
  print 'tanh of 'dang' = 'tanhf

  demo = dang
  ang = math_acosh(cosh)
  dang = ang * r2d
  dangf = math_format(fmt,dang)
  print 'acosh of 'coshf' = 'dangf' <-- math_acosh() may have a bug ... this value should be 'demo

  ang = math_asinh(sinh)
  dang = ang * r2d
  dangf = math_format(fmt,dang)
  print 'asinh of 'sinhf' = 'dangf

  ang = math_atanh(tanh)
  dang = ang * r2d
  dangf = math_format(fmt,dang)
  print 'atanh of 'tanhf' = 'dangf
  print ' '
endif

* VALNUM
if (dovalnum)
  print 'Evaluating strings to see if they are numbers:'
  print '0 = not a number'
  print '1 = integer'
  print '2 = not an integer'
  num = '3.1455'
  rc = valnum(num)
  print 'valnum of 'num' = 'rc

  num = '3'
  rc = valnum(num)
  print 'valnum of 'num' = 'rc

  num = '3e10'
  rc = valnum(num)
  print 'valnum of 'num' = 'rc

  num = '390210'
  rc = valnum(num)
  print 'valnum of 'num' = 'rc

  num = 'a string'
  rc = valnum(num)
  print 'valnum of 'num' = 'rc
  print ' '
endif

* INTEGER CONVERSION
if (doint) 
  print 'Real-to-integer conversion:'
  v = 3.0
  while(v < 4.0) 
    rc1 = math_nint(v)
    rc2 = math_int(v)
    print 'nint of 'v' = 'rc1'    int of 'v' = 'rc2
    v = v + 0.1
  endwhile
  print ' '
endif

* MISCELLANEOUS FUNCTIONS
if (domisc)
  print 'Exponents:'
  pow = math_pow(2,0.5);
  print '2 raised to the power 0.5 = 'pow
  print ' '

  print 'Exponential function:'
  num = math_exp(1)
  print 'exp(1) = 'num
  print ' '

  print 'Natural logarithm:'
  num = math_log(2.71828)
  print 'log(2.71828) = 'num
  print ' '

  print 'Base 10 logarithm:'
  num = math_log10(2)
  print 'log10(2) = 'num
  pow = math_pow(10,0.30103);
  print '10 raised to the power 0.3103 = 'pow
  print ' '

  print 'Modulo operator fmod:'
  fmod = math_fmod(6,2.25);
  print '6 fmod 2.25 (the remainder when 6 is divided by 2.25) = 'fmod
  print ' '

  print 'Modulo operator mod:'
  fmod = math_mod(6,2.25);
  print '6 mod 2.25 (the integer part of the remainder when 6 is divided by 2.25) = 'fmod
  print ' '

  print 'String operations:'
  s = 'this is a test'
  rc = math_strlen(s)
  print 'length of the string "'s'" = 'rc

  p = 2
  rc = wrdpos(s,p)
  print 'word 'p' of the string "'s'" starts at character 'rc
  print ' '
endif


return

#! /bin/csh

cd /nnas8/hyolee/uahnc/grib

set wgrib=/opt/grads-2.0.a3/grads-1.9b4/bin/wgrib

set year=1979
echo $year

while ($year <= 1979)
  set mon=1
  echo $mon
while ($mon <= 12)
  if ($mon <= 9) then
      set dum=0
  else
      set dum=''
  endif
  echo grib2ctl.pl  ./uahnca.pm${year}${dum}${mon}_pres.grb > ${year}${dum}${mon}.ctl
  grib2ctl.pl ./uahnca.pm${year}${dum}${mon}_pres.grb > ${year}${dum}${mon}.ctl
  gribmap -0 -i  ${year}${dum}${mon}.ctl
  echo gribmap -0 -i  ${year}${dum}${mon}.ctl
@ mon++
end
@ year++
end

#! /bin/csh

         set ifile=/nnas8/hyolee/uahnc/pm/uahnca.pmh9dec
         echo pmfile=`ls -ltr ${ifile}`
         set pmfile=`ls -ltr ${ifile}`
         set pmsize=$pmfile[5]
         echo ${pmsize}
         if ( $pmsize != 693518336 ) then
            echo "File size of ${ifile} is not equal 693518336 "
            stop
         endif