Changeset 2297 for LMDZ5/trunk/libf/phylmd/cosp

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/cosp

Files:

• cosp_output_mod.F90 (modified) (3 diffs)
• cosp_output_write_mod.F90 (modified) (8 diffs)

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