| 1 | c------------------------------------------------------ |
|---|
| 2 | c Fill some parameters in the 'include' files |
|---|
| 3 | c >> Do part of the job previously done by phyetat0.F |
|---|
| 4 | c >> Complete list of parameters is found in tabfi.F |
|---|
| 5 | c------------------------------------------------------ |
|---|
| 6 | c |
|---|
| 7 | c Values are defined in the module_model_constants.F WRF routine |
|---|
| 8 | c |
|---|
| 9 | ! in 'comcstfi.h' |
|---|
| 10 | omeg=womeg |
|---|
| 11 | mugaz=wmugaz |
|---|
| 12 | print*,"check: rad,cpp,g,r,rcp,daysec,omeg,mugaz" |
|---|
| 13 | print*,rad,cpp,g,r,rcp,daysec,omeg,mugaz |
|---|
| 14 | |
|---|
| 15 | ! in 'planet.h' |
|---|
| 16 | year_day=wyear_day |
|---|
| 17 | periheli=wperiheli |
|---|
| 18 | aphelie=waphelie |
|---|
| 19 | peri_day=wperi_day |
|---|
| 20 | obliquit=wobliquit |
|---|
| 21 | emin_turb=wemin_turb |
|---|
| 22 | lmixmin=wlmixmin |
|---|
| 23 | print*,"check: year_day,periheli,aphelie,peri_day,obliquit" |
|---|
| 24 | print*,year_day,periheli,aphelie,peri_day,obliquit |
|---|
| 25 | print*,"check: emin_turb,lmixmin" |
|---|
| 26 | print*,emin_turb,lmixmin |
|---|
| 27 | |
|---|
| 28 | ! in 'surfdat.h' |
|---|
| 29 | emissiv=wemissiv |
|---|
| 30 | emisice(1)=wemissiceN |
|---|
| 31 | emisice(2)=wemissiceS |
|---|
| 32 | albedice(1)=walbediceN |
|---|
| 33 | albedice(2)=walbediceS |
|---|
| 34 | iceradius(1)=wiceradiusN |
|---|
| 35 | iceradius(2)=wiceradiusS |
|---|
| 36 | dtemisice(1)=wdtemisiceN |
|---|
| 37 | dtemisice(2)=wdtemisiceS |
|---|
| 38 | z0_default=wz0 |
|---|
| 39 | print*,"check: z0def,emissiv,emisice,albedice,iceradius,dtemisice" |
|---|
| 40 | print*,z0_default,emissiv,emisice,albedice,iceradius,dtemisice |
|---|
| 41 | |
|---|
| 42 | c |
|---|
| 43 | c Values are defined in the WPS processing |
|---|
| 44 | c |
|---|
| 45 | albedodat(:)=walbedodat(:) |
|---|
| 46 | !!!!! ***WRF inertiedat was moved, new physics !! |
|---|
| 47 | !inertiedat(:)=winertiedat(:) |
|---|
| 48 | phisfi(:)=wphisfi(:) |
|---|
| 49 | print*,"check: albedodat(1),phisfi(1)" |
|---|
| 50 | print*,albedodat(1),phisfi(1) |
|---|
| 51 | print*,"check: albedodat(end),phisfi(end)" |
|---|
| 52 | print*,albedodat(ngrid),phisfi(ngrid) |
|---|
| 53 | |
|---|
| 54 | ! NB: usually, gravity wave scheme is useless in mesoscale modeling |
|---|
| 55 | ! NB: we however keep the option for coarse grid case ... |
|---|
| 56 | zmea(:)=wzmea(:) |
|---|
| 57 | zstd(:)=wzstd(:) |
|---|
| 58 | zsig(:)=wzsig(:) |
|---|
| 59 | zgam(:)=wzgam(:) |
|---|
| 60 | zthe(:)=wzthe(:) |
|---|
| 61 | print*,"check: gw param" |
|---|
| 62 | print*,zmea(1),zmea(ngrid) |
|---|
| 63 | print*,zstd(1),zstd(ngrid) |
|---|
| 64 | print*,zsig(1),zsig(ngrid) |
|---|
| 65 | print*,zgam(1),zgam(ngrid) |
|---|
| 66 | print*,zthe(1),zthe(ngrid) |
|---|
| 67 | |
|---|
| 68 | ! |
|---|
| 69 | ! in slope.h |
|---|
| 70 | ! |
|---|
| 71 | theta_sl(:)=wtheta(:) |
|---|
| 72 | psi_sl(:)=wpsi(:) |
|---|
| 73 | print*,"check: theta_sl(1),psi_sl(1)" |
|---|
| 74 | print*,theta_sl(1),psi_sl(1) |
|---|
| 75 | print*,"check: theta_sl(end),psi_sl(end)" |
|---|
| 76 | print*,theta_sl(ngrid),psi_sl(ngrid) |
|---|
| 77 | |
|---|
| 78 | ! |
|---|
| 79 | ! in comsoil.h |
|---|
| 80 | ! |
|---|
| 81 | volcapa=wvolcapa |
|---|
| 82 | print*,"check: volcapa" |
|---|
| 83 | print*,volcapa |
|---|
| 84 | |
|---|
| 85 | ! |
|---|
| 86 | ! in surfdat.h |
|---|
| 87 | ! |
|---|
| 88 | z0(:) = wz0tab(:) |
|---|
| 89 | print*,"check: z0" |
|---|
| 90 | print*,z0(1),z0(ngrid) |
|---|
| 91 | |
|---|
| 92 | c***************************************************** |
|---|
| 93 | c Since it comes from WRF settings, we have to |
|---|
| 94 | c fill dtphys in the include file |
|---|
| 95 | c It must be set now, because it is used afterwards |
|---|
| 96 | c***************************************************** |
|---|
| 97 | dtphys=wdt*ptimestep |
|---|
| 98 | print*,'Physical timestep (s) ',dtphys |
|---|
