Changeset 4263 for LMDZ6/trunk/libf/dyn3dmem

Timestamp:

Sep 21, 2022, 9:35:28 AM (3 years ago)

Author:

dcugnet

Message:

Fixes for INCO, CO2i AND REPROBUS, mostly because some sections are specific to type_trac=="lmdz",
which is not always equivalent to ANY(types_trac=='lmdz).
Also force the water phases to get tracers(*)%component='lmdz' so that nqo can be correctly computed.

File:

• LMDZ6/trunk/libf/dyn3dmem/dynetat0_loc.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/dyn3dmem/dynetat0_loc.F90 ¶

@@ -46 +46 @@
   REAL,             ALLOCATABLE :: ucov_glo(:,:),    q_glo(:,:), phis_glo(:)
   REAL,             ALLOCATABLE :: teta_glo(:,:)
+  LOGICAL :: lOldHNO3
 !-------------------------------------------------------------------------------
   modname="dynetat0_loc"
@@ -151 +152 @@
 !--- Tracers
   ALLOCATE(q_glo(ip1jmp1,llm))
+#ifdef REPROBUS
+!--- Old convention: HNO3 = g+l and HNO3_g = g ; new convention: HNO3tot = g+l and HNO3 = g
+  lOldHNO3 = NF90_INQ_VARID(fID,var,vID), "inq", 'HNO3tot') /= NF90_NoErr
+#endif
   DO iq=1,nqtot
     var = tracers(iq)%name
@@ -166 +171 @@
       CALL msg('Tracer <O3> is missing => initialized to <OX>', modname)
       CALL get_var2( 'OX' , q_glo); q(ijb_u:ije_u,:,iq)=q_glo(ijb_u:ije_u,:)
-    !--------------------------------------------------------------------------------------------------------------------------
+#endif
+    !--------------------------------------------------------------------------------------------------------------------------
+#ifdef REPROBUS
+    ELSE IF(lOldHNO3 .AND. var == 'HNO3') THEN                                           !=== REPROBUS: HNO3_g INSTEAD OF HNO3
+      CALL msg('Tracer <HNO3_g> is used for <HNO3> (old convention)', modname)
+      CALL get_var2( 'HNO3_g', q_glo); q(ijb_u:ije_u,:,iq)=q_glo(ijb_u:ije_u,:)
+    ELSE IF(lOldHNO3 .AND. var == 'HNO3tot') THEN                                        !=== REPROBUS: HNO3 INSTEAD OF HNO3tot
+      CALL msg('Tracer <HNO3> is used for <HNO3tot> (old convention)', modname)
+      CALL get_var2( 'HNO3', q_glo); q(ijb_u:ije_u,:,iq)=q_glo(ijb_u:ije_u,:)
 #endif
     ELSE IF(tracers(iq)%iso_iGroup == iH2O .AND. niso > 0) THEN                          !=== WATER ISOTOPES