Changeset 1319 for LMDZ4/trunk/libf/phylmd/readaerosol.F90

Timestamp:

Feb 23, 2010, 10:29:54 PM (14 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

File:

• LMDZ4/trunk/libf/phylmd/readaerosol.F90 (modified) (2 diffs)

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.