Hello all,

 

I am running RegCM model ,after executing  terrain.x ,I find the file size of DOMAIN.INFO and DOMAIN.CTL in Input directory  is zero byte.

Below attach is logfile terrain scripts.  

LAMCON Projection Map ( 25S-45N & 40E-115E) 

I may want to know the procedure of calculating truelatL and truelatH .

 

Please help me to come over this issue.

 

rm -f *.o terrain
ifort -tpp7 -g -cm -w -w90 -w95 -convert big_endian -c -L../../Commons/env/liblinuxamd64/ -lnetcdf -Vaxlib terrain.f
ifort -tpp7 -g -cm -w -w90 -w95 -convert big_endian -o terrain terrain.o -L../../Commons/env/liblinuxamd64/ -lnetcdf -Vaxlib
 ntypec=           2
 iy=         259
 jx=         278
 ds=   30.00000   
 clat=   10.00000   
 clon=   77.50000   
 rin=   1.500000   
 iproj=LAMCON
 
 
 **************************************************
 FILE ALREADY EXISTS:  ../../Input/DOMAIN.INFO                          
 Do you want to overwrite the existing file? [y/n/q]
y
 
 
 **************************************************
 FILE ALREADY EXISTS:  ../../Input/DOMAIN.CTL                           
 Do you want to overwrite the existing file? [y/n/q]
y
 ***** Terrain resolution (min):    2.000000   
 after calling SETUP
 PSX,PSI1 =    1.396263      0.3490658   
 PSX,PSI1 =    1.396263      0.3490658   
 XN,TRUELATL,TRUELATH =   0.8067775       35.00000       70.00000   
 after calling MAP PROJECTION
 xminlat,xmaxlat,xminlon,xmaxlon,ntypec=     -17.16     42.73     32.43    122.81         2
 after calling MXMNLL
    4877089 terrain heights read from land use volume
 after calling RDLDTR
 after calling XYOBSLL
  rin,ds(m) =  -0.116E+11  -0.718E+10
  rin,ds(m) =  -0.116E+11  -0.718E+10
 after calling ANAL2
 after calling SURF
 LANDUSE MASK           0 exceed the limit

 

domain.param file
 parameter(iproj='LAMCON') ! Map projection (LAMCON,NORMER,POLSTR,ROTMER)
      parameter(iy=259)          ! number of grid points in y direction (i)
      parameter(jx=278)          ! number of grid points in x direction (j)
      parameter(kz=23)          ! number of vertical levels (14, 18 or 23)
      parameter(nsg=1)          ! number of subgrids in 1 grid (one direction)

      parameter(ds=30.0)        ! grid point separation in km
      parameter(ptop= 5.0)      ! pressure of model top in cb.

      parameter(clat= 10.00)    ! Central latitude (North positive)
      parameter(clon= 77.50)    ! Central longitude (East positive)

      parameter(plat=clat)      ! Pole latitude  (used for ROTMER)
      parameter(plon=clon)      ! Pole longitude (used for ROTMER)

      parameter(truelatL= 35.)  ! LAMCON true latitude (low  latitude side)
      parameter(truelatH= 70.)  ! LAMCON true latitude (high latitude side)

      parameter(ntypec=2)      ! Global topo/landuse data resolution (min)
      parameter(ntypec_s=2)    ! Global topo/landuse data resolution (min)
      parameter(h2opct=75.)     ! Minimum H2O percent to be considered water

      parameter(ifanal =.true.) ! Type of interpolation for topography
      parameter(smthbdy=.false.)! Smoothing of boundaries
      parameter(lakadj=.false.) ! Adjust Great Lakes levels (for US domain)

      parameter(igrads=1)       ! Output GrADS control file?
      parameter(ibigend=1)      ! 1=big-endian; (always 1)
      parameter(ibyte=1)        ! For direct access open statements (1 or 4)

      parameter(FUDGE_LND=.false.)   ! LANDUSE adjustments?
      parameter(FUDGE_LND_s=.false.) ! Subgrid LANDUSE adjustments?

      parameter(FUDGE_TEX=.false.)   ! TEXTURE adjustments?
      parameter(FUDGE_TEX_s=.false.) ! Subgrid TEXTURE adjustments?

      parameter(filout='../../Input/DOMAIN.INFO') ! Location/name of output
      parameter(filctl='../../Input/DOMAIN.CTL')  ! Location/name of GrADS ctl

Thanks
--
Rajesh Kumar Maurya
M.Sc(Tech)- Geophysics(Meteorology); Junior Research Fellow & Research Scholar
Department of Geophysics
Banaras Hindu University
Varanasi 221005