Ignore:
Timestamp:
Aug 28, 2007, 5:10:14 PM (17 years ago)
Author:
lmdzadmin
Message:

Correction bug sorties ISCCP (7pres X 7tau) + ENSEMBLES + ajout freq_ISCCP/ecrit_ISCCP
IM

Location:
LMDZ4/trunk/libf/phytherm
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/trunk/libf/phytherm/calcul_simulISCCP.h

    r814 r832  
    1010c
    1111      nbapp_isccp=30 !appel toutes les 15h
    12       isccppas=NINT((itap*dtime)/3600.) !Nb. d'heures de la physique
     12cIM 170107      isccppas=NINT((itap*dtime)/3600.) !Nb. d'heures de la physique
    1313      freqin_pdt(n)=ifreq_isccp(n)
    1414c
  • LMDZ4/trunk/libf/phytherm/ini_histISCCP.h

    r814 r832  
    1515cIM 300505     zstophy = dtime
    1616c appel du simulateur toutes les 3heures
    17           zcals(1) = dtime *6.  !toutes les 3h (en s)
     17!IM on lit la frequence d'appel dans physiq.def
     18!         zcals(1) = dtime *6.  !toutes les 3h (en s)
     19          zcals(1) = freq_ISCCP !toutes les freq_ISCCP secondes
    1820        DO n=1, napisccp
    1921          zcalh(n) = zcals(n)/3600. !stoutes les Xh (en heures)
     
    2931c       IF(freqout_isccp.EQ.1.) THEN
    3032c ecriture jounaliere
    31           zout_isccp(1) = ecrit_day !(en s)
     33!IM on ecrit les resultats du simulateur ISCCP toutes les
     34! ecrit_ISCCP secondes      zout_isccp(1) = ecrit_day !(en s)
     35          zout_isccp(1) = ecrit_ISCCP !(en s)
    3236c ecriture mensuelle
    3337c         zout = dtime * ecrit_mth !(en s)
     
    3741c le nombre de sous-colonnes ncol : ncol=(100.*zcalh)/zoutd
    3842          ncol(n)=NINT((100.*zcalh(n))/zoutj(n))
    39           PRINT*,'n ncol',n,ncol(n)
     43          IF(ncol(n).GT.ncolmx) THEN
     44           PRINT*,'Warning: Augmenter le nombre colonnes du simulateur'
     45           PRINT*,'         ISCCP ncol=', ncol,' ncolmx=',ncolmx
     46c          PRINT*,'n ncol',n,ncol(n)
     47           CALL abort
     48          ENDIF
    4049c
    4150        DO l=1, ncol(n)
     
    8089cym          enddo
    8190c
    82           print*,'On stoke le fichier histISCCP instantanne sur ',
     91          print*,'On stoke le fichier histISCCP sur ',
    8392     s   imin_ins,imax_ins,jmin_ins,jmax_ins
    8493cym          print*,'On stoke le fichier histISCCP instantanne sur ',
     
    94103         CALL histbeg_phy("histISCCP.nc", itau_phy, zjulian, dtime,
    95104     .                 nhori, nid_isccp)
    96 
    97105        ENDIF !(1.EQ.0) THEN
    98106c
     
    106114         CALL histbeg_phy("histISCCP.nc", itau_phy, zjulian, dtime,
    107115     .                 nhori, nid_isccp)
    108 
    109116c
    110117        IF(type_run.EQ."ENSP".OR.type_run.EQ."CLIM") THEN
     
    128135c
    129136         DO k=1, kmaxm1
    130           CALL histdef(nid_isccp, "cldISCCP_"//taulev(k)//typinout(n),
     137          CALL histdef(nid_isccp, "cldISCCP_"//taulev(k)//verticaxe(n),
    131138     .                "LMDZ ISCCP cld", "%",
    132139     .                iim, jj_nb,nhori,lmaxm1,1,lmaxm1,nvert,32,
     
    134141         ENDDO
    135142c
    136          CALL histdef(nid_isccp, "nsunlit"//typinout(n),
     143         CALL histdef(nid_isccp, "nsunlit"//verticaxe(n),
    137144     .                "Nb of calls with sunlit ", "%",
    138145     .                iim, jj_nb,nhori,1,1,1,-99,32,
     
    150157           DO l=1, lmaxm1
    151158c
    152            CALL histdef(nid_isccp, pclev(l)//taulev(k)//typinout(n),
     159           CALL histdef(nid_isccp, pclev(l)//taulev(k)//verticaxe(n),
    153160     .                "LMDZ ISCCP cld "//cnameisccp(l,k), "%",
    154161     .                iim, jj_nb,nhori,1,1,1,-99,32,
     
    159166c
    160167c         print*,'n=',n,' avant histdef(..Nb of calls sunlit'
    161           CALL histdef(nid_isccp, "nsunlit"//typinout(n),
     168          CALL histdef(nid_isccp, "nsunlit"//verticaxe(n),
    162169     .                "Nb of calls with sunlit ", "%",
    163170     .                iim, jj_nb,nhori,1,1,1,-99,32,
  • LMDZ4/trunk/libf/phytherm/write_histISCCP.h

    r814 r832  
    2020cIM: champ 3d : (lon,lat,pres) pour un tau fixe
    2121c
    22        CALL histwrite_phy(nid_isccp,"cldISCCP_"//taulev(k)//typinout(n),
     22       CALL histwrite_phy(nid_isccp,"cldISCCP_"//taulev(k)//verticaxe(n),
    2323     .                  itau_w,zx_tmp_fi3d)
    2424        ENDDO !k
    2525c
    2626cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1,nbsunlit(1,:,n),zx_tmp_2d)
    27         CALL histwrite_phy(nid_isccp,"nsunlit"//typinout(n),itau_w,
     27        CALL histwrite_phy(nid_isccp,"nsunlit"//verticaxe(n),itau_w,
    2828     .                 nbsunlit(1,:,n))
    2929c
     
    5555cIM: champ 2d : (lon,lat) pour un tau et une pc fixes
    5656c
    57          CALL histwrite_phy(nid_isccp,pclev(l)//taulev(k)//typinout(n),
    58      .                  itau_w,fq_is_true(:,:,l,n))
     57         CALL histwrite_phy(nid_isccp,pclev(l)//taulev(k)//verticaxe(n),
     58     .                  itau_w,fq_is_true(:,k,l,n))
    5959         ENDDO !l
    6060        ENDDO !k
     
    6262c       print*,'n=',n,' write_ISCCP avant nbsunlit'
    6363cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1,nbsunlit(1,:,n),zx_tmp_2d)
    64         CALL histwrite_phy(nid_isccp,"nsunlit"//typinout(n),
     64        CALL histwrite_phy(nid_isccp,"nsunlit"//verticaxe(n),
    6565     .                 itau_w,nbsunlit(1,:,n))
    6666c
Note: See TracChangeset for help on using the changeset viewer.