Changeset 4458 for LMDZ6/trunk/libf/phylmd/physiq_mod.F90
- Timestamp:
- Mar 8, 2023, 3:54:13 PM (15 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
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.