Changeset 1319 for LMDZ4/trunk/libf/dyn3d/ce0l.F90

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

File:

• LMDZ4/trunk/libf/dyn3d/ce0l.F90 (moved) (moved from LMDZ4/trunk/libf/dyn3d/ce0l.F) (1 diff)

LMDZ4/trunk/libf/dyn3d/ce0l.F90

@@ -2 +2 @@
 ! $Id$
 !
-       PROGRAM create_etat0_limit
+!-------------------------------------------------------------------------------
+!
+PROGRAM ce0l
+!
+!-------------------------------------------------------------------------------
+! Purpose: Calls etat0, creates initial states and limit_netcdf
+!
+!     interbar=.T. for barycentric interpolation inter_barxy
+!     extrap  =.T. for data extrapolation, like for the SSTs when file does not
+!                  contain ocean points only.
+!     oldice  =.T. for old-style ice, obtained using grille_m (grid_atob).
+!     masque is created in etat0, passed to limit to ensure consistancy.
+!-------------------------------------------------------------------------------
 #ifdef CPP_EARTH
 ! This prog. is designed to work for Earth
-       USE dimphy
-       USE comgeomphy
-       USE infotrac
+  USE dimphy
+  USE comgeomphy
+  USE infotrac
+
 #ifdef CPP_IOIPSL
-       use ioipsl, only: ioconf_calendar
+  USE ioipsl, ONLY: ioconf_calendar
 #endif
-       IMPLICIT NONE
-c
-c
-c     Programme d'appel a etat0, creation des etats initiaux et limit_netcdf
-c   
-c
-c     interbar = .T . si appel a  interpol. barycentrique inter_barxy
-c
-c     extrap   = .T . si on fait une extrapolation de donnees , comme pour
-c       les  SST  lorsque  le fichier ne contient pas uniquement  des points 
-c     oceaniques .
-c
-c     oldice   = .T. si l'on veut garder les anciennes glaces , obtenues
-c     par  grille_m  ( grid_atob ) .
-c
-c     on cree le masque dans etat0 que l'on passe ensuite dans limit pour 
-c     garder les coherences
 
-      LOGICAL interbar, extrap , oldice
-      PARAMETER ( interbar = .true. , extrap = .FALSE. , oldice=.false.)
+#endif
+  IMPLICIT NONE
+#ifndef CPP_EARTH
+  WRITE(lunout,*)'limit_netcdf: Earth-specific routine, needs Earth physics'
+#else
+!-------------------------------------------------------------------------------
+! Local variables:
+  LOGICAL, PARAMETER :: interbar=.TRUE., extrap=.FALSE., oldice=.FALSE.
 #include "dimensions.h"
 #include "paramet.h"
 #include "indicesol.h"
-#include  "control.h"
-      REAL :: masque(iip1,jjp1)
-!      REAL :: pctsrf(iim*(jjm-1)+2, nbsrf)
+#include "iniprint.h"
+#include "control.h"
+#include "temps.h"
+#include "logic.h"
+  INTEGER, PARAMETER            :: longcles=20
+  REAL,    DIMENSION(longcles)  :: clesphy0
+  REAL,    DIMENSION(iip1,jjp1) :: masque
+  CHARACTER(LEN=15)             :: calnd
+!-------------------------------------------------------------------------------
+  CALL conf_gcm( 99, .TRUE. , clesphy0 )
 
-      IF (config_inca /= 'none') THEN
-#ifdef INCA
-         call init_const_lmdz(
-     $        nbtr,anneeref,dayref,
-     $        iphysiq, day_step,nday)
-#endif
-         print *, 'nbtr =' , nbtr 
-      END IF
-
-      CALL Init_Phys_lmdz(iim,jjp1,llm,1,(/(jjm-1)*iim+2/))
-      PRINT *,'---> klon=',klon
-      call InitComgeomphy
+  CALL Init_Phys_lmdz(iim,jjp1,llm,1,(/(jjm-1)*iim+2/))
+  WRITE(lunout,*)'---> klon=',klon
+  CALL InitComgeomphy
 
 #ifdef CPP_IOIPSL
-      call ioconf_calendar('360d')
+  SELECT CASE(calend)
+    CASE('earth_360d');CALL ioconf_calendar('360d');      calnd='a 360 jours/an'
+    CASE('earth_365d');CALL ioconf_calendar('noleap');    calnd='a 365 jours/an'
+    CASE('earth_366d');CALL ioconf_calendar('366d');      calnd='bissextile'
+    CASE('gregorian'); CALL ioconf_calendar('gregorian'); calnd='gregorien'
+    CASE('standard');  CALL ioconf_calendar('gregorian'); calnd='gregorien'
+    CASE('julian');    CALL ioconf_calendar('julian');    calnd='julien'
+    CASE('proleptic_gregorian'); CALL ioconf_calendar('gregorian')
+  !--- DC Bof...  => IOIPSL a mettre a jour: proleptic_gregorian /= gregorian
+    CASE DEFAULT
+      CALL abort_gcm('ce0l','Mauvais choix de calendrier',1)
+  END SELECT
+  WRITE(lunout,*)'CALENDRIER CHOISI: Terrestre '//TRIM(calnd)
 #endif
 
-      WRITE(6,*) '  *********************  '
-      WRITE(6,*) ' interbar = ',interbar
-      CALL etat0_netcdf ( interbar, masque )
-c
-      WRITE(6,1)
-      WRITE(6,*) '  *********************  '
-      WRITE(6,*) '  ***  Limit_netcdf ***  '
-      WRITE(6,*) '  *********************  '
-      WRITE(6,1)
-      
-c     
-      CALL  limit_netcdf ( interbar, extrap , oldice, masque)
+  IF (config_inca /= 'none') THEN
+#ifdef INCA
+    CALL init_const_lmdz(nbtr,anneeref,dayref,iphysiq,day_step,nday)
+    CALL init_inca_para(iim,jjm+1,klon,1,klon_mpi_para_nb,0)
+    WRITE(lunout,*)'nbtr =' , nbtr 
+#endif
+  END IF
 
-1     FORMAT(//)
+  WRITE(lunout,'(//)')
+  WRITE(lunout,*) '  *********************  '
+  WRITE(lunout,*) '  ***  etat0_netcdf ***  '
+  WRITE(lunout,*) '  *********************  '
+  WRITE(lunout,'(//)')
+  WRITE(lunout,*) ' interbar = ',interbar
+  CALL etat0_netcdf(interbar,masque,ok_etat0)
+
+  IF(ok_limit) THEN
+  WRITE(lunout,'(//)')
+  WRITE(lunout,*) '  *********************  '
+  WRITE(lunout,*) '  ***  Limit_netcdf ***  '
+  WRITE(lunout,*) '  *********************  '
+  WRITE(lunout,'(//)')
+  CALL limit_netcdf(interbar,extrap,oldice,masque)
+  END IF
 
 #endif
-! of #ifdef CPP_EARTH
-      STOP
-      END
+! of #ifndef CPP_EARTH #else
+  STOP
+
+END PROGRAM ce0l
+!
+!-------------------------------------------------------------------------------