Changeset 2955 for LMDZ5/trunk/libf/phylmd/cosp/cosp_stats.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_stats.F90 (modified) (5 diffs)

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

@@ -121 +121 @@
                                            Nlr,vgrid%zl,vgrid%zu,temp_out)
             temp_c(:,:) = temp_out(:,1,:)
+            stlidar%proftemp = temp_c                                     !TIBO
+            where (stlidar%proftemp  < 150.) stlidar%proftemp   = R_UNDEF !TIBO
+            where (stlidar%proftemp  > 350.) stlidar%proftemp   = R_UNDEF !TIBO
 
             call cosp_change_vertical_grid(Npoints,1,Nlevels,gbx%zlev,gbx%zlev_half,ph_in, &
@@ -131 +134 @@
                             ,LIDAR_UNDEF,ok_lidar_cfad &
                             ,stlidar%cfad_sr,stlidar%srbval &
-                            ,LIDAR_NCAT,stlidar%lidarcld,stlidar%lidarcldphase &
-                            ,stlidar%cldlayer,stlidar%cldlayerphase,stlidar%lidarcldtmp &
-                            ,stlidar%parasolrefl)
+                            ,LIDAR_NCAT,LIDAR_NTYPE,stlidar%lidarcld,stlidar%lidarcldtype & !OPAQ
+                            ,stlidar%lidarcldphase,stlidar%cldlayer,stlidar%cldtype &       !OPAQ
+                            ,stlidar%cldlayerphase,stlidar%lidarcldtmp &                    !OPAQ
+                            ,stlidar%parasolrefl,vgrid%z,stlidar%profSR)                    !OPAQ !TIBO
         endif
 
@@ -140 +144 @@
                                     temp_c,betatot_out,betaperptot_out,betamol_c,Ze_out, &
                                     stradar%lidar_only_freq_cloud,stradar%radar_lidar_tcc)
-        deallocate(temp_out,temp_c,betaperptot_out)
+        deallocate(temp_in,temp_out,temp_c,betaperptot_out) !TIBO +temp_in
 
         ! Deallocate arrays at coarse resolution
@@ -154 +158 @@
                         ,LIDAR_UNDEF,ok_lidar_cfad &
                         ,stlidar%cfad_sr,stlidar%srbval &
-                        ,LIDAR_NCAT,stlidar%lidarcld,stlidar%lidarcldphase,stlidar%cldlayer,stlidar%cldlayerphase &
-                        ,stlidar%lidarcldtmp,stlidar%parasolrefl)
+                        ,LIDAR_NCAT,LIDAR_NTYPE,stlidar%lidarcld,stlidar%lidarcldtype & !OPAQ
+                        ,stlidar%lidarcldphase,stlidar%cldlayer,stlidar%cldtype &       !OPAQ
+                        ,stlidar%cldlayerphase,stlidar%lidarcldtmp &                    !OPAQ
+                        ,stlidar%parasolrefl,vgrid%z,stlidar%profSR)                    !OPAQ !TIBO
         !++++++++++++ Lidar-only cloud amount and lidar&radar total cloud mount ++++++++++++++++
         if (cfg%Lradar_sim.and.cfg%Llidar_sim) call cosp_lidar_only_cloud(Npoints,Ncolumns,Nlr, &
@@ -163 +169 @@
    ! Replace undef
    where (stlidar%cfad_sr   == LIDAR_UNDEF) stlidar%cfad_sr   = R_UNDEF
+   where (stlidar%profSR   == LIDAR_UNDEF) stlidar%profSR   = R_UNDEF !TIBO
    where (stlidar%lidarcld  == LIDAR_UNDEF) stlidar%lidarcld  = R_UNDEF
+   where (stlidar%lidarcldtype  == LIDAR_UNDEF) stlidar%lidarcldtype  = R_UNDEF !OPAQ
    where (stlidar%cldlayer  == LIDAR_UNDEF) stlidar%cldlayer  = R_UNDEF
+   where (stlidar%cldtype  == LIDAR_UNDEF) stlidar%cldtype  = R_UNDEF           !OPAQ
    where (stlidar%parasolrefl == LIDAR_UNDEF) stlidar%parasolrefl = R_UNDEF
    where (stlidar%cldlayerphase  == LIDAR_UNDEF) stlidar%cldlayerphase  = R_UNDEF