Changeset 2428 for LMDZ5/trunk/libf/phylmd/cosp/cosp_lidar.F90

Timestamp:

Jan 27, 2016, 10:42:32 AM (8 years ago)

Author:

idelkadi

Message:

Mise a jour du simulateur COSP (passage de la version v3.2 a la version v1.4) :

• mise a jour des sources pour ISCCP, CALIPSO et PARASOL
• prise en compte des changements de phases pour les nuages (Calipso)
• rajout de plusieurs diagnostiques (fraction nuageuse en fonction de la temperature, ...)

http://lmdz.lmd.jussieu.fr/Members/aidelkadi/cosp

File:

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

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

@@ -1 +1 @@
 ! (c) British Crown Copyright 2008, the Met Office.
 ! All rights reserved.
+! $Revision: 88 $, $Date: 2013-11-13 15:08:38 +0100 (mer. 13 nov. 2013) $
+! $URL: http://cfmip-obs-sim.googlecode.com/svn/stable/v1.4.0/cosp_lidar.F90 $
 ! 
 ! Redistribution and use in source and binary forms, with or without modification, are permitted 
@@ -29 +31 @@
 !                               Call lidar_simulator changed (lsca, gbx%cca and depol removed; 
 !                               frac_out changed in sgx%frac_out)
+! Jun 2011 - G. Cesana        - Added betaperp_tot argument
 !
 ! 
@@ -43 +46 @@
 !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 SUBROUTINE COSP_LIDAR(gbx,sgx,sghydro,y)
-  
+
   ! Arguments
   type(cosp_gridbox),intent(in) :: gbx  ! Gridbox info
@@ -49 +52 @@
   type(cosp_sghydro),intent(in) :: sghydro  ! Subgrid info for hydrometeors
   type(cosp_sglidar),intent(inout) :: y ! Subgrid output
-  
+
   ! Local variables 
   integer :: i
   real :: presf(sgx%Npoints, sgx%Nlevels + 1)
-  real :: frac_out(sgx%Npoints, sgx%Nlevels)
   real,dimension(sgx%Npoints, sgx%Nlevels) :: lsca,mr_ll,mr_li,mr_cl,mr_ci
   real,dimension(sgx%Npoints, sgx%Nlevels) :: beta_tot,tau_tot
+  real,dimension(sgx%Npoints, sgx%Nlevels) :: betaperp_tot
   real,dimension(sgx%Npoints, PARASOL_NREFL)  :: refle
-  
-  
+
   presf(:,1:sgx%Nlevels) = gbx%ph
   presf(:,sgx%Nlevels + 1) = 0.0
-!   presf(:,sgx%Nlevels + 1) = gbx%p(:,sgx%Nlevels) - (presf(:,sgx%Nlevels) - gbx%p(:,sgx%Nlevels)) 
   lsca = gbx%tca-gbx%cca
   do i=1,sgx%Ncolumns
@@ -69 +70 @@
       mr_cl(:,:) = sghydro%mr_hydro(:,i,:,I_CVCLIQ)
       mr_ci(:,:) = sghydro%mr_hydro(:,i,:,I_CVCICE)
-      frac_out(:,:) = sgx%frac_out(:,i,:)    
-      call lidar_simulator(sgx%Npoints, sgx%Nlevels, 4 &
-                 , PARASOL_NREFL, LIDAR_UNDEF  &
-                 , gbx%p, presf, gbx%T &
-                 , mr_ll, mr_li, mr_cl, mr_ci &
-                 , gbx%Reff(:,:,I_LSCLIQ), gbx%Reff(:,:,I_LSCICE), gbx%Reff(:,:,I_CVCLIQ), gbx%Reff(:,:,I_CVCICE) &
-                 , frac_out, gbx%lidar_ice_type, y%beta_mol, beta_tot, tau_tot  &
-                 , refle ) ! reflectance
-      
+      call lidar_simulator(sgx%Npoints, sgx%Nlevels, 4, PARASOL_NREFL, LIDAR_UNDEF  &
+                 , gbx%p, presf, gbx%T, mr_ll, mr_li, mr_cl, mr_ci &
+                 , gbx%Reff(:,:,I_LSCLIQ), gbx%Reff(:,:,I_LSCICE) &
+                 , gbx%Reff(:,:,I_CVCLIQ), gbx%Reff(:,:,I_CVCICE) &
+                 , gbx%lidar_ice_type, y%beta_mol, beta_tot &
+                 , betaperp_tot, tau_tot, refle )
+
+      y%betaperp_tot(:,i,:) = betaperp_tot(:,:)
       y%beta_tot(:,i,:) = beta_tot(:,:)
       y%tau_tot(:,i,:)  = tau_tot(:,:)