This block data set contains definitions of most physical parameters used in model.
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c---> igrid=0 for A grid (default); c
c---> igrid=1 for B grid (special option) c
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
data igrid/0/
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c I/O OPTIONS for different types of binary formats c
c-----------------------------------------------------------------c
c---> ioptr=0 read data from HP/Convex format c
C on Cray set correct assign options c
c---> ioptr=1 read data from Cray MPP (T3D,T3E) format c
c---> ioptr=2 read data from Cray PVP (C90,J90,Y-MP) format c
c c
c---> ioptw=0 write data to HP/Convex format c
C on Cray set correct assign options c
c---> ioptw=1 write data to Cray MPP (T3D,T3E) format c
c---> ioptw=2 write data to Cray PVP (C90,J90,Y-MP) format c
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
data ioptw,ioptr/0,0/
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
create initialization of environmental state,
c see comments in subroutine tinit_i or tinit_r
c-----------------------------------------------------------------c
c--> initi=1 idealized initialization (tinit_i)
c--> initi=0 realistic initialization from sounding (tinit_r)
c--> lipps=0 Boussinesq
c--> lipps=1 Anelastic, isentropic
c--> lipps=2 Anelastic, Clark-Farley
c--> lipps=3 Anelastic, isothermal
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
data initi/1/
data lipps/2/
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
create model boundaries c
c-----------------------------------------------------------------c
c---> iab,iabth,iabqw - absorber flags for velocity, theta, water c
c---> dxabL- thickensses of left lateral absorbers in x [m] c
c---> dxabR- thickensses of right lateral absorbers in x [m] c
c---> dyab - thickensses of lateral absorbers in y [m] c
c---> zab - base of sponge in vertical (z) [m] c
c---> towxL- absorber time scale in x (left side) [s] c
c---> towxR- absorber time scale in x (right side) [s] c
c---> towy - absorber time scale in y [s] c
c---> towz - absorber time scale in z [s] c
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
data iab,iabth,iabqw/1,1,0/
data dxabL,towxL/1.e0,1.e10/
data dxabR,towxR/1.e0,1.e10/
data dyab ,towy /1.e0,1.e10/
data zab ,towz /12.e3,864000./
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
common block for pressure boundary treatment
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c data e1,e2,e3/1.,0.,0./ ! 0_th order
cc data e1,e2,e3/1.5,-0.5,0./ ! 1_th order
data e1,e2,e3/1.875,-1.250,0.375/ ! 2_nd order; default
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
common block for grag coefficients
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
common/dragc/ dragx(nth),dragy(nth),drgnorm, itd, idrag
data dragx/nth*0./,dragy/nth*0./
data idrag/0/
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
constants for pressure solver; c
c-----------------------------------------------------------------c
c---> itp_,epp_ are max number of iterations and constraint of c
c---> accuracy max|1/rho*div(rho*v)*dt| and run, correspondingly for itp0,epp0 and itp1,epp1; c
c---> see string "epp1=" for other opt. c
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
data itp0,epp0,itp1,epp1/300,1.e-6,100,1.e-5/
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c---> rds - sphere radius (m) c
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c data rds/0.183/ ! cylindrical annulus (inner radius)
data rds/ 6371.22e3/ ! Earth
c data rds/ 637.122e3/ ! 0.1 Earth
C data rds/ 30.e3/ ! small Earth
c data rds/ 696.e6/ ! Sun
c data rds/ 6070.e3/ ! Venus
c data rds/ 3380.e3/ ! Mars
c data rds/66900.e3/ ! Jupiter
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
create basic constants (tt00,th00,pr00,rh00 must be consistent): c
c-----------------------------------------------------------------c
c---> g - gravity (m/s^2) c
c---> rg - dry gas constant (J/kgK) c
c---> tt00 - ground level temperature (K) c
c---> th00 - ground level potential temperature (K) c
c---> pr00 - ground level pressure (Pa) c
c---> rh00 - ground level density (kg/m^3) c
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
common/ctherm/ rg,cp,cap,st,g,th00,tt00,pr00,rh00,u00,v00,u0z,v0z
data g,rg,tt00,th00,pr00,rh00 /9.80665,287.04,
c $ 1025.0,1025.0,1.e5,1025.0/ ! WATER
C $ 1025.0,1025.0,1.e5, 1.0/ ! WATER
$ 303.16,303.16,1.03e5,1.225/ ! Earth
c $ /3.73,189.0 ,190. ,190. ,0.8e3 ,0.022/ ! Mars
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
create environmental profiles (use subroutine tinit_i) c
c-----------------------------------------------------------------c
c---> u00 - wind constant; x component, profile U (m/s) c
c---> v00 - wind constant; y component, profile V (m/s) c
c---> st - stability parameter; inverse scale height (m^-1) c
c---> u0z - vertical wind shear dU/dz (s^-1) c
c---> v0z - vertical wind shear dV/dz (s^-1) c
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
data u00, v00, st, u0z, v0z
$ /20.0, 0.0, 1.0e-5, 0.e-3, 0.e-3/
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
create time-dependent lower boundary; specifics in topolog routine
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
data tt,tend/0.,-1.e10/
data time/0./
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
coriolis force specification c
c-----------------------------------------------------------------c
c---> fcr0 - coriolis parameter (s^-1) c
c---> ang - f=fcr0*sin(ang) for f- or beta-plane approximation c
c---> btpl - beta-plane approximation on or off c
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
common/corio/ ang,fcr0,btpl,fcr2(np,mp),fcr3(np,mp),initprs
c data fcr0,ang,initprs/5.8e-6, 0.,0/ ! Sun
c data fcr0,ang,initprs/0.5985e-6,0.,0/ ! Venus
data fcr0,ang,initprs/1.4584e-4,45.,0/ ! Earth
c data fcr0,ang,initprs/1.4176e-4,0.,0/ ! Mars
c data fcr0,ang,initprs/3.5047e-4,0.,0/ ! Jupiter
c data fcr0,ang,initprs/0.,0.,0/ ! tank
data btpl/0./
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
create algorithm for sperical metric forces c
c-----------------------------------------------------------------c
c---> mtrord - order of approximation (1 or 2), mtrord=3 beta plane
c---> mtrimx - No of iterations (ge.2) for mtrord=2 c
c---> mtrord=mtrord*isphere+1-isphere c
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
common/metrif/ mtrord,mtrimx
data mtrord,mtrimx/2,2/
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
create variables for SGS model c
c-----------------------------------------------------------------c
c turbulent scaling parameters as given by Schumann in c
c Theor. Comp. Fluid Dyn (1991) 2:279-290; note that Prndt=cm/ch c
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c salt in water
c data ceps/.845/,cL/0.845/,cm/.0856/,css/.165/,prndt/685.3973/ !Gill
c salty water + heat
c (salt diffusivity neglected, added viscosity taken into account)
c data ceps/.845/,cL/0.845/,cm/.0856/,css/.165/,prndt/7.458/ !Schumann
data ceps/.845/,cL/0.845/,cm/.0856/,css/.165/,prndt/.42/ !Schumann
c data ceps/.930/,cL/1.e15/,cm/.1000/,css/.180/,prndt/.33/ !Moeng-Nieuwstadt
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c surface layer fluxes for heat/moisture/momentum c
c-----------------------------------------------------------------c
c---> hf00 - sensible heat flux [K m/s] c
c---> qf00 - latent heat flux [kg/kg m/s] c
c---> cdrg - drag coefficient c
c---> inorm- normalization switch c
c---> zi - height of lowest inversion (min Richardson number) c
c---> wstr - scaling vertical velocity in mixin-layer theory c
c---> tstr - sclaing temperature in mixin-layer theory c
c---> qstr - scaling water-vapor mixing ratio in mixin-layer theory
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
data inorm/0/,zi,wstr,tstr,qstr/4*1./
c data hf00/20.e-2/,qf00/12.e-5/,cdrg/0.001/
data hf00/ 0.e-2/,qf00/ 0.e-5/,cdrg/0.000/
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
create HF filter from the upwind option in transport schemes c
c-----------------------------------------------------------------c
c---> liner - flag for activating mpfl'th smoothing
c---> mpfl - decrease to first order every mpfl'th timestep
c---> ampd - advection order operator for mpfl'th smoothing
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
data mpfl,ampd/nth,0.00/
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
construct moist model
c-----------------------------------------------------------------c
create condensation/evaporation constants
common/cmoist/ rv,t00,ee0,hlat
data rv,t00,ee0,hlat/461.,273.16,611.,2.53e6/
create rain constants (ice=0 case)
common/rain/ rac,qctr,rc
data rac,qctr,rc/1.e-3, .005e-3, 2.2/
c data rac,qctr,rc/1.e-3, .5e-1, 2.2/
common/dftemp/dftmx,dftmn,dftav,dftsd
data dftmx,dftmn,dftav,dftsd/4*0./
c time step of microphysical model (no more than 30s for explicit)
common/tmoist/ dtm
data dtm/10./ ! Wojtek recommends dtm < 300 s
cc initialization for ice=1 case
common/rain_p0/ ar,br,cr,dr,er,alphr,betr,gamb1r,gambd1r,anor
common/rain_p1/ dconc,ddisp
common/snow_p0/ as,bs,cs,ds,es,alphs,bets,gamb1s,gambd1s,anos
common/temp_p/ tup,tdn
common/latent/ hlatv,hlats
create parametrs for general precip formulation:
cc mass, terminal velocity, diameter
data ar,br,cr,dr /5.2e2,3.,130.,0.5/
data as,bs,cs,ds /2.5e-2,2.,4.,0.25/
cc collection ef., alpha, beta
data er,alphr,betr /0.8, 1., 2./
data es,alphs,bets /0.2, .3, 3./
cc No
data anor,anos /2*1.e7/
cc latent heats:
data hlatv,hlats /2.53e6,2.84e6/
cc cloud droplet concentration (per cc)
data dconc /200./ ! must be between 50 and 2000
cc limiting temperatures
data tup,tdn /273.,263./
cc gammas:
data gamb1r,gambd1r /6.0,11.7/
data gamb1s,gambd1s /2.0,2.56/
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
create mountain
c-----------------------------------------------------------------c
c amp - maximum height of mountain [m]
c xml - zonal scale of mountain [m] or [radians]
c xml0 - zonal shift from the domain center [m] or [radians]
c yml - meridional scale of mountain [m] or [radians]
c yml0 - meridional shift from the domain center [m] or [radians]
c angle- angle to skew the mountain ridge [radians]
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
data amp,xml,yml,xml0,yml0,angle/2.e3,2*.34906585,3*0./