Ignore:
Timestamp:
Oct 6, 2006, 1:06:22 PM (18 years ago)
Author:
lmdzadmin
Message:

On passe a des ecrit_ins, ecrit_day, etc en nombre de jours (REAL)
On lit frequence ecriture traceurs ecrit_trac dans physiq.def
Correction petits pbs ini_histrac.h
IM

Location:
LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd
Files:
15 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/clesphys.h

    r652 r723  
    2626       LOGICAL ok_isccp, ok_regdyn
    2727       REAL lonmin_ins, lonmax_ins, latmin_ins, latmax_ins
    28        INTEGER ecrit_ins, ecrit_hf, ecrit_hf2mth, ecrit_day
    29        INTEGER ecrit_mth, ecrit_tra, ecrit_reg
     28       REAL ecrit_ins, ecrit_hf, ecrit_hf2mth, ecrit_day
     29       REAL ecrit_mth, ecrit_tra, ecrit_reg
    3030
    3131       COMMON/clesphys/cycle_diurne, soil_model, new_oliq,

  • LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_bilKP_ave.h

    r644 r723  
    55c
    66         zsto = dtime
    7          zout = dtime * FLOAT(ecrit_day)
     7         zout = dtime * ecrit_day
    88         typeval=tave
    99c

  • LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histISCCP.h

    r644 r723  
    4747     .                 lmaxm1, cldtopres, nvert,'down')
    4848c
    49         IF(type_run.EQ."ENSP".OR.type_run.EQ."CLIM") THEN
     49      IF(type_run.EQ."ENSP".OR.type_run.EQ."CLIM") THEN
    5050c
    5151c variables a ecrire
    … …  
    6363     .                "ave(X)", zstophy,zout)
    6464c
    65         ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
     65      ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
    6666c
    6767         DO k=1, kmaxm1
    … …  
    8181     .                "ave(X)", zstophy,zout)
    8282c
    83         ENDIF
     83      ENDIF
    8484c
    8585        CALL histend(nid_isccp)

  • LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histday.h

    r644 r723  
    55c
    66       zstophy = dtime
    7        zstoday= dtime * FLOAT(ecrit_day)
    8        zout = dtime * FLOAT(ecrit_day)
     7       zstoday= dtime * ecrit_day
     8       zout = dtime * ecrit_day
    99c
    1010         idayref = day_ref

  • LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histday_seri.h

    r644 r723  
    22c $Header$
    33c
    4       IF (type_run.EQ."AMIP") THEN
     4       IF (type_run.EQ."AMIP") THEN
    55c
    66       zstophy = dtime
    7        zout = dtime * FLOAT(ecrit_day)
     7       zout = dtime * ecrit_day
    88c
    99         idayref = day_ref

  • LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histhf.h

    r644 r723  
    55      IF (ok_hf) THEN
    66c
    7        zstohf = dtime * FLOAT(ecrit_hf)
     7       zstohf = dtime * ecrit_hf
    88       zstophy = dtime
    9        zstoday = dtime * FLOAT(ecrit_day)
    10        zout = dtime * FLOAT(ecrit_hf)
     9       zstoday = dtime * ecrit_day
     10       zout = dtime * ecrit_hf
    1111c
    1212         idayref = day_ref
    … …  
    190190      ENDIF !lev_histhf.GE.4
    191191c
     192c
    192193      ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
    193194c

  • LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histhf3d.h

    r644 r723  
    88cIM     zstohf = dtime * FLOAT(ecrit_hf)
    99cIM     zout = dtime * FLOAT(ecrit_hf)
    10         zstohf = dtime * FLOAT(ecrit_hf)
    11         zout = dtime * FLOAT(ecrit_hf)
     10        zstohf = dtime * ecrit_hf
     11        zout = dtime * ecrit_hf
    1212c
    1313c       PRINT*, 'La frequence de sortie hf3d est de ', ecrit_hf

  • LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histmth.h

    r653 r723  
    55c
    66       zstophy = dtime
    7        zstomth = dtime * FLOAT(ecrit_mth) !valeurs mensuelles "instantannees"
     7       zstomth = dtime * ecrit_mth !valeurs mensuelles "instantannees"
    88       zstorad = dtime * FLOAT(radpas)    !variables stockees toutes les 2h
    9        zstohf = dtime *FLOAT(ecrit_hf)   !variables stockees toutes les 6h
    10        zout = dtime * FLOAT(ecrit_mth)    !frequence des sorties : 1 val. par mois
     9       zstohf = dtime * ecrit_hf   !variables stockees toutes les 6h
     10       zout = dtime * ecrit_mth    !frequence des sorties : 1 val. par mois
    1111c
    1212cIM    PRINT*,' La zstophy,zstorad,zstohf,zstomth,zout=',

  • LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histmthNMC.h

    r644 r723  
    55c
    66       zstophy = dtime
    7        zstohf = dtime * FLOAT(ecrit_hf)
    8        zstomth = dtime * FLOAT(ecrit_mth)
    9        zout = dtime * FLOAT(ecrit_mth)
     7       zstohf = dtime * ecrit_hf
     8       zstomth = dtime * ecrit_mth
     9       zout = dtime * ecrit_mth
    1010c
    1111         idayref = day_ref

  • LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histrac.h

    r646 r723  
    3939
    4040         zsto = pdtphys
    41          zout = pdtphys * FLOAT(ecrit_tra)
     41         zout = pdtphys * ecrit_tra
    4242c
    4343         CALL histdef(nid_tra, "phis", "Surface geop. height", "-",

  • LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/moy_undefSTD.F

    r644 r723  
    11      SUBROUTINE moy_undefSTD(nlevSTD,itap,
    2      $           ecrit_day,ecrit_mth,ecrit_hf2mth,
     2     $           dtime,ecrit_day,ecrit_mth,ecrit_hf2mth,
    33     $           tnondef,tsumSTD)
    44      IMPLICIT none
    … …  
    2828      INTEGER nlevSTD, klevSTD, itap
    2929      PARAMETER(klevSTD=17)
    30       INTEGER ecrit_day, ecrit_mth, ecrit_hf2mth
     30      REAL dtime, ecrit_day, ecrit_mth, ecrit_hf2mth
    3131c
    3232c variables locales
    … …  
    4040c calcul 1 fois par jour
    4141c
    42       IF(MOD(itap,ecrit_day).EQ.0) THEN
     42      IF(MOD(itap,NINT(ecrit_day/dtime)).EQ.0) THEN
    4343       DO k=1, nlevSTD
    4444        DO i=1, klon
    45          IF(tnondef(i,k,1).NE.FLOAT(ecrit_day)) THEN
     45         IF(tnondef(i,k,1).NE.(ecrit_day/dtime)) THEN
    4646          tsumSTD(i,k,1)=tsumSTD(i,k,1)/
    47      $    (FLOAT(ecrit_day)-tnondef(i,k,1))
     47     $    (ecrit_day/dtime-tnondef(i,k,1))
    4848         ELSE
    4949          tsumSTD(i,k,1)=1.E+20
    … …  
    5555c calcul 1 fois par mois
    5656c
    57       IF(MOD(itap,ecrit_mth).EQ.0) THEN
     57      IF(MOD(itap,NINT(ecrit_mth/dtime)).EQ.0) THEN
    5858       DO k=1, nlevSTD
    5959        DO i=1, klon
    60          IF(tnondef(i,k,2).NE.FLOAT(ecrit_mth)) THEN
     60         IF(tnondef(i,k,2).NE.ecrit_mth/dtime) THEN
    6161          tsumSTD(i,k,2)=tsumSTD(i,k,2)/
    62      $    (FLOAT(ecrit_mth)-tnondef(i,k,2))
     62     $    (ecrit_mth/dtime-tnondef(i,k,2))
    6363         ELSE
    6464          tsumSTD(i,k,2)=1.E+20
    6565         ENDIF !tnondef
    6666c
    67          IF(tnondef(i,k,3).NE.FLOAT(ecrit_hf2mth)) THEN
     67         IF(tnondef(i,k,3).NE.NINT(ecrit_hf2mth)) THEN
    6868          tsumSTD(i,k,3)=tsumSTD(i,k,3)/
    69      $    (FLOAT(ecrit_hf2mth)-tnondef(i,k,3))
     69     $    (ecrit_hf2mth-tnondef(i,k,3))
    7070         ELSE
    7171          tsumSTD(i,k,3)=1.E+20

  • LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/physiq.F

    r656 r723  
    3838c   CLEFS CPP POUR LES IO
    3939c   =====================
     40#define histins
    4041#define histhf
    4142#define histday
    4243#define histmth
    43 #define histins
    4444#define histISCCP
    4545#define histREGDYN
    … …  
    456456      INTEGER nbapp_isccp,isccppas
    457457
    458 #undef histISCCP
    459 #define histISCCP
    460 c     data ok_isccp,ecrit_isccp/.true.,0.125/     
    461 c     data ok_isccp,ecrit_isccp/.true.,1./     
    462 cIM 190504     data ok_isccp/.true./     
    463 cIM 190504 #else
    464 cIM 190504     data ok_isccp/.false./
    465 cIM 190504 #endif
    466 
    467458c sorties statistiques regime dynamique
    468459c     logical ok_regdyn
    … …  
    471462c     save ok_regdyn, ecrit_regdyn, nid_regdyn
    472463      save nid_regdyn
    473 
    474 #undef histREGDYN
    475 #define histREGDYN
    476 cIM 190504 #ifdef histREGDYN
    477 c     data ok_regdyn,ecrit_regdyn/.true.,0.125/
    478 c     data ok_regdyn,ecrit_regdyn/.true.,1./
    479 cIM 190504    data ok_regdyn/.true./
    480 cIM 190504 #else
    481 cIM 190504   data ok_regdyn/.false./
    482 cIM 190504 #endif
    483464
    484465      REAL zx_tau(kmaxm1), zx_pc(lmaxm1), zx_o500(iwmax)
    … …  
    15201501#endif
    15211502
    1522 c#undef histmthNMC
    1523 c#define histmthNMC
    15241503#ifdef histmthNMC
    15251504#include "ini_histmthNMC.h"
    … …  
    27472726c====================================================================
    27482727C
    2749       call phytrac (iflag_con,
    2750      I                   rnpb,
     2728      call phytrac (rnpb,
    27512729     I                   itap,
    27522730     I                   julien,

  • LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/phytrac.F

    r644 r723  
    44c
    55c
    6       SUBROUTINE phytrac (iflag_con,
    7      I                    rnpb,
     6      SUBROUTINE phytrac (rnpb,
    87     I                    nstep,
    98     I                    julien,
    … …  
    156155      REAL flxmass_w(klon,klev)
    157156#endif
    158       integer iflag_con
     157cIM   integer iflag_con
     158#include "clesphys.h"
    159159
    160160cAA Rem : nbtr : nombre de vrais traceurs est defini dans dimphy.h
    … …  
    266266C Variables liees a l'ecriture de la bande histoire : phytrac.nc
    267267c
    268       INTEGER ecrit_tra
    269       SAVE ecrit_tra   
     268cIM   INTEGER ecrit_tra
     269cIM   SAVE ecrit_tra   
    270270      logical ok_sync
    271271      parameter (ok_sync = .true.)
    … …  
    387387         if (debutphy) then
    388388
    389           ecrit_tra = NINT(86400./pdtphys *ecritphy)
    390           print*,'dans phytrac ',pdtphys,ecritphy,ecrit_tra
     389cIM       ecrit_tra = NINT(86400./pdtphys *ecritphy)
     390cIM       print*,'dans phytrac ',pdtphys,ecritphy,ecrit_tra
    391391
    392392         if(nbtr.lt.nqmax) then

  • LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/write_histmthNMC.h

    r644 r723  
    8484        DO i=1, klon
    8585         IF(tnondef(i,k,3).NE.1.E+20) THEN
    86           zx_tmp_fi3d(i,k) = (100.*tnondef(i,k,3))/FLOAT(ecrit_hf2mth)
     86          zx_tmp_fi3d(i,k) = (100.*tnondef(i,k,3))/ecrit_hf2mth
    8787         ELSE
    8888          zx_tmp_fi3d(i,k) = 1.E+20

  • LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/write_histrac.h

    r646 r723  
    66      ndex3d = 0
    77c
    8       zsto = pdtphys
    9       zout = pdtphys * FLOAT(ecrit_tra)
    108      itau_w = itau_phy + nstep
    11 
    129     
    13       i = NINT(zout/zsto)
    1410      CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d)
    1511      CALL histwrite(nid_tra,"phis",itau_w,zx_tmp_2d,iim*(jjm+1),ndex2d)
Note: See TracChangeset for help on using the changeset viewer.