Changeset 1912 for LMDZ5/trunk


Ignore:
Timestamp:
Dec 5, 2013, 6:32:35 PM (11 years ago)
Author:
musat
Message:

1) Modifications pour faire des simulations par an avec un calendrier realiste (365 jours ou autre).

Il faut mettre une frequence de sortie de -1 (variable phys_out_filetimesteps dans config.def) pour
que IOIPSL calcule les moyennes mensuels en prenant en compte des longuers variables de chaque
mois. Par exemple, pour le fichier histmth (1er fichier) et histmthNMC (7eme).
phys_out_filetimesteps= -1 1day 6hr 6hr 6hr 1d -1 1day 6hr

2) Corrections titres variables niveaux de pression des fichiers histmth, histday, etc

Location:
LMDZ5/trunk/libf/phylmd
Files:
8 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/phylmd/calcul_STDlev.h

    r1907 r1912  
    66        missing_val=nf90_fill_real
    77!
    8 !IM freq_moyNMC = frequences auxquelles on moyenne les champs accumules
    9 !IM               sur les niveaux de pression standard du NMC
    10       DO n=1, nout
    11        freq_moyNMC(n)=freq_outNMC(n)/freq_calNMC(n)
    12       ENDDO
    13 !
    14         CALL ini_undefSTD(itap,freq_outNMC)
     8        CALL ini_undefSTD(itap,itapm1)
    159!
    1610!IM on interpole les champs sur les niveaux STD de pression
    … …  
    135129!IM on somme les valeurs toutes les freq_calNMC secondes
    136130!
    137        CALL undefSTD(itap,freq_calNMC, read_climoz)
    138 !
    139 !IM on moyenne a la fin du mois ou du jour (toutes les freq_outNMC secondes)
    140 !
    141        CALL moy_undefSTD(itap,freq_outNMC,freq_moyNMC)
     131       CALL undefSTD(itap, read_climoz)
     132!
     133!IM on moyenne a la fin du mois, du jour ou toutes les 6h
     134!
     135       CALL moy_undefSTD(itap,itapm1)
    142136!
    143137       CALL plevel(klon,klev,.true.,pplay,50000., &

  • LMDZ5/trunk/libf/phylmd/calcul_divers.h

    r1907 r1912  
    22! $Header$
    33!
    4 
    5 !     Initialisations diverses au "debut" du mois
    6       IF(debut) THEN
    7          nday_rain(:)=0.
    8 
     4! Initialisations diverses au tout debut
     5      IF(itap.EQ.1) THEN
     6         itapm1=0
    97!        surface terre
    10          paire_ter(:)=0.
    118         DO i=1, klon
    129            IF(pctsrf(i,is_ter).GT.0.) THEN
    … …  
    1613      ENDIF
    1714
    18 !IM   Calcul une fois par jour : total_rain, nday_rain
    19       IF(MOD(itap,INT(un_jour/dtime)).EQ.0) THEN
     15! Initialisation debut de mois
     16      IF(itap.EQ.itapm1+1) THEN
     17        nday_rain(:)=0.
     18!       print*,'initialisation mois suivants day_rain itap',itap
     19      ENDIF
     20
     21! Calcul fin de journee : total_rain, nday_rain
     22      IF(MOD(itap,NINT(un_jour/dtime)).EQ.0) THEN
     23!        print*,'calcul nday_rain itap ',itap
    2024         DO i = 1, klon
    2125            total_rain(i)=rain_fall(i)+snow_fall(i) 
    … …  
    2327         ENDDO
    2428      ENDIF
     29
     30! Initialisation fin de mois
     31      IF(MOD(itap-itapm1,NINT(mth_len*un_jour/dtime)).EQ.0) THEN
     32        itapm1=itapm1+NINT(mth_len*un_jour/dtime)
     33!       print*,'initialisation itapm1 ',itapm1
     34      ENDIF

  • LMDZ5/trunk/libf/phylmd/ini_undefSTD.F

    r1907 r1912  
    33!
    44
    5       SUBROUTINE ini_undefSTD(itap,
    6      $           freq_outNMC)
     5      SUBROUTINE ini_undefSTD(itap,itapm1)
    76      USE dimphy
    87      USE phys_state_var_mod ! Variables sauvegardees de la physique
     8      USE phys_cal_mod, only : mth_len
    99      IMPLICIT none
     10      include "clesphys.h"
    1011c
    1112c====================================================================
    … …  
    3637c variables Input/Output
    3738c     INTEGER nlevSTD, klevSTD, itap
    38       INTEGER itap
     39      INTEGER itap, itapm1
    3940c     PARAMETER(klevSTD=17)
    4041c     REAL dtime
    … …  
    4445      INTEGER i, k, n
    4546c     PARAMETER(nout=3) !nout=1 day/nout=2 mth/nout=3 NMC
    46       REAL freq_outNMC(nout)
     47c     REAL freq_outNMC(nout)
     48      REAL un_jour
     49      PARAMETER(un_jour=86400.)
    4750c
    4851c variables Output
    … …  
    5255      DO n=1, nout
    5356c
    54 c initialisation variables en debut de la journee ou du mois
     57c initialisation variables en debut du mois, de la journee ou des 6h
    5558c
    56        IF(MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.1.) THEN
     59       IF(n.EQ.1.AND.itap-itapm1.EQ.1.OR.
     60     &n.GT.1.AND.MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.1) THEN
     61c       print*,'ini_undefSTD n itap',n,itap
    5762        DO k=1, nlevSTD
    5863         DO i=1, klon
    … …  
    7883         ENDDO !i
    7984        ENDDO !k
    80 c
    81        ENDIF !MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.1.
    82 c
     85       ENDIF !
    8386      ENDDO !n
    84 c
    8587      RETURN
    8688      END 

  • LMDZ5/trunk/libf/phylmd/moy_undefSTD.F

    r1907 r1912  
    22! $Id$
    33!
    4       SUBROUTINE moy_undefSTD(itap,freq_outNMC,freq_moyNMC)
     4      SUBROUTINE moy_undefSTD(itap,itapm1)
    55      USE netcdf
    66      USE dimphy
    7       USE phys_state_var_mod ! Variables sauvegardees de la physique
     7      USE phys_state_var_mod
     8      USE phys_cal_mod, only : mth_len
    89      IMPLICIT none
     10      include "clesphys.h"
    911c
    1012c====================================================================
    … …  
    3335c     INTEGER nlevSTD, klevSTD, itap
    3436c     PARAMETER(klevSTD=17)
    35       INTEGER itap
     37      INTEGER itap, itapm1
    3638c
    3739c variables locales
    … …  
    4042      INTEGER i, k, n
    4143c     REAL dtime, freq_outNMC(nout), freq_moyNMC(nout)
    42       REAL freq_outNMC(nout), freq_moyNMC(nout)
     44c     REAL freq_outNMC(nout), freq_calNMC(nout)
     45      REAL freq_moyNMC(nout)
    4346c
    4447c variables Output
    … …  
    4649c     REAL tsumSTD(klon,klevSTD,nout)
    4750c
     51      REAL un_jour
     52      PARAMETER(un_jour=86400.)
    4853      REAL missing_val
    4954c
    … …  
    5156c
    5257      DO n=1, nout
     58       IF(freq_outNMC(n).LT.0) THEN
     59         freq_moyNMC(n)=(mth_len*un_jour)/freq_calNMC(n)
     60c        print*,'moy_undefSTD n freq_out freq_moy =',
     61c    $n,freq_moyNMC(n)
     62       ENDIF         
    5363c
    54 c calcul 1 fois par jour
     64c calcul 1 fois pas mois, 1 fois par jour ou toutes les 6h
    5565c
    56        IF(MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.0) THEN
     66       IF(n.EQ.1.AND.itap.EQ.itapm1.OR.
     67     $n.GT.1.AND.MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.0) THEN
     68c
     69c       print*,'moy_undefSTD n itap itapm1',n,itap,itapm1
    5770c
    5871        DO k=1, nlevSTD
    … …  
    6174           tsumSTD(i,k,n)=tsumSTD(i,k,n)/
    6275     $     (freq_moyNMC(n)-tnondef(i,k,n))
    63 cIM BEG
    6476          usumSTD(i,k,n)=usumSTD(i,k,n)/
    6577     $     (freq_moyNMC(n)-tnondef(i,k,n))
    … …  
    96108          O3daysumSTD(i,k,n)=O3daysumSTD(i,k,n)/
    97109     $     (freq_moyNMC(n)-tnondef(i,k,n))
    98 cIM END
    99110          ELSE
    100111           tsumSTD(i,k,n)=missing_val

  • LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

    r1907 r1912  
    500500
    501501  TYPE(ctrl_out), SAVE, DIMENSION(7) :: o_uSTDlevs     = (/                    &
    502       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u850', "Zonal wind 1hPa", "m/s",     &
    503       (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    504       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u700', "Zonal wind 2hPa", "m/s",     &
    505       (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    506       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u500', "Zonal wind 3hPa", "m/s",     &
    507       (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    508       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u200', "Zonal wind 4hPa", "m/s",     &
    509       (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    510       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u100', "Zonal wind 5hPa", "m/s",     &
    511       (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    512       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u50', "Zonal wind 6hPa", "m/s",     &
    513       (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    514       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u10', "Zonal wind 7hPa", "m/s",     &
     502      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u850', "Zonal wind 850hPa", "m/s",     &
     503      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
     504      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u700', "Zonal wind 700hPa", "m/s",     &
     505      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
     506      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u500', "Zonal wind 500hPa", "m/s",     &
     507      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
     508      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u200', "Zonal wind 200hPa", "m/s",     &
     509      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
     510      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u100', "Zonal wind 100hPa", "m/s",     &
     511      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
     512      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u50', "Zonal wind 50hPa", "m/s",     &
     513      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
     514      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'u10', "Zonal wind 10hPa", "m/s",     &
    515515      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)) /)
    516516
    517517  TYPE(ctrl_out), SAVE, DIMENSION(7) :: o_vSTDlevs     = (/                     &
    518       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v850', "Meridional wind 1hPa", "m/s", &
     518      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v850', "Meridional wind 850hPa", "m/s", &
    519519      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)),  &
    520       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v700', "Meridional wind 2hPa", "m/s", &
     520      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v700', "Meridional wind 700hPa", "m/s", &
    521521      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)),  &
    522       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v500', "Meridional wind 3hPa", "m/s", &
     522      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v500', "Meridional wind 500hPa", "m/s", &
    523523      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)),  &
    524       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v200', "Meridional wind 4hPa", "m/s", &
     524      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v200', "Meridional wind 200hPa", "m/s", &
    525525      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)),  &
    526       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v100', "Meridional wind 5hPa", "m/s", &
     526      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v100', "Meridional wind 100hPa", "m/s", &
    527527      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)),  &
    528       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v50', "Meridional wind 6hPa", "m/s",  &
     528      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v50', "Meridional wind 50hPa", "m/s",  &
    529529      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)),  &
    530       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v10', "Meridional wind 7hPa", "m/s",  &
     530      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'v10', "Meridional wind 10hPa", "m/s",  &
    531531      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)) /)
    532532
    … …  
    534534      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w850', "Vertical wind 1hPa", "Pa/s", &
    535535      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    536       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w700', "Vertical wind 2hPa", "Pa/s", &
    537       (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    538       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w500', "Vertical wind 3hPa", "Pa/s", &
    539       (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    540       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w200', "Vertical wind 4hPa", "Pa/s", &
    541       (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    542       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w100', "Vertical wind 5hPa", "Pa/s", &
    543       (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    544       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w50', "Vertical wind 6hPa", "Pa/s",  &
    545       (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    546       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w10', "Vertical wind 7hPa", "Pa/s",  &
     536      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w700', "Vertical wind 700hPa", "Pa/s", &
     537      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
     538      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w500', "Vertical wind 500hPa", "Pa/s", &
     539      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
     540      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w200', "Vertical wind 200hPa", "Pa/s", &
     541      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
     542      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w100', "Vertical wind 100hPa", "Pa/s", &
     543      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
     544      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w50', "Vertical wind 50hPa", "Pa/s",  &
     545      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
     546      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'w10', "Vertical wind 10hPa", "Pa/s",  &
    547547      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)) /)
    548548
    … …  
    550550      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t850', "Temperature 1hPa", "K",      &
    551551      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    552       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t700', "Temperature 2hPa", "K",      &
    553       (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    554       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t500', "Temperature 3hPa", "K",      &
    555       (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    556       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t200', "Temperature 4hPa", "K",      &
    557       (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    558       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t100', "Temperature 5hPa", "K",      &
    559       (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    560       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t50',  "Temperature 6hPa", "K",      &
    561       (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    562       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t10',  "Temperature 7hPa", "K",      &
     552      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t700', "Temperature 700hPa", "K",      &
     553      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
     554      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t500', "Temperature 500hPa", "K",      &
     555      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
     556      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t200', "Temperature 200hPa", "K",      &
     557      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
     558      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t100', "Temperature 100hPa", "K",      &
     559      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
     560      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t50',  "Temperature 50hPa", "K",      &
     561      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
     562      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'t10',  "Temperature 10hPa", "K",      &
    563563      (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)) /)
    564564
    … …  
    566566      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q850', "Specific humidity 1hPa", &
    567567      "kg/kg", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    568       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q700', "Specific humidity 2hPa", &
     568      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q700', "Specific humidity 700hPa", &
    569569      "kg/kg", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    570       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q500', "Specific humidity 3hPa", &
     570      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q500', "Specific humidity 500hPa", &
    571571      "kg/kg", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    572       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q200', "Specific humidity 4hPa", &
     572      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q200', "Specific humidity 200hPa", &
    573573      "kg/kg", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    574       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q100', "Specific humidity 5hPa", &
     574      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q100', "Specific humidity 100hPa", &
    575575      "kg/kg", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    576       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q50', "Specific humidity 6hPa",  &
     576      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q50', "Specific humidity 50hPa",  &
    577577      "kg/kg", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    578       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q10', "Specific humidity 7hPa", &
     578      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'q10', "Specific humidity 10hPa", &
    579579      "kg/kg", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)) /)
    580580
    … …  
    582582      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z850', "Geopotential height 1hPa",        &
    583583      "m", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    584       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z700', "Geopotential height 2hPa",        &
     584      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z700', "Geopotential height 700hPa",        &
    585585      "m", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    586       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z500', "Geopotential height 3hPa",        &
     586      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z500', "Geopotential height 500hPa",        &
    587587      "m", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    588       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z200', "Geopotential height 4hPa",        &
     588      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z200', "Geopotential height 200hPa",        &
    589589      "m", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    590       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z100', "Geopotential height 5hPa",        &
     590      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z100', "Geopotential height 100hPa",        &
    591591      "m", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    592       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z50', "Geopotential height 6hPa",         &
     592      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z50', "Geopotential height 50hPa",         &
    593593      "m", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)), &
    594       ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z10', "Geopotential height 7hPa",         &
     594      ctrl_out((/ 1, 7, 7, 10, 10, 10, 11, 11, 11 /),'z10', "Geopotential height 10hPa",         &
    595595      "m", (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)) /)
    596596

  • LMDZ5/trunk/libf/phylmd/phys_output_mod.F90

    r1907 r1912  
    4444    USE surface_data, ONLY : ok_snow
    4545    USE phys_output_ctrlout_mod
     46    USE mod_grid_phy_lmdz, only: klon_glo
    4647
    4748#ifdef CPP_XIOS
    … …  
    207208    ecrit_files(5) = ecrit_LES
    208209    ecrit_files(6) = ecrit_ins
    209     ecrit_files(7) = freq_outNMC(1)/86400.
    210     ecrit_files(8) = freq_outNMC(2)/86400.
    211     ecrit_files(9) = freq_outNMC(3)/86400.
     210    ecrit_files(7) = freq_outNMC(1)
     211    ecrit_files(8) = freq_outNMC(2)
     212    ecrit_files(9) = freq_outNMC(3)
    212213
    213214    !! Lectures des parametres de sorties dans physiq.def
    … …  
    248249    ! ug Réglage du calendrier xios
    249250    !Temps julian => an, mois, jour, heure
    250     CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
     251    IF (klon_glo==1) THEN                                       
     252        CALL ymds2ju(annee_ref, 1, day_ref, hour, zjulian)
     253    ELSE
     254        CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
     255    END IF
    251256    CALL ju2ymds(zjulian, x_an, x_mois, x_jour, x_heure)
    252257    CALL wxios_set_cal(dtime, calend, x_an, x_mois, x_jour, x_heure)
    … …  
    278283          ! ...)*86400.
    279284          ecrit_files(iff)=ecrit_files(iff)*86400.
     285      ELSE IF (chtimestep(iff).eq.'-1') then
     286          PRINT*,'ecrit_files(',iff,') < 0 so IOIPSL work on different'
     287          PRINT*,'months length'
     288          ecrit_files(iff)=-1.
    280289      else
    281290       CALL convers_timesteps(chtimestep(iff),dtime,ecrit_files(iff))
    … …  
    304313
    305314          idayref = day_ref
    306 !          CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)       
    307 ! correction pour l heure initiale                               !jyg
    308 !                                                                !jyg
    309           CALL ymds2ju(annee_ref, 1, idayref, hour, zjulian)         !jyg
    310 ! correction pour l heure initiale                               !jyg
    311 !                                                                !jyg
    312 !!!      CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)       !jyg
    313 ! correction pour l heure initiale                               !jyg
    314 !                                                                !jyg
    315 !      CALL ymds2ju(annee_ref, 1, idayref, hour, zjulian)         !jyg
    316 
     315! A noter pour
     316! l heure initiale - dans les fichiers histoire hist* - on met comme
     317! heure de debut soit la vraie heure (pour le 1D) soit 0h (pour le 3D)
     318! afin d avoir une seule sortie mensuelle par mois lorsque l on tourne
     319! par annee (IM).
     320!
     321     IF (klon_glo==1) THEN                                     
     322         CALL ymds2ju(annee_ref, 1, idayref, hour, zjulian)     !jyg
     323     ELSE
     324         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
     325     END IF                                                     
    317326!!!!!!!!!!!!!!!!! Traitement dans le cas ou l'on veut stocker sur un domaine limite !!
    318327!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    … …  
    479488    ecrit_LES = ecrit_files(5)
    480489    ecrit_ins = ecrit_files(6)
    481 
     490    freq_outNMC(1) = ecrit_files(7)
     491    freq_outNMC(2) = ecrit_files(8)
     492    freq_outNMC(3) = ecrit_files(9)
     493
     494    WRITE(lunout,*)'freq_outNMC(1)=',freq_outNMC(1)
     495    WRITE(lunout,*)'freq_outNMC(2)=',freq_outNMC(2)
     496    WRITE(lunout,*)'freq_outNMC(3)=',freq_outNMC(3)
    482497    WRITE(lunout,*)'swaero_diag=',swaero_diag
    483498    WRITE(lunout,*)'Fin phys_output_mod.F90'

  • LMDZ5/trunk/libf/phylmd/physiq.F90

    r1907 r1912  
    493493  REAL un_jour
    494494  PARAMETER(un_jour=86400.)
     495  INTEGER itapm1
     496  SAVE itapm1
    495497  !======================================================================
    496498  !
    … …  
    34313433  !   SORTIES
    34323434  !=======================================================================
    3433 
     3435  !
     3436  !IM initialisation + calculs divers diag AMIP2
     3437  !
     3438  include "calcul_divers.h"
     3439  !
    34343440  !IM Interpolation sur les niveaux de pression du NMC
    34353441  !   -------------------------------------------------
    … …  
    34483454     ENDDO
    34493455  ENDDO
    3450   !
    3451   !IM initialisation + calculs divers diag AMIP2
    3452   !
    3453   include "calcul_divers.h"
    34543456  !
    34553457  IF (type_trac == 'inca') THEN

  • LMDZ5/trunk/libf/phylmd/undefSTD.F

    r1907 r1912  
    22! $Id$
    33!
    4       SUBROUTINE undefSTD(itap,freq_calNMC, read_climoz)
     4      SUBROUTINE undefSTD(itap, read_climoz)
    55      USE netcdf
    66      USE dimphy
    77      USE phys_state_var_mod ! Variables sauvegardees de la physique
    88      IMPLICIT none
     9      include "clesphys.h"
    910c
    1011c====================================================================
    … …  
    4243c     PARAMETER(nout=3) !nout=1 : day; =2 : mth; =3 : NMC
    4344      INTEGER i, k, n
    44       REAL freq_calNMC(nout)
     45c     REAL freq_calNMC(nout)
    4546      INTEGER read_climoz
    4647c
Note: See TracChangeset for help on using the changeset viewer.