Changeset 2720 for LMDZ5/branches/testing/libf/phylmd/cosp/cosp_output_write_mod.F90

Timestamp:

Nov 30, 2016, 1:28:41 PM (8 years ago)

Author:

Laurent Fairhead

Message:

Merged trunk changes r2664:2719 into testing branch

Location:

LMDZ5/branches/testing

Files:

• . (modified) (1 prop)
• libf/phylmd/cosp/cosp_output_write_mod.F90 (modified) (2 diffs)

LMDZ5/branches/testing/libf/phylmd/cosp/cosp_output_write_mod.F90

@@ -372 +372 @@
  endif
 
+    where(modis%Optical_Thickness_vs_ReffIce == R_UNDEF) &
+          modis%Optical_Thickness_vs_ReffIce = Cosp_fill_value
+
+    where(modis%Optical_Thickness_vs_ReffLiq == R_UNDEF) &
+          modis%Optical_Thickness_vs_ReffLiq = Cosp_fill_value
+
+    do icl=1,7
+   CALL histwrite3d_cosp(o_crimodis, &
+     modis%Optical_Thickness_vs_ReffIce(:,icl,:),nvertReffIce,icl)
+   CALL histwrite3d_cosp(o_crlmodis, &
+     modis%Optical_Thickness_vs_ReffLiq(:,icl,:),nvertReffLiq,icl)
+    enddo
+
  IF(.NOT.cosp_varsdefined) THEN
 !$OMP MASTER
@@ -521 +534 @@
       ELSE IF (nvertsave.eq.nvertmisr(iff)) THEN
           klevs=MISR_N_CTH
-           nam_axvert="cth16"
+          nam_axvert="cth16"
+      ELSE IF (nvertsave.eq.nvertReffIce(iff)) THEN
+          klevs= numMODISReffIceBins
+          nam_axvert="ReffIce"
+      ELSE IF (nvertsave.eq.nvertReffLiq(iff)) THEN
+          klevs= numMODISReffLiqBins
+          nam_axvert="ReffLiq"
       ELSE
            klevs=Nlevout