Changeset 1928 for LMDZ4/branches/LMDZ4_AR5/libf/cosp/phys_cosp.F90

Timestamp:

Jan 16, 2014, 1:43:34 PM (10 years ago)

Author:

idelkadi

Message:

• Correction du bug dans l'interface avec la physique : On ne veut pas que la distinction entre les nuages convectifs et stratiformes soit prise en compte dans les calculs Cosp. Dans l'interface avec la physique, sont passees en entree pour Cosp, les quantites totales (stratiforme + convective) du contenus en eau et d'autres variables. La fraction nuageuse convective calculee dans la physique est passee en entree pour Cosp dans la version buggee. Elle est remise a 0 dans cette version corrigee.

• Mise a jour pour ISCCP : la fraction d'ensoleillement calculee par LMDZ est passee en argument pour Isccp

• Rajouts d'autre variables diagnostiques : sunlit (=0 si nuit et =1 si jour) parasol_crefl (reflectance integree)

File:

• LMDZ4/branches/LMDZ4_AR5/libf/cosp/phys_cosp.F90 (modified) (6 diffs)

LMDZ4/branches/LMDZ4_AR5/libf/cosp/phys_cosp.F90

@@ -1 +1 @@
 ! Simulateur COSP : Cfmip Observation Simulator Package
+
 ! ISCCP, Radar (QuickBeam), Lidar et Parasol (ACTSIM), MISR, RTTOVS
 !Idelkadi Abderrahmane Aout-Septembre 2009
-
 
   subroutine phys_cosp( itap,dtime,freq_cosp, &
                         ok_mensuelCOSP,ok_journeCOSP,ok_hfCOSP, &
                         ecrit_mth,ecrit_day,ecrit_hf, &
-                        Nptslmdz,Nlevlmdz,lon,lat, presnivs,overlaplmdz, &
+                        Nptslmdz,Nlevlmdz,lon,lat, presnivs,overlaplmdz,sunlit, &
                         ref_liq,ref_ice,fracTerLic,u_wind,v_wind,phis,phi,ph,p,skt,t, &
                         sh,rh,tca,cca,mr_lsliq,mr_lsice,fl_lsrainI,fl_lssnowI, &
@@ -144 +144 @@
                                      zlev,zlev_half,mr_ozone,radliq,radice,dtau_s,dem_s,ref_liq,ref_ice
   real,dimension(Nptslmdz,Nlevlmdz) ::  fl_lsrainI,fl_lssnowI,fl_ccrainI,fl_ccsnowI
-  real,dimension(Nptslmdz)        :: lon,lat,skt,fracTerLic,u_wind,v_wind,phis
+  real,dimension(Nptslmdz)        :: lon,lat,skt,fracTerLic,u_wind,v_wind,phis,sunlit    
   real,dimension(Nlevlmdz)        :: presnivs
   integer                         :: itap,k,ip
   real                            :: dtime,freq_cosp
   logical, parameter              :: lCOSP=.FALSE.
+
+  real, dimension(Nptslmdz,PARASOL_NREFL) :: parasolcrefl, Ncref
 
 !
@@ -237 +239 @@
         gbx%q = rh*100.
         gbx%sh = sh
-        gbx%cca = cca !convective_cloud_amount (1)
+        gbx%cca = 0.
         gbx%tca = tca ! total_cloud_amount (1)
         gbx%psfc = ph(:,1) !pression de surface
@@ -253 +255 @@
         gbx%u_wind  = u_wind !eastward_wind (m s-1)
         gbx%v_wind  = v_wind !northward_wind
-! Attention
-        gbx%sunlit  = 1
+
+      do ip = 1, Npoints
+        if (sunlit(ip).le.0.) then
+           gbx%sunlit(ip)=0
+        else
+           gbx%sunlit(ip)=1
+        endif
+      enddo
 
 ! A voir l equivalent LMDZ
@@ -280 +288 @@
      gbx%Reff(:,:,I_CVCLIQ) = ref_liq*1e-6
      gbx%Reff(:,:,I_CVCICE) = ref_ice*1e-6
-
         ! ISCCP simulator
         gbx%dtau_s   = dtau_s
@@ -423 +430 @@
      if(isccp%boxptop(ip,ii).eq.-1.E+30)then
        isccp%boxptop(ip,ii)=0.
+     endif
+    enddo
+   enddo
+
+   do k=1,PARASOL_NREFL
+    do ii=1, Npoints
+     if (stlidar%cldlayer(ii,4).gt.0.01) then
+        parasolcrefl(ii,k)=(stlidar%parasolrefl(ii,k)-0.03*(1.-stlidar%cldlayer(ii,4)))/ &
+                   stlidar%cldlayer(ii,4)
+        Ncref(ii,k) = 1.
+     else
+        parasolcrefl(ii,k)=0.
+        Ncref(ii,k) = 0.
      endif
     enddo