Changeset 4458 for LMDZ6/trunk/libf/phylmd
- Timestamp:
- Mar 8, 2023, 3:54:13 PM (16 months ago)
- Location:
- LMDZ6/trunk/libf/phylmd
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/clesphys.h
r4448 r4458 8 8 !..include cles_phys.h 9 9 10 ! threshold on zrel_mountto activate SSO schemes11 REAL zrel_ mount_t10 ! threshold on to activate SSO schemes 11 REAL zrel_oro_t, zpmm_orodr_t, zpmm_orolf_t, zstd_orodr_t 12 12 INTEGER iflag_cycle_diurne 13 13 LOGICAL soil_model,new_oliq,ok_orodr,ok_orolf … … 123 123 & , co2_ppm0 & 124 124 !FC 125 & , Cd_frein, zrel_mount_t&125 & , Cd_frein,zrel_oro_t,zpmm_orodr_t,zpmm_orolf_t,zstd_orodr_t & 126 126 & , ecrit_LES & 127 127 & , ecrit_ins, ecrit_hf, ecrit_day & -
LMDZ6/trunk/libf/phylmd/conf_phys_m.F90
r4448 r4458 225 225 LOGICAL, SAVE :: ok_lic_cond_omp 226 226 ! 227 REAL, SAVE :: zrel_mount_t_omp 227 REAL, SAVE :: zrel_oro_t_omp, zstd_orodr_t_omp 228 REAL, SAVE :: zpmm_orodr_t_omp, zpmm_orolf_t_omp 228 229 INTEGER, SAVE :: iflag_cycle_diurne_omp 229 230 LOGICAL, SAVE :: soil_model_omp,new_oliq_omp … … 900 901 901 902 902 !Config Key = zrel_ mount_t903 !Config Desc = zrel_ mount_t903 !Config Key = zrel_oro_t 904 !Config Desc = zrel_oro_t 904 905 !Config Def = 0. 905 906 !Config Help = Connais pas ! 906 zrel_mount_t_omp = 0. 907 CALL getin('zrel_mount_t', zrel_mount_t_omp) 908 907 zrel_oro_t_omp = 0. 908 CALL getin('zrel_oro_t', zrel_oro_t_omp) 909 910 !Config Key = zstd_orodr_t 911 !Config Desc = zstd_orodr_t 912 !Config Def = 0. 913 !Config Help = Connais pas ! 914 zstd_orodr_t_omp = 10. 915 CALL getin('zstd_orodr_t', zstd_orodr_t_omp) 916 917 918 !Config Key = zpmm_orodr_t 919 !Config Desc = zpmm_orodr_t 920 !Config Def = 0. 921 !Config Help = Connais pas ! 922 zpmm_orodr_t_omp = 100. 923 CALL getin('zpmm_orodr_t', zpmm_orodr_t_omp) 924 925 926 !Config Key = zpmm_orolf_t 927 !Config Desc = zpmm_orolf_t 928 !Config Def = 0. 929 !Config Help = Connais pas ! 930 zpmm_orolf_t_omp = 100. 931 CALL getin('zpmm_orolf_t', zpmm_orolf_t_omp) 909 932 910 933 !Config Key = ok_limitvrai … … 2412 2435 ok_orodr = ok_orodr_omp 2413 2436 ok_orolf = ok_orolf_omp 2414 zrel_mount_t=zrel_mount_t_omp 2437 zrel_oro_t=zrel_oro_t_omp 2438 zstd_orodr_t=zstd_orodr_t_omp 2439 zpmm_orodr_t=zpmm_orodr_t_omp 2440 zpmm_orolf_t=zpmm_orolf_t_omp 2415 2441 ok_limitvrai = ok_limitvrai_omp 2416 2442 nbapp_rad = nbapp_rad_omp … … 2835 2861 WRITE(lunout,*) ' ok_orodr=',ok_orodr 2836 2862 WRITE(lunout,*) ' ok_orolf=',ok_orolf 2837 WRITE(lunout,*) ' zrel_mount_t=',zrel_mount_t 2863 WRITE(lunout,*) ' zrel_oro_t=',zrel_oro_t 2864 WRITE(lunout,*) ' zstd_orodr_t=',zstd_orodr_t 2865 WRITE(lunout,*) ' zpmm_orodr_t=',zpmm_orodr_t 2866 WRITE(lunout,*) ' zpmm_orolf_t=',zpmm_orolf_t 2838 2867 WRITE(lunout,*) ' ok_limitvrai=',ok_limitvrai 2839 2868 WRITE(lunout,*) ' nbapp_rad=',nbapp_rad -
LMDZ6/trunk/libf/phylmd/physiq_mod.F90
r4449 r4458 975 975 976 976 REAL picefra(klon,klev) 977 REAL zrel_ mount(klon)977 REAL zrel_oro(klon) 978 978 !IM cf. AM 081204 END 979 979 ! … … 4531 4531 DO i=1,klon 4532 4532 itest(i)=0 4533 zrel_ mount(i)=zstd(i)/(max(zsig(i),1.E-8)*sqrt(cell_area(i)))4534 !zrel_ mount: relative mountain height wrt relief explained by mean slope4535 ! -> condition on zrel_ mountcan deactivate the drag on tilted planar terrains4533 zrel_oro(i)=zstd(i)/(max(zsig(i),1.E-8)*sqrt(cell_area(i))) 4534 !zrel_oro: relative mountain height wrt relief explained by mean slope 4535 ! -> condition on zrel_oro can deactivate the drag on tilted planar terrains 4536 4536 ! such as ice sheets (work by V. Wiener) 4537 IF (((zpic(i)-zmea(i)).GT.100.).AND.(zstd(i).GT.10.0).AND.(zrel_mount(i).GE.zrel_mount_t)) THEN 4537 ! zpmm_orodr_t and zstd_orodr_t are activation thresholds set by F. Lott to 4538 ! earn computation time but they are not physical. 4539 IF (((zpic(i)-zmea(i)).GT.zpmm_orodr_t).AND.(zstd(i).GT.zstd_orodr_t).AND.(zrel_oro(i).GE.zrel_oro_t)) THEN 4538 4540 itest(i)=1 4539 4541 igwd=igwd+1 … … 4584 4586 DO i=1,klon 4585 4587 itest(i)=0 4586 !zrel_ mount: relative mountain height wrt relief explained by mean slope4587 ! -> condition on zrel_ mountcan deactivate the lifting on tilted planar terrains4588 !zrel_oro: relative mountain height wrt relief explained by mean slope 4589 ! -> condition on zrel_oro can deactivate the lifting on tilted planar terrains 4588 4590 ! such as ice sheets (work by V. Wiener) 4589 zrel_ mount(i)=zstd(i)/(max(zsig(i),1.E-8)*sqrt(cell_area(i)))4590 IF (((zpic(i)-zmea(i)).GT. 100.).AND.(zrel_mount(i).GE.zrel_mount_t)) THEN4591 zrel_oro(i)=zstd(i)/(max(zsig(i),1.E-8)*sqrt(cell_area(i))) 4592 IF (((zpic(i)-zmea(i)).GT.zpmm_orolf_t).AND.(zrel_oro(i).GE.zrel_oro_t)) THEN 4591 4593 itest(i)=1 4592 4594 igwd=igwd+1 … … 4778 4780 ! car on peut s'attendre a ce que les petites echelles produisent aussi de la TKE 4779 4781 ! Mais attention, cela ne va pas dans le sens de la conservation de l'energie! 4780 IF ( zstd(i).GT.1.0) THEN4782 IF ((zstd(i).GT.1.0) .AND.(zrel_oro(i).GE.zrel_oro_t)) THEN 4781 4783 itest(i)=1 4782 4784 igwd=igwd+1 … … 4790 4792 DO i=1,klon 4791 4793 itest(i)=0 4792 IF (((zpic(i)-zmea(i)).GT. 100.).AND.(zstd(i).GT.10.0)) THEN4794 IF (((zpic(i)-zmea(i)).GT.zpmm_orodr_t).AND.(zstd(i).GT.zstd_orodr_t).AND.(zrel_oro(i).GE.zrel_oro_t)) THEN 4793 4795 itest(i)=1 4794 4796 igwd=igwd+1
Note: See TracChangeset
for help on using the changeset viewer.