Changeset 3331 for LMDZ5/branches/IPSLCM5A2.1_ISO/libf/dyn3dmem/dynetat0_loc.f90

Timestamp:

May 22, 2018, 4:30:21 PM (7 years ago)

Author:

acozic

Message:

Add modification for isotopes

File:

• LMDZ5/branches/IPSLCM5A2.1_ISO/libf/dyn3dmem/dynetat0_loc.f90 (modified) (3 diffs)

LMDZ5/branches/IPSLCM5A2.1_ISO/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
@@ -171 +159 @@
     WRITE(lunout,*)"         It is hence initialized to zero"
     q(ijb_u:ije_u,:,iq)=0.
+  END DO !DO iq=1,nqtot
+
    !--- CRisi: for isotops, theoretical initialization using very simplified
-   !           Rayleigh distillation las.
+   !           Rayleigh distillation.
+  DO iq=1,nqtot     
+    if (ok_iso_verif) then
+       write(*,*) 'dynetat0_loc 167: iq,iso_num(iq),iqpere(iq)=',iq,iso_num(iq),iqpere(iq)
+    endif
     IF(ok_isotopes.AND.iso_num(iq)>0) THEN
       IF(zone_num(iq)==0) q(:,:,iq)=q(:,:,iqpere(iq))*tnat(iso_num(iq))        &
@@ -178 +172 @@
       IF(zone_num(iq)==1) q(:,:,iq)=q(:,:,iqiso(iso_indnum(iq),phase_num(iq)))
     END IF
-  END DO
+  END DO !DO iq=1,nqtot
+
+  if (ok_iso_verif) then
+         call check_isotopes(q,ijb_u,ije_u,'dynetat0_loc 171')
+  endif !if (ok_iso_verif) then
+
   DEALLOCATE(q_glo)
   CALL err(NF90_CLOSE(fID),"close",fichnom)