Changeset 3229

Timestamp:

Feb 21, 2018, 3:21:00 PM (7 years ago)

Author:

acozic

Message:

Add a patch to solve a problem when we use inca and we don't find O3 in start.nc file

File:

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

LMDZ5/branches/IPSLCM5A2.1/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