Changeset 1319 for LMDZ4/trunk/libf/phylmd

Timestamp:

Feb 23, 2010, 10:29:54 PM (15 years ago)

Author:

Laurent Fairhead

Message:

• Modifications to the start and limit creation routines to account for different

calendars

• Modification to phyetat0 to force the mask read in the start files to match the

surface fractions read in the limit file

• Force readaerosol.F90 to read in aerosols file with 12 timesteps

---

• Modifications aux routines de créations des fichiers start et limit pour prendre

en compte différents calendriers

• Modification à phyetat0 pour forcer le masque lu dans le fichier start à être

compatible avec les fractions de surface lu dans le fichier limit

• Forcer readaerosol à ne lire que des fichiers à 12 pas de temps

Location:

LMDZ4/trunk/libf/phylmd

Files:

• phyetat0.F (modified) (2 diffs)
• physiq.F (modified) (1 diff)
• readaerosol.F90 (modified) (2 diffs)

LMDZ4/trunk/libf/phylmd/phyetat0.F

@@ -228 +228 @@
      $          'coherente ', i, zmasq(i), pctsrf(i, is_ter)
      $          ,pctsrf(i, is_lic)
+            zmasq(i) = fractint(i)
         ENDIF 
       END DO 
@@ -237 +238 @@
      $          'coherente ', i, zmasq(i) , pctsrf(i, is_oce)
      $          ,pctsrf(i, is_sic)
+            zmasq(i) = fractint(i)
         ENDIF 
       END DO 

LMDZ4/trunk/libf/phylmd/physiq.F

@@ -2548 +2548 @@
 
 !   les ratqs sont une combinaison de ratqss et ratqsc
-       print*,'PHYLMD NOUVEAU TAU_RATQS ',tau_ratqs
+!       print*,'PHYLMD NOUVEAU TAU_RATQS ',tau_ratqs
 
          if (tau_ratqs>1.e-10) then

LMDZ4/trunk/libf/phylmd/readaerosol.F90

@@ -200 +200 @@
     REAL, DIMENSION(klon,12), INTENT(OUT) :: psurf_out    ! Surface pression for 12 months
     REAL, DIMENSION(klon,12), INTENT(OUT) :: load_out     ! Aerosol mass load in each column
+    INTEGER                               :: nbr_tsteps   ! number of month in file read
 
 ! Local variables
@@ -273 +274 @@
           CALL abort_gcm('get_aero_fromfile', 'latitudes do not correspond between file and model',1)
        END IF
+
+! 1.5) Check number of month in file opened
+!
+!**************************************************************************************************
+       ierr = nf90_inq_dimid(ncid, 'TIME',dimid)
+       CALL check_err( nf90_inquire_dimension(ncid, dimid, len = nbr_tsteps) )
+!       IF (nbr_tsteps /= 12 .AND. nbr_tsteps /= 14) THEN
+       IF (nbr_tsteps /= 12 ) THEN
+         CALL abort_gcm('get_aero_fromfile', 'not the right number of months in aerosol file read (should be 12 for the moment)',1)
+       ENDIF
+
 
 ! 2) Check if old or new file is avalabale.