Changeset 2297 for LMDZ5/trunk/libf/phylmd

Timestamp:

Jun 12, 2015, 11:38:57 PM (9 years ago)

Author:

musat

Message:

RRTM :
Correction bug interface avec rrtm (radlwsw_m.F90)
Ajouter la possibilite d'utilsation d'anciennes proprites optique dans le cas RRTM 2bandes (readaerosol_optic_rrtm.F90)
Correction rrtm (rrtm_rtrn1a_140gp.F90)

COSP :
Modifications pour traiter les valeurs indefines :
mise a 0 en cas d'utilisation de IOIPSL
mise a la valeur recuperee dans .xml en cas d'utilisation de XIOS

Location:

LMDZ5/trunk/libf/phylmd

Files:

• cosp/cosp_output_mod.F90 (modified) (3 diffs)
• cosp/cosp_output_write_mod.F90 (modified) (8 diffs)
• radlwsw_m.F90 (modified) (1 diff)
• rrtm/readaerosol_optic_rrtm.F90 (modified) (2 diffs)
• rrtm/rrtm_rtrn1a_140gp.F90 (modified) (1 diff)

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

@@ -80 +80 @@
            by the ISCCP Simulator","K", (/ ('', i=1, 3) /))
 
-   LOGICAL, SAVE :: cosp_varsdefined = .FALSE. ! ug PAS THREADPRIVATE ET C'EST NORMAL
+  LOGICAL, SAVE :: cosp_varsdefined = .FALSE. ! ug PAS THREADPRIVATE ET C'EST NORMAL
+  REAL, SAVE  :: Cosp_fill_value
+!$OMP THREADPRIVATE(Cosp_fill_value)
+ 
 
 CONTAINS
@@ -172 +175 @@
 
 #ifdef CPP_XIOS
+    
+! recuperer la valeur indefine Xios
+    CALL xios_get_field_attr("clcalipso",default_value=Cosp_fill_value)
     ! ug R\'eglage du calendrier xios
     !Temps julian => an, mois, jour, heure
@@ -230 +236 @@
 
       CALL histvert(cosp_nidfiles(iff),"column","column","count",Ncolumns,column_ax(1:Ncolumns),nvertcol(iff))
+
+!!! Valeur indefinie en cas IOIPSL
+     Cosp_fill_value=0.
+
 #endif
 

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

@@ -92 +92 @@
      do ip = 1,Npoints
      if(stlidar%lidarcld(ip,k).eq.R_UNDEF)then
-      stlidar%lidarcld(ip,k)=0.
+      stlidar%lidarcld(ip,k)=Cosp_fill_value
      endif
      enddo
@@ -99 +99 @@
       do ip = 1,Npoints
        if(stlidar%cfad_sr(ip,ii,k).eq.R_UNDEF)then
-        stlidar%cfad_sr(ip,ii,k)=0.
+        stlidar%cfad_sr(ip,ii,k)=Cosp_fill_value
        endif
       enddo
@@ -108 +108 @@
    do k = 1,Nlevlmdz
      if(sglidar%beta_mol(ip,k).eq.R_UNDEF)then
-      sglidar%beta_mol(ip,k)=0.
+      sglidar%beta_mol(ip,k)=Cosp_fill_value
      endif
 
      do ii= 1,Ncolumns
        if(sglidar%beta_tot(ip,ii,k).eq.R_UNDEF)then
-        sglidar%beta_tot(ip,ii,k)=0.
+        sglidar%beta_tot(ip,ii,k)=Cosp_fill_value
        endif
      enddo
@@ -123 +123 @@
     do ip = 1,Npoints
      if(stlidar%cldlayer(ip,k).eq.R_UNDEF)then
-      stlidar%cldlayer(ip,k)=0.
+      stlidar%cldlayer(ip,k)=Cosp_fill_value
      endif
     enddo
@@ -167 +167 @@
    do ip = 1,Npoints
     if(isccp%totalcldarea(ip).eq.-1.E+30)then
-      isccp%totalcldarea(ip)=0.
+      isccp%totalcldarea(ip)=Cosp_fill_value
     endif
     if(isccp%meanptop(ip).eq.-1.E+30)then
-      isccp%meanptop(ip)=0.
+      isccp%meanptop(ip)=Cosp_fill_value
     endif
     if(isccp%meantaucld(ip).eq.-1.E+30)then
-      isccp%meantaucld(ip)=0.
+      isccp%meantaucld(ip)=Cosp_fill_value
     endif
     if(isccp%meanalbedocld(ip).eq.-1.E+30)then
-      isccp%meanalbedocld(ip)=0.
+      isccp%meanalbedocld(ip)=Cosp_fill_value
     endif
     if(isccp%meantb(ip).eq.-1.E+30)then
-      isccp%meantb(ip)=0.
+      isccp%meantb(ip)=Cosp_fill_value
     endif
     if(isccp%meantbclr(ip).eq.-1.E+30)then
-      isccp%meantbclr(ip)=0.
+      isccp%meantbclr(ip)=Cosp_fill_value
     endif
 
@@ -188 +188 @@
      do ii=1,7
      if(isccp%fq_isccp(ip,ii,k).eq.-1.E+30)then
-      isccp%fq_isccp(ip,ii,k)=0.
+      isccp%fq_isccp(ip,ii,k)=Cosp_fill_value
      endif
      enddo
@@ -195 +195 @@
     do ii=1,Ncolumns
      if(isccp%boxtau(ip,ii).eq.-1.E+30)then
-       isccp%boxtau(ip,ii)=0.
+       isccp%boxtau(ip,ii)=Cosp_fill_value
      endif
     enddo
@@ -201 +201 @@
     do ii=1,Ncolumns
      if(isccp%boxptop(ip,ii).eq.-1.E+30)then
-       isccp%boxptop(ip,ii)=0.
+       isccp%boxptop(ip,ii)=Cosp_fill_value
      endif
     enddo

LMDZ5/trunk/libf/phylmd/radlwsw_m.F90

@@ -1000 +1000 @@
           zalbpla(i)   = ZFSUP(i,klev+1)/ZFSDN(i,klev+1)
          ENDIF
-         zsollwdown(i)= ZFLDN(i,1)
+!!! 5 juin 2015
+!!! Correction MP bug RRTM
+         zsollwdown(i)= -1.*ZFLDN(i,1)
       ENDDO
 !     print*,'OK2'

LMDZ5/trunk/libf/phylmd/rrtm/readaerosol_optic_rrtm.F90

@@ -303 +303 @@
   END DO
 
+! AI 12 juin 2015
+! Modif Olivier pour prendre en compte des prop optiq aero ancienne pour RRTM 2bds 
+
+IF (NSW.EQ.2) THEN
+!--old aerosol properties
+
+   ! old aeropt_2bands for rrtm
+   CALL aeropt_2bands( &
+         pdel, m_allaer, pdtphys, rhcl, &
+         tau_aero, piz_aero, cg_aero,   &
+         m_allaer_pi, flag_aerosol, &
+         pplay, t_seri, presnivs)
+
+   ! aeropt_5wv only for validation and diagnostics.
+   CALL aeropt_5wv(                    &
+        pdel, m_allaer,                &
+        pdtphys, rhcl, aerindex,       &
+        flag_aerosol, pplay, t_seri,   &
+        tausum_aero, tau3d_aero, presnivs)
+
+ELSEIF (NSW.EQ.6) THEN
+!--new aerosol propertie
   ! aeropt_6bands for rrtm
   CALL aeropt_6bands_rrtm( &
@@ -316 +338 @@
        flag_aerosol, pplay, t_seri,   &
        tausum_aero, tau3d_aero )
+
+ELSE
+
+    PRINT *,'Cas NSW non prevu pour RRTM - NSW=',NSW
+    STOP
+
+ENDIF
 
   ! Diagnostics calculation for CMIP5 protocol

LMDZ5/trunk/libf/phylmd/rrtm/rrtm_rtrn1a_140gp.F90

@@ -496 +496 @@
 !    Z_FACCMB2D(I_LEV-1) = Z_FACCLD1D(I_LEV-1) * Z_FACCLR2D(I_LEV) *&
 !     & (1.0_JPRB - Z_CLDFRAC(I_LEV+1))  
-    if (istcldd(i_lev).ne.1) then
+    if (istcldd(i_lev).ne.1.and.i_lev.ne.0) then
        z_faccmb1d(i_lev-1) = max(0.,min(z_cldfrac(i_lev+1)-z_cldfrac(i_lev), &
                             z_cldfrac(i_lev-1)-z_cldfrac(i_lev)))