Changeset 4352 for LMDZ6/trunk
- Timestamp:
- Nov 20, 2022, 12:01:18 PM (3 years ago)
- Location:
- LMDZ6/trunk/libf/phylmd
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/clesphys.h
r4245 r4352 7 7 ! 8 8 !..include cles_phys.h 9 ! 9 10 ! threshold on zrel_mount to activate SSO schemes 11 REAL zrel_mount_t 10 12 INTEGER iflag_cycle_diurne 11 13 LOGICAL soil_model,new_oliq,ok_orodr,ok_orolf … … 122 124 & , co2_ppm0 & 123 125 !FC 124 & , Cd_frein 126 & , Cd_frein, zrel_mount_t & 125 127 & , ecrit_LES & 126 128 & , ecrit_ins, ecrit_hf, ecrit_day & -
LMDZ6/trunk/libf/phylmd/conf_phys_m.F90
r4250 r4352 226 226 LOGICAL, SAVE :: ok_lic_cond_omp 227 227 ! 228 REAL, SAVE :: zrel_mount_t_omp 228 229 INTEGER, SAVE :: iflag_cycle_diurne_omp 229 230 LOGICAL, SAVE :: soil_model_omp,new_oliq_omp … … 900 901 CALL getin('ok_orolf', ok_orolf_omp) 901 902 903 904 !Config Key = zrel_mount_t 905 !Config Desc = zrel_mount_t 906 !Config Def = 0. 907 !Config Help = Connais pas ! 908 zrel_mount_t_omp = 0. 909 CALL getin('zrel_mount_t', zrel_mount_t_omp) 910 911 902 912 !Config Key = ok_limitvrai 903 913 !Config Desc = Force la lecture de la bonne annee … … 2426 2436 ok_orodr = ok_orodr_omp 2427 2437 ok_orolf = ok_orolf_omp 2438 zrel_mount_t=zrel_mount_t_omp 2428 2439 ok_limitvrai = ok_limitvrai_omp 2429 2440 nbapp_rad = nbapp_rad_omp … … 2851 2862 WRITE(lunout,*) ' ok_orodr=',ok_orodr 2852 2863 WRITE(lunout,*) ' ok_orolf=',ok_orolf 2864 WRITE(lunout,*) ' zrel_mount_t=',zrel_mount_t 2853 2865 WRITE(lunout,*) ' ok_limitvrai=',ok_limitvrai 2854 2866 WRITE(lunout,*) ' nbapp_rad=',nbapp_rad -
LMDZ6/trunk/libf/phylmd/physiq_mod.F90
r4298 r4352 972 972 973 973 REAL picefra(klon,klev) 974 REAL zrel_mount(klon) 974 975 !IM cf. AM 081204 END 975 976 ! … … 4547 4548 DO i=1,klon 4548 4549 itest(i)=0 4549 ! IF ((zstd(i).gt.10.0)) THEN 4550 IF (((zpic(i)-zmea(i)).GT.100.).AND.(zstd(i).GT.10.0)) THEN 4550 zrel_mount(i)=zstd(i)/(max(zsig(i),1.E-8)*sqrt(cell_area(i))) 4551 !zrel_mount: relative mountain height wrt relief explained by mean slope 4552 ! -> condition on zrel_mount can deactivate the drag on tilted planar terrains 4553 ! such as ice sheets (work by V. Wiener) 4554 IF (((zpic(i)-zmea(i)).GT.100.).AND.(zstd(i).GT.10.0).AND.(zrel_mount(i).GE.zrel_mount_t)) THEN 4551 4555 itest(i)=1 4552 4556 igwd=igwd+1 … … 4597 4601 DO i=1,klon 4598 4602 itest(i)=0 4599 IF ((zpic(i)-zmea(i)).GT.100.) THEN 4603 !zrel_mount: relative mountain height wrt relief explained by mean slope 4604 ! -> condition on zrel_mount can deactivate the lifting on tilted planar terrains 4605 ! such as ice sheets (work by V. Wiener) 4606 zrel_mount(i)=zstd(i)/(max(zsig(i),1.E-8)*sqrt(cell_area(i))) 4607 IF (((zpic(i)-zmea(i)).GT.100.).AND.(zrel_mount(i).GE.zrel_mount_t)) THEN 4600 4608 itest(i)=1 4601 4609 igwd=igwd+1
Note: See TracChangeset
for help on using the changeset viewer.