Changeset 2955 for LMDZ5/trunk/libf/phylmd/cosp/cosp_output_write_mod.F90

Timestamp:

Jul 17, 2017, 5:52:31 PM (7 years ago)

Author:

musat

Message:

Ajout modifs de Rodrigo Guzman et Marine Bonazolla dans la version COSP/CMIP6
IM

File:

• LMDZ5/trunk/libf/phylmd/cosp/cosp_output_write_mod.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/cosp/cosp_output_write_mod.F90

@@ -104 +104 @@
       stlidar%lidarcld(ip,k)=missing_val
      endif
+     if(stlidar%proftemp(ip,k).eq.R_UNDEF)then !TIBO
+      stlidar%proftemp(ip,k)=missing_val       !TIBO
+     endif                                     !TIBO
      enddo
 
@@ -113 +116 @@
       enddo
      enddo
-   enddo
+
+     do ii= 1,Ncolumns                               !TIBO
+      do ip = 1,Npoints                              !TIBO
+!       if(stlidar%profSR(ip,ii,k).eq.R_UNDEF)then    !TIBO
+!        stlidar%profSR(ip,ii,k)=missing_val          !TIBO
+       if(stlidar%profSR(ip,k,ii).eq.R_UNDEF)then    !TIBO2
+        stlidar%profSR(ip,k,ii)=missing_val          !TIBO2
+       endif                                         !TIBO
+      enddo                                          !TIBO
+     enddo                                           !TIBO
+  enddo
 
   do ip = 1,Npoints
@@ -138 +151 @@
    enddo
 
+   do k = 1,LIDAR_NTYPE                       !OPAQ
+    do ip = 1,Npoints                         !OPAQ
+     if(stlidar%cldtype(ip,k).eq.R_UNDEF)then !OPAQ
+       stlidar%cldtype(ip,k)=missing_val      !OPAQ
+     endif                                    !OPAQ
+    enddo                                     !OPAQ
+   enddo                                      !OPAQ
+
 ! AI 11 / 2015
 
@@ -144 +165 @@
    where(stlidar%cldlayerphase == R_UNDEF) stlidar%cldlayerphase = missing_val
    where(stlidar%lidarcldphase == R_UNDEF) stlidar%lidarcldphase = missing_val
+   where(stlidar%lidarcldtype == R_UNDEF) stlidar%lidarcldtype = missing_val   !OPAQ
    where(stlidar%lidarcldtmp == R_UNDEF) stlidar%lidarcldtmp = missing_val
    
@@ -176 +198 @@
    CALL histwrite3d_cosp(o_clcalipsotmpun,stlidar%lidarcldtmp(:,:,4),nverttemp)
 
+   CALL histwrite2d_cosp(o_clopaquecalipso,stlidar%cldtype(:,1))               !OPAQ
+   CALL histwrite2d_cosp(o_clthincalipso,stlidar%cldtype(:,2))                 !OPAQ
+   CALL histwrite2d_cosp(o_clzopaquecalipso,stlidar%cldtype(:,3))              !OPAQ
+
+   CALL histwrite3d_cosp(o_clcalipsoopaque,stlidar%lidarcldtype(:,:,1),nvert)  !OPAQ
+   CALL histwrite3d_cosp(o_clcalipsothin,stlidar%lidarcldtype(:,:,2),nvert)    !OPAQ
+   CALL histwrite3d_cosp(o_clcalipsozopaque,stlidar%lidarcldtype(:,:,3),nvert) !OPAQ
+   CALL histwrite3d_cosp(o_clcalipsoopacity,stlidar%lidarcldtype(:,:,4),nvert) !OPAQ
+
+   CALL histwrite3d_cosp(o_proftemp,stlidar%proftemp,nvert)                    !TIBO
+
 #ifdef CPP_XIOS
    CALL histwrite4d_cosp(o_cfad_lidarsr532,stlidar%cfad_sr)
+   CALL histwrite4d_cosp(o_profSR,stlidar%profSR)                              !TIBO
 #else
    do icl=1,SR_BINS
       CALL histwrite3d_cosp(o_cfad_lidarsr532,stlidar%cfad_sr(:,icl,:),nvert,icl)
    enddo
-#endif
-
-   CALL histwrite3d_cosp(o_parasol_refl,stlidar%parasolrefl,nvertp)
+   do icl=1,Ncolumns                                                              !TIBO
+      CALL histwrite3d_cosp(o_profSR,stlidar%profSR(:,icl,:),nvert,icl)           !TIBO
+   enddo                                                                          !TIBO
+     CALL histwrite3d_cosp(o_parasol_refl,stlidar%parasolrefl,nvertp)
+#endif
+
    do k=1,PARASOL_NREFL
     do ip=1, Npoints