Changeset 2254 for trunk/LMDZ.GENERIC/libf/phystd/inifis_mod.F90

Timestamp:

Mar 10, 2020, 4:59:49 PM (5 years ago)

Author:

jvatant

Message:

Fix a bug in the 2-layers aerosol scheme in aeropacity.F90
-> It appeared that the normalization for total column optical depth was a

wrong way to go, the actual strato and tropo values weren't matcihing the
input values, and one of the main consequence was to strongly underestimate
the actual strato optical depth compared to input value in most cases.

-> Now we ensure the normalization to the input values for both layers, while

keeping a dependence in dP for each GCM layer.

-> Also add a sanity check ensuring pres_top_tropo > pres_bottom_strato.
--JVO

File:

• trunk/LMDZ.GENERIC/libf/phystd/inifis_mod.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phystd/inifis_mod.F90

@@ -594 +594 @@
      write(*,*)" pres_bottom_strato = ",pres_bottom_strato
 
+     ! Sanity check
+     if (pres_bottom_strato .gt. pres_top_tropo) then
+       print*,'Error : TWOLAY AEROSOL, Please ensure that in callphys.def'
+       print*,'you have pres_top_tropo > pres_bottom_strato !'
+       stop
+     endif
+
      write(*,*)"TWOLAY AEROSOL: pres_top_strato? in pa"
      pres_top_strato=100.0