Changeset 2790

Timestamp:

Feb 8, 2017, 4:05:38 PM (8 years ago)

Author:

acozic

Message:

Modification for INCA
If there is no value for O3 concentration in start.nc, we decide to read value of OX (OX=O3 + O1D + O) instead of initialized to zero.

File:

• LMDZ5/trunk/libf/dyn3dmem/dynetat0_loc.f90 (modified) (1 diff)

LMDZ5/trunk/libf/dyn3dmem/dynetat0_loc.f90

@@ -153 +153 @@
   DO iq=1,nqtot
     var=tname(iq)
+#ifdef INCA
+    IF (var .eq. "O3" ) THEN 
+       IF(NF90_INQ_VARID(fID,var,vID) == NF90_NoErr) THEN
+          CALL get_var2(var,q_glo); q(ijb_u:ije_u,:,iq)=q_glo(ijb_u:ije_u,:); CYCLE
+       ELSE
+          WRITE(lunout,*) 'Tracer O3 is missing - it is initialized to OX'
+          IF(NF90_INQ_VARID(fID,"OX",vID) == NF90_NoErr) THEN
+             CALL get_var2("OX",q_glo); q(ijb_u:ije_u,:,iq)=q_glo(ijb_u:ije_u,:); CYCLE
+          ENDIF
+       ENDIF
+    ENDIF
+#endif
     IF(NF90_INQ_VARID(fID,var,vID)==NF90_NoErr) THEN
       CALL get_var2(var,q_glo); q(ijb_u:ije_u,:,iq)=q_glo(ijb_u:ije_u,:); CYCLE