Changeset 106


Ignore:
Timestamp:
Mar 28, 2011, 3:18:39 PM (14 years ago)
Author:
slebonnois
Message:

SL : mise a jour de phytitan pour etre conforme aux sources actuelles
utilisees sur gnome.

Location:
trunk/libf/phytitan
Files:
12 edited

Legend:

Unmodified
Added
Removed

  • trunk/libf/phytitan/ini_histday.h

    r97 r106  
    104104     .                "ave(X)", zsto,zout)
    105105c
     106      ENDIF !lev_histday.GE.2
     107c
     108c-------------------------------------------------------
     109      IF(lev_histday.GE.3) THEN
     110c
     111cccccccccccccccccc  Tracers
     112c
    106113         if (iflag_trac.eq.1) THEN
    107114          if (microfi.eq.1) then
    … …  
    121128         endif
    122129c
    123       ENDIF !lev_histday.GE.2
    124 c
    125 c-------------------------------------------------------
    126       IF(lev_histday.GE.3) THEN
    127 c
    128130cccccccccccccccccc  Radiative transfer
    129131c
    … …  
    152154     .                32, "ave(X)", zsto1,zout)
    153155c
     156         DO k=7,NSPECV,10
     157           write(str1,'(i2.2)') k
     158         CALL histdef(nid_day,"thv"//str1,"Haze Opa Vis",
     159     .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     160     .                "ave(X)",zsto1,zout)
     161         ENDDO
     162c
     163         DO k=7,NSPECV,10
     164           write(str1,'(i2.2)') k
     165         CALL histdef(nid_day,"khv"//str1,"Haze ext Vis ",
     166     .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     167     .                "ave(X)",zsto1,zout)
     168         ENDDO
     169c
     170         DO k=7,NSPECV,10
     171           write(str1,'(i2.2)') k
     172         CALL histdef(nid_day,"tgv"//str1,"Haze Opa Vis",
     173     .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     174     .                "ave(X)",zsto1,zout)
     175         ENDDO
     176c
     177         DO k=7,NSPECV,10
     178           write(str1,'(i2.2)') k
     179         CALL histdef(nid_day,"kgv"//str1,"Haze ext Vis ",
     180     .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     181     .                "ave(X)",zsto1,zout)
     182         ENDDO
     183c
     184         DO k=8,NSPECI,10
     185           write(str1,'(i2.2)') k
     186         CALL histdef(nid_day,"thi"//str1,"Haze Opa IR",
     187     .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     188     .                "ave(X)",zsto1,zout)
     189         ENDDO
     190c
     191         DO k=8,NSPECI,10
     192           write(str1,'(i2.2)') k
     193         CALL histdef(nid_day,"khi"//str1,"Haze ext IR ",
     194     .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     195     .                "ave(X)",zsto1,zout)
     196         ENDDO
     197c
     198         DO k=8,NSPECI,10
     199           write(str1,'(i2.2)') k
     200         CALL histdef(nid_day,"tgi"//str1,"Haze Opa IR",
     201     .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     202     .                "ave(X)",zsto1,zout)
     203         ENDDO
     204c
     205         DO k=8,NSPECI,10
     206           write(str1,'(i2.2)') k
     207         CALL histdef(nid_day,"kgi"//str1,"Haze ext IR ",
     208     .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     209     .                "ave(X)",zsto1,zout)
     210         ENDDO
     211c
    154212      ENDIF !lev_histday.GE.3
    155213c
    … …  
    181239     .                "ave(X)", zsto,zout)
    182240c
    183          CALL histdef(nid_day, "dtec", "Cinetic dissip dT", "K/s",
    184      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    185      .                "ave(X)", zsto,zout)
     241c        CALL histdef(nid_day, "dtec", "Cinetic dissip dT", "K/s",
     242c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     243c    .                "ave(X)", zsto,zout)
    186244c
    187245      ENDIF !lev_histday.GE.4

  • trunk/libf/phytitan/ini_histins.h

    r97 r106  
    9696     .                "inst(X)", zsto,zout)
    9797c
    98          CALL histdef(nid_ins, "duvdf", "Boundary-layer dU", "m/s2",
    99      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    100      .                "inst(X)", zsto,zout)
    101 c
    102          CALL histdef(nid_ins, "dudyn", "Dynamics dU", "m/s2",
    103      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    104      .                "inst(X)", zsto,zout)
     98c        CALL histdef(nid_ins, "duvdf", "Boundary-layer dU", "m/s2",
     99c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     100c    .                "inst(X)", zsto,zout)
     101c
     102c        CALL histdef(nid_ins, "dudyn", "Dynamics dU", "m/s2",
     103c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     104c    .                "inst(X)", zsto,zout)
    105105c
    106106      ENDIF !lev_histday.GE.2
    … …  
    154154     .                32, "inst(X)", zsto,zout)
    155155c
     156         DO k=7,NSPECV,10
     157           write(str1,'(i2.2)') k
     158         CALL histdef(nid_ins,"thv"//str1,"Haze Opa Vis",
     159     .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     160     .                "ave(X)",zsto1,zout)
     161         ENDDO
     162c
     163         DO k=7,NSPECV,10
     164           write(str1,'(i2.2)') k
     165         CALL histdef(nid_ins,"khv"//str1,"Haze ext Vis ",
     166     .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     167     .                "ave(X)",zsto1,zout)
     168         ENDDO
     169c
     170         DO k=7,NSPECV,10
     171           write(str1,'(i2.2)') k
     172         CALL histdef(nid_ins,"tgv"//str1,"Haze Opa Vis",
     173     .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     174     .                "ave(X)",zsto1,zout)
     175         ENDDO
     176c
     177         DO k=7,NSPECV,10
     178           write(str1,'(i2.2)') k
     179         CALL histdef(nid_ins,"kgv"//str1,"Haze ext Vis ",
     180     .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     181     .                "ave(X)",zsto1,zout)
     182         ENDDO
     183c
     184         DO k=8,NSPECI,10
     185           write(str1,'(i2.2)') k
     186         CALL histdef(nid_ins,"thi"//str1,"Haze Opa IR",
     187     .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     188     .                "ave(X)",zsto1,zout)
     189         ENDDO
     190c
     191         DO k=8,NSPECI,10
     192           write(str1,'(i2.2)') k
     193         CALL histdef(nid_ins,"khi"//str1,"Haze ext IR ",
     194     .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     195     .                "ave(X)",zsto1,zout)
     196         ENDDO
     197c
     198         DO k=8,NSPECI,10
     199           write(str1,'(i2.2)') k
     200         CALL histdef(nid_ins,"tgi"//str1,"Haze Opa IR",
     201     .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     202     .                "ave(X)",zsto1,zout)
     203         ENDDO
     204c
     205         DO k=8,NSPECI,10
     206           write(str1,'(i2.2)') k
     207         CALL histdef(nid_ins,"kgi"//str1,"Haze ext IR ",
     208     .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     209     .                "ave(X)",zsto1,zout)
     210         ENDDO
     211c
    156212      ENDIF !lev_histday.GE.3
    157213c
    … …  
    183239     .                "inst(X)", zsto,zout)
    184240c
    185          CALL histdef(nid_ins, "dtec", "Cinetic dissip dT", "K/s",
    186      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    187      .                "inst(X)", zsto,zout)
     241c        CALL histdef(nid_ins, "dtec", "Cinetic dissip dT", "K/s",
     242c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     243c    .                "inst(X)", zsto,zout)
    188244c
    189245c        CALL histdef(nid_ins, "dvvdf", "Boundary-layer dV", "m/s2",

  • trunk/libf/phytitan/ini_histmth.h

    r97 r106  
    5050     .                "ave(X)", zsto,zout)
    5151c
    52          CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
    53      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    54      .                "ave(X)", zsto,zout)
    55 c
    56          CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
    57      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    58      .                "ave(X)", zsto,zout)
     52c        CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
     53c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     54c    .                "ave(X)", zsto,zout)
     55c
     56c        CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
     57c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     58c    .                "ave(X)", zsto,zout)
     59c
     60c        CALL histdef(nid_mth, "cdragh", "Drag coef on T", "-",
     61c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     62c    .                "ave(X)", zsto,zout)
     63c
     64c        CALL histdef(nid_mth, "cdragm", "Drag coef on U", "-",
     65c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     66c    .                "ave(X)", zsto,zout)
    5967c
    6068      ENDIF !lev_histmth.GE.1
    … …  
    8896     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    8997     .                "ave(X)", zsto,zout)
     98c
     99c        CALL histdef(nid_mth, "Kz", "vertical diffusion coef", "m2/s",
     100c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     101c    .                "ave(X)", zsto,zout)
    90102c
    91103         CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
    … …  
    147159     .                32, "ave(X)", zsto1,zout)
    148160c
    149          DO k=3,nspecv,5
     161         DO k=7,NSPECV,10
    150162           write(str1,'(i2.2)') k
    151163         CALL histdef(nid_mth,"thv"//str1,"Haze Opa Vis",
    … …  
    154166         ENDDO
    155167c
    156          DO k=3,nspecv,5
     168         DO k=7,NSPECV,10
    157169           write(str1,'(i2.2)') k
    158170         CALL histdef(nid_mth,"khv"//str1,"Haze ext Vis ",
     171     .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     172     .                "ave(X)",zsto1,zout)
     173         ENDDO
     174c
     175         DO k=7,NSPECV,10
     176           write(str1,'(i2.2)') k
     177         CALL histdef(nid_mth,"tgv"//str1,"Haze Opa Vis",
     178     .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     179     .                "ave(X)",zsto1,zout)
     180         ENDDO
     181c
     182         DO k=7,NSPECV,10
     183           write(str1,'(i2.2)') k
     184         CALL histdef(nid_mth,"kgv"//str1,"Haze ext Vis ",
     185     .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     186     .                "ave(X)",zsto1,zout)
     187         ENDDO
     188c
     189         DO k=8,NSPECI,10
     190           write(str1,'(i2.2)') k
     191         CALL histdef(nid_mth,"thi"//str1,"Haze Opa IR",
     192     .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     193     .                "ave(X)",zsto1,zout)
     194         ENDDO
     195c
     196         DO k=8,NSPECI,10
     197           write(str1,'(i2.2)') k
     198         CALL histdef(nid_mth,"khi"//str1,"Haze ext IR ",
     199     .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     200     .                "ave(X)",zsto1,zout)
     201         ENDDO
     202c
     203         DO k=8,NSPECI,10
     204           write(str1,'(i2.2)') k
     205         CALL histdef(nid_mth,"tgi"//str1,"Haze Opa IR",
     206     .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     207     .                "ave(X)",zsto1,zout)
     208         ENDDO
     209c
     210         DO k=8,NSPECI,10
     211           write(str1,'(i2.2)') k
     212         CALL histdef(nid_mth,"kgi"//str1,"Haze ext IR ",
    159213     .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
    160214     .                "ave(X)",zsto1,zout)

  • trunk/libf/phytitan/lell.F

    r3 r106  
    105105     &  71.2,  71.4,  71.5,  71.8,  72.3,  73.1,  73.9,  74.7,  75.7,
    106106     &  76.7,  78.0,  79.2,  80.6,  82.1,  83.6,  85.5,  87.3,  88.5,
    107      &  89.5,  90.5,  91.5,  92.1,  92.7,  93.3,  93.9/
     107     &  89.5,  90.5,  91.5,  92.1,  92.9,  93.7,  94.5/
     108!    &  89.5,  90.5,  91.5,  92.1,  92.7,  93.3,  93.9/
    108109! PRESSURE IN MILLIBARS
    109110      DATA PLE/
    … …  
    132133     & 4.87E+02, 5.47E+02, 6.14E+02, 6.88E+02, 7.70E+02, 8.59E+02,
    133134     & 9.57E+02, 1.06E+03, 1.12E+03, 1.18E+03, 1.24E+03, 1.30E+03,
    134      & 1.34E+03, 1.37E+03, 1.40E+03, 1.44E+03/
     135     & 1.34E+03, 1.38E+03, 1.43E+03, 1.48E+03/
     136!    & 1.34E+03, 1.37E+03, 1.40E+03, 1.44E+03/
    135137! NUMBER DENSITY
    136138      DATA DLE/
    … …  
    159161     & 4.68E+19, 5.19E+19, 5.75E+19, 6.34E+19, 6.97E+19, 7.66E+19,
    160162     & 8.36E+19, 9.11E+19, 9.48E+19, 9.88E+19, 1.03E+20, 1.07E+20,
    161      & 1.09E+20, 1.11E+20, 1.13E+20, 1.15E+20/
     163     & 1.09E+20, 1.11E+20, 1.14E+20, 1.17E+20/
     164!    & 1.09E+20, 1.11E+20, 1.13E+20, 1.15E+20/
    162165! MEAN MOLECULAR WEIGHT
    163166      DATA XMULE/

  • trunk/libf/phytitan/lell_light.F

    r3 r106  
    111111     &  71.2,  71.4,  71.5,  71.8,  72.3,  73.1,  73.9,  74.7,  75.7,
    112112     &  76.7,  78.0,  79.2,  80.6,  82.1,  83.6,  85.5,  87.3,  88.5,
    113      &  89.5,  90.5,  91.5,  92.1,  92.7,  93.3,  93.9/
     113     &  89.5,  90.5,  91.5,  92.1,  92.9,  93.7,  94.5/
     114!    &  89.5,  90.5,  91.5,  92.1,  92.7,  93.3,  93.9/
    114115! PRESSURE IN MILLIBARS
    115116      DATA PLE/
    … …  
    138139     & 4.87E+02, 5.47E+02, 6.14E+02, 6.88E+02, 7.70E+02, 8.59E+02,
    139140     & 9.57E+02, 1.06E+03, 1.12E+03, 1.18E+03, 1.24E+03, 1.30E+03,
    140      & 1.34E+03, 1.37E+03, 1.40E+03, 1.44E+03/
     141     & 1.34E+03, 1.38E+03, 1.43E+03, 1.48E+03/
     142!    & 1.34E+03, 1.37E+03, 1.40E+03, 1.44E+03/
    141143! NUMBER DENSITY
    142144      DATA DLE/
    … …  
    165167     & 4.68E+19, 5.19E+19, 5.75E+19, 6.34E+19, 6.97E+19, 7.66E+19,
    166168     & 8.36E+19, 9.11E+19, 9.48E+19, 9.88E+19, 1.03E+20, 1.07E+20,
    167      & 1.09E+20, 1.11E+20, 1.13E+20, 1.15E+20/
     169     & 1.09E+20, 1.11E+20, 1.14E+20, 1.17E+20/
     170!    & 1.09E+20, 1.11E+20, 1.13E+20, 1.15E+20/
    168171! MEAN MOLECULAR WEIGHT
    169172      DATA XMULE/

  • trunk/libf/phytitan/pg3.F

    r97 r106  
    12041204         
    12051205
    1206          p0=.3
    1207 c        p0=1.
     1206c        p0=.3
     1207         p0=1.
    12081208         do i=1,nz
    12091209          if (pb(i).lt.p0.and.pb(i+1).gt.p0) zalt0=(z(i)+z(i+1))/2.

  • trunk/libf/phytitan/physiq.F

    r105 r106  
    418418      ok_sync=.TRUE.
    419419
    420       bilansmc = 1
     420      bilansmc = 0
    421421      ballons  = 0
    422422
    … …  
    748748      DO l=1,klev
    749749         DO i=1,klon
    750             zzlay(i,l)=zphi(i,l)/RG
     750c           zzlay(i,l)=zphi(i,l)/RG
     751c SI ON TIENT COMPTE DE LA VARIATION DE G AVEC L'ALTITUDE:
     752            zzlay(i,l)=RG*RA*RA/(RG*RA-zphi(i,l))-RA
    751753         ENDDO
    752754      ENDDO
    753755      DO i=1,klon
    754          zzlev(i,1)=0.
     756c        zzlev(i,1)=0.
     757c CORRECTION 13/01/2011 
     758c (correspond a la position de la surface en ce point vs RA)
     759         zzlev(i,1)=pphis(i)/RG
    755760      ENDDO
    756761      DO l=2,klev
    … …  
    764769         zzlev(i,klev+1)=zzlay(i,klev)+(zzlay(i,klev)-zzlev(i,klev))
    765770      ENDDO
     771
     772c- - - - - - - - - - - - - - - -
     773c DIAGNOSTIQUE GRILLE VERTICALE
     774c- - - - - - - - - - - - - - - -
     775c     print*,"DIAGNOSTIQUE GRILLE VERTICALE"
     776c     i=klon/2
     777c     print*,"Niveau  Pression  Altitude    (lev puis lay)"
     778c     do l=1,klev
     779c      print*,l,paprs(i,l),zzlev(i,l)
     780c      print*,l,pplay(i,l),zzlay(i,l)
     781c     enddo
     782c     print*,klev+1,paprs(i,klev+1),zzlev(i,klev+1)
     783c     stop
    766784
    767785c====================================================================

  • trunk/libf/phytitan/radtitan.F

    r104 r106  
    253253210   CONTINUE
    254254
    255       IF(press(nlevel-1).GE.1.44) then
     255c  a cause du tableau predefini dans lell.F (et lell_light.F)
     256c     IF(press(nlevel-1).GE.1.44) then
     257      IF(press(nlevel-1).GE.1.48) then
    256258           STOP'pression au sol trop grande'
    257259          PRINT*,'pression au sol trop grande'
    258260      endif
    259261
    260 c      PRESS(nlevel)=1.44
    261 c      XCORR=1.44/PRESS(nlevel)
     262c      PRESS(nlevel)=1.48
     263c      XCORR=1.48/PRESS(nlevel)
    262264c     DO 211 J=1,NLEVEL
    263265c        PRESS(J)=XCORR*PRESS(J)

  • trunk/libf/phytitan/suphec.F

    r97 r106  
    5555c TITAN
    5656      RSIYEA=9.28e8      ! 673 jTitan
    57       RSIDAY=1.37889e6   ! 15.96 j
     57c LEQUEL DES 2 ?? A VERIFIER !!
     58c     RSIDAY=1.37889e6   ! 15.96 j
     59      RDAY=1.37889e6   ! 15.96 j
    5860      ROMEGA=2.*RPI/RSIDAY
    5961c 1/(duree du jour) = 1/(periode rotation) - 1/(periode revolution)
    60       RDAY=RSIDAY/(1.-RSIDAY/RSIYEA)
     62c     RDAY=RSIDAY/(1.-RSIDAY/RSIYEA)
     63      RSIDAY=RDAY/(1.+RDAY/RSIYEA)
    6164      REA=1.5e12
    6265      REPSM=0.  ! 0. veut dire qu'on commence au point vernal

  • trunk/libf/phytitan/write_histday.h

    r97 r106  
    9494     .                                   iim*jjmp1*klev,ndex3d)
    9595c
     96      ENDIF !lev_histday.GE.2
     97c
     98c-------------------------------------------------------
     99      IF(lev_histday.GE.3) THEN
     100c
     101cccccccccccccccccc  Tracers
     102c
    96103         if (iflag_trac.eq.1) THEN
    97104          if (microfi.eq.1) then
    98105           DO iq=1,nmicro
    99106       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qaer(1,1,iq), zx_tmp_3d)
    100        CALL histwrite(nid_day,tname(iq),itau_w,zx_tmp_3d,
     107       z4_tmp_3d=zx_tmp_3d
     108       CALL histwrite(nid_day,tname(iq),itau_w,z4_tmp_3d,
    101109     .                                   iim*jjmp1*klev,ndex3d)
    102110           ENDDO
    … …  
    105113           DO iq=nmicro+1,nqmax
    106114       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,tr_seri(1,1,iq),zx_tmp_3d)
    107        CALL histwrite(nid_day,tname(iq),itau_w,zx_tmp_3d,
     115       z4_tmp_3d=zx_tmp_3d
     116       CALL histwrite(nid_day,tname(iq),itau_w,z4_tmp_3d,
    108117     .                                   iim*jjmp1*klev,ndex3d)
    109118           ENDDO
    … …  
    111120         endif
    112121c
    113       ENDIF !lev_histday.GE.2
    114 c
    115 c-------------------------------------------------------
    116       IF(lev_histday.GE.3) THEN
    117 c
    118122cccccccccccccccccc  Radiative transfer
    119123c
    … …  
    141145     .                                   iim*jjmp1*klev,ndex3d)
    142146c
     147c 3D adding Tau and k  (31/08/10)
     148c
     149       do k=7,NSPECV,10
     150         do i=1,klon
     151         do l=1,klev
     152           t_tauhvd(i,l)=TAUHVD(i,klev-l+1,k)
     153         enddo
     154         enddo
     155         write(str1,'(i2.2)') k
     156      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
     157      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     158      z4_tmp_3d=zx_tmp_3d
     159      CALL histwrite(nid_day,"thv"//str1,itau_w,z4_tmp_3d,
     160     .                                   iim*jjmp1*klev,ndex3d)
     161       enddo      ! fin boucle NSPECV
     162
     163       do k=7,NSPECV,10
     164         do i=1,klon
     165         do l=1,klev
     166         if(l.ne.klev)
     167     s    t_khvd(i,l)=TAUHVD(i,klev-l+1,k)
     168     s    -TAUHVD(i,klev-l+1-1,k)
     169
     170          if(l.eq.klev)
     171     s    t_khvd(i,l)=TAUHVD(i,klev-l+1,k)
     172
     173         t_khvd(i,l)=t_khvd(i,l)/(zzlev(i,l+1)-zzlev(i,l))
     174         enddo
     175         enddo
     176         write(str1,'(i2.2)') k
     177      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
     178      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     179      z4_tmp_3d=zx_tmp_3d
     180      CALL histwrite(nid_day,"khv"//str1,itau_w,z4_tmp_3d,
     181     .                                   iim*jjmp1*klev,ndex3d)
     182       enddo      ! fin boucle NSPECV
     183
     184       do k=7,NSPECV,10
     185         do i=1,klon
     186         do l=1,klev
     187           t_tauhvd(i,l)=TAUGVD(i,klev-l+1,k)
     188         enddo
     189         enddo
     190         write(str1,'(i2.2)') k
     191      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
     192      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     193      z4_tmp_3d=zx_tmp_3d
     194      CALL histwrite(nid_day,"tgv"//str1,itau_w,z4_tmp_3d,
     195     .                                   iim*jjmp1*klev,ndex3d)
     196       enddo      ! fin boucle NSPECV
     197
     198       do k=7,NSPECV,10
     199         do i=1,klon
     200         do l=1,klev
     201         if(l.ne.klev)
     202     s    t_khvd(i,l)=TAUGVD(i,klev-l+1,k)
     203     s    -TAUGVD(i,klev-l+1-1,k)
     204
     205          if(l.eq.klev)
     206     s    t_khvd(i,l)=TAUGVD(i,klev-l+1,k)
     207
     208         t_khvd(i,l)=t_khvd(i,l)/(zzlev(i,l+1)-zzlev(i,l))
     209         enddo
     210         enddo
     211         write(str1,'(i2.2)') k
     212      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
     213      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     214      z4_tmp_3d=zx_tmp_3d
     215      CALL histwrite(nid_day,"kgv"//str1,itau_w,z4_tmp_3d,
     216     .                                   iim*jjmp1*klev,ndex3d)
     217       enddo      ! fin boucle NSPECV
     218
     219       do k=8,NSPECI,10
     220         do i=1,klon
     221         do l=1,klev
     222           t_tauhvd(i,l)=TAUHID(i,klev-l+1,k)
     223         enddo
     224         enddo
     225         write(str1,'(i2.2)') k
     226      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
     227      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     228      z4_tmp_3d=zx_tmp_3d
     229      CALL histwrite(nid_day,"thi"//str1,itau_w,z4_tmp_3d,
     230     .                                   iim*jjmp1*klev,ndex3d)
     231       enddo      ! fin boucle NSPECI
     232
     233       do k=8,NSPECI,10
     234         do i=1,klon
     235         do l=1,klev
     236         if(l.ne.klev)
     237     s    t_khvd(i,l)=TAUHID(i,klev-l+1,k)
     238     s    -TAUHID(i,klev-l+1-1,k)
     239
     240          if(l.eq.klev)
     241     s    t_khvd(i,l)=TAUHID(i,klev-l+1,k)
     242
     243         t_khvd(i,l)=t_khvd(i,l)/(zzlev(i,l+1)-zzlev(i,l))
     244         enddo
     245         enddo
     246         write(str1,'(i2.2)') k
     247      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
     248      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     249      z4_tmp_3d=zx_tmp_3d
     250      CALL histwrite(nid_day,"khi"//str1,itau_w,z4_tmp_3d,
     251     .                                   iim*jjmp1*klev,ndex3d)
     252       enddo      ! fin boucle NSPECI
     253
     254       do k=8,NSPECI,10
     255         do i=1,klon
     256         do l=1,klev
     257           t_tauhvd(i,l)=TAUGID(i,klev-l+1,k)
     258         enddo
     259         enddo
     260         write(str1,'(i2.2)') k
     261      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
     262      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     263      z4_tmp_3d=zx_tmp_3d
     264      CALL histwrite(nid_day,"tgi"//str1,itau_w,z4_tmp_3d,
     265     .                                   iim*jjmp1*klev,ndex3d)
     266       enddo      ! fin boucle NSPECI
     267
     268       do k=8,NSPECI,10
     269         do i=1,klon
     270         do l=1,klev
     271         if(l.ne.klev)
     272     s    t_khvd(i,l)=TAUGID(i,klev-l+1,k)
     273     s    -TAUGID(i,klev-l+1-1,k)
     274
     275          if(l.eq.klev)
     276     s    t_khvd(i,l)=TAUGID(i,klev-l+1,k)
     277
     278         t_khvd(i,l)=t_khvd(i,l)/(zzlev(i,l+1)-zzlev(i,l))
     279         enddo
     280         enddo
     281         write(str1,'(i2.2)') k
     282      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
     283      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     284      z4_tmp_3d=zx_tmp_3d
     285      CALL histwrite(nid_day,"kgi"//str1,itau_w,z4_tmp_3d,
     286     .                                   iim*jjmp1*klev,ndex3d)
     287       enddo      ! fin boucle NSPECI
     288
    143289      ENDIF !lev_histday.GE.3
    144290c
    … …  
    177323     .                                   iim*jjmp1*klev,ndex3d)
    178324c K/s     
    179       zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
    180       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    181       CALL histwrite(nid_day,"dtec",itau_w,zx_tmp_3d,
    182      .                                   iim*jjmp1*klev,ndex3d)
     325c     zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
     326c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     327c     CALL histwrite(nid_day,"dtec",itau_w,zx_tmp_3d,
     328c    .                                   iim*jjmp1*klev,ndex3d)
    183329c
    184330      ENDIF !lev_histday.GE.4

  • trunk/libf/phytitan/write_histins.h

    r97 r106  
    8686      CALL histwrite(nid_ins,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    8787c
    88       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d)
    89       CALL histwrite(nid_ins,"duvdf",itau_w,zx_tmp_3d,
    90      .                                   iim*jjmp1*klev,ndex3d)
    91 c
    92       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_dyn, zx_tmp_3d)
    93       CALL histwrite(nid_ins,"dudyn",itau_w,zx_tmp_3d,
    94      .                                   iim*jjmp1*klev,ndex3d)
     88c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d)
     89c     CALL histwrite(nid_ins,"duvdf",itau_w,zx_tmp_3d,
     90c    .                                   iim*jjmp1*klev,ndex3d)
     91c
     92c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_dyn, zx_tmp_3d)
     93c     CALL histwrite(nid_ins,"dudyn",itau_w,zx_tmp_3d,
     94c    .                                   iim*jjmp1*klev,ndex3d)
    9595c
    9696      ENDIF !lev_histday.GE.2
    … …  
    143143     .                                   iim*jjmp1*klev,ndex3d)
    144144c
     145c 3D adding Tau and k  (31/08/10)
     146c
     147       do k=7,NSPECV,10
     148         do i=1,klon
     149         do l=1,klev
     150           t_tauhvd(i,l)=TAUHVD(i,klev-l+1,k)
     151         enddo
     152         enddo
     153         write(str1,'(i2.2)') k
     154      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
     155      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     156      z4_tmp_3d=zx_tmp_3d
     157      CALL histwrite(nid_ins,"thv"//str1,itau_w,z4_tmp_3d,
     158     .                                   iim*jjmp1*klev,ndex3d)
     159       enddo      ! fin boucle NSPECV
     160
     161       do k=7,NSPECV,10
     162         do i=1,klon
     163         do l=1,klev
     164         if(l.ne.klev)
     165     s    t_khvd(i,l)=TAUHVD(i,klev-l+1,k)
     166     s    -TAUHVD(i,klev-l+1-1,k)
     167
     168          if(l.eq.klev)
     169     s    t_khvd(i,l)=TAUHVD(i,klev-l+1,k)
     170
     171         t_khvd(i,l)=t_khvd(i,l)/(zzlev(i,l+1)-zzlev(i,l))
     172         enddo
     173         enddo
     174         write(str1,'(i2.2)') k
     175      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
     176      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     177      z4_tmp_3d=zx_tmp_3d
     178      CALL histwrite(nid_ins,"khv"//str1,itau_w,z4_tmp_3d,
     179     .                                   iim*jjmp1*klev,ndex3d)
     180       enddo      ! fin boucle NSPECV
     181
     182       do k=7,NSPECV,10
     183         do i=1,klon
     184         do l=1,klev
     185           t_tauhvd(i,l)=TAUGVD(i,klev-l+1,k)
     186         enddo
     187         enddo
     188         write(str1,'(i2.2)') k
     189      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
     190      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     191      z4_tmp_3d=zx_tmp_3d
     192      CALL histwrite(nid_ins,"tgv"//str1,itau_w,z4_tmp_3d,
     193     .                                   iim*jjmp1*klev,ndex3d)
     194       enddo      ! fin boucle NSPECV
     195
     196       do k=7,NSPECV,10
     197         do i=1,klon
     198         do l=1,klev
     199         if(l.ne.klev)
     200     s    t_khvd(i,l)=TAUGVD(i,klev-l+1,k)
     201     s    -TAUGVD(i,klev-l+1-1,k)
     202
     203          if(l.eq.klev)
     204     s    t_khvd(i,l)=TAUGVD(i,klev-l+1,k)
     205
     206         t_khvd(i,l)=t_khvd(i,l)/(zzlev(i,l+1)-zzlev(i,l))
     207         enddo
     208         enddo
     209         write(str1,'(i2.2)') k
     210      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
     211      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     212      z4_tmp_3d=zx_tmp_3d
     213      CALL histwrite(nid_ins,"kgv"//str1,itau_w,z4_tmp_3d,
     214     .                                   iim*jjmp1*klev,ndex3d)
     215       enddo      ! fin boucle NSPECV
     216
     217       do k=8,NSPECI,10
     218         do i=1,klon
     219         do l=1,klev
     220           t_tauhvd(i,l)=TAUHID(i,klev-l+1,k)
     221         enddo
     222         enddo
     223         write(str1,'(i2.2)') k
     224      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
     225      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     226      z4_tmp_3d=zx_tmp_3d
     227      CALL histwrite(nid_ins,"thi"//str1,itau_w,z4_tmp_3d,
     228     .                                   iim*jjmp1*klev,ndex3d)
     229       enddo      ! fin boucle NSPECI
     230
     231       do k=8,NSPECI,10
     232         do i=1,klon
     233         do l=1,klev
     234         if(l.ne.klev)
     235     s    t_khvd(i,l)=TAUHID(i,klev-l+1,k)
     236     s    -TAUHID(i,klev-l+1-1,k)
     237
     238          if(l.eq.klev)
     239     s    t_khvd(i,l)=TAUHID(i,klev-l+1,k)
     240
     241         t_khvd(i,l)=t_khvd(i,l)/(zzlev(i,l+1)-zzlev(i,l))
     242         enddo
     243         enddo
     244         write(str1,'(i2.2)') k
     245      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
     246      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     247      z4_tmp_3d=zx_tmp_3d
     248      CALL histwrite(nid_ins,"khi"//str1,itau_w,z4_tmp_3d,
     249     .                                   iim*jjmp1*klev,ndex3d)
     250       enddo      ! fin boucle NSPECI
     251
     252       do k=8,NSPECI,10
     253         do i=1,klon
     254         do l=1,klev
     255           t_tauhvd(i,l)=TAUGID(i,klev-l+1,k)
     256         enddo
     257         enddo
     258         write(str1,'(i2.2)') k
     259      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
     260      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     261      z4_tmp_3d=zx_tmp_3d
     262      CALL histwrite(nid_ins,"tgi"//str1,itau_w,z4_tmp_3d,
     263     .                                   iim*jjmp1*klev,ndex3d)
     264       enddo      ! fin boucle NSPECI
     265
     266       do k=8,NSPECI,10
     267         do i=1,klon
     268         do l=1,klev
     269         if(l.ne.klev)
     270     s    t_khvd(i,l)=TAUGID(i,klev-l+1,k)
     271     s    -TAUGID(i,klev-l+1-1,k)
     272
     273          if(l.eq.klev)
     274     s    t_khvd(i,l)=TAUGID(i,klev-l+1,k)
     275
     276         t_khvd(i,l)=t_khvd(i,l)/(zzlev(i,l+1)-zzlev(i,l))
     277         enddo
     278         enddo
     279         write(str1,'(i2.2)') k
     280      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
     281      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     282      z4_tmp_3d=zx_tmp_3d
     283      CALL histwrite(nid_ins,"kgi"//str1,itau_w,z4_tmp_3d,
     284     .                                   iim*jjmp1*klev,ndex3d)
     285       enddo      ! fin boucle NSPECI
     286
    145287      ENDIF !lev_histday.GE.3
    146288c
    … …  
    179321     .                                   iim*jjmp1*klev,ndex3d)
    180322c K/s     
    181       zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
    182       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    183       CALL histwrite(nid_ins,"dtec",itau_w,zx_tmp_3d,
    184      .                                   iim*jjmp1*klev,ndex3d)
     323c     zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
     324c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     325c     CALL histwrite(nid_ins,"dtec",itau_w,zx_tmp_3d,
     326c    .                                   iim*jjmp1*klev,ndex3d)
    185327c
    186328c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_vdf, zx_tmp_3d)

  • trunk/libf/phytitan/write_histmth.h

    r105 r106  
    4848c     CALL histwrite(nid_mth,"ve",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    4949c
     50c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
     51c     z4_tmp_2d=zx_tmp_2d
     52c     CALL histwrite(nid_mth,"cdragh",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     53c
     54c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
     55c     z4_tmp_2d=zx_tmp_2d
     56c     CALL histwrite(nid_mth,"cdragm",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     57c
    5058      ENDIF !lev_histmth.GE.1
    5159c
    … …  
    7886      CALL histwrite(nid_mth,"vitw",itau_w,zx_tmp_3d,
    7987     .                                   iim*jjmp1*klev,ndex3d)
     88c
     89c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ycoefh, zx_tmp_3d)
     90c     z4_tmp_3d=zx_tmp_3d
     91c     CALL histwrite(nid_mth,"Kz",itau_w,z4_tmp_3d,
     92c    .                                   iim*jjmp1*klev,ndex3d)
    8093c
    8194      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
    … …  
    136149c
    137150c 3D adding Tau and k  (31/08/10)
    138 c
    139        do k=3,nspecv,5
    140          do i=1,klon
    141          do l=1,klev
    142            t_tauhvd(i,j)=TAUHVD(i,klev-l+1,k)
     151c 
     152       do k=7,NSPECV,10
     153         do i=1,klon
     154         do l=1,klev
     155           t_tauhvd(i,l)=TAUHVD(i,klev-l+1,k)
    143156         enddo
    144157         enddo
    … …  
    146159      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
    147160      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    148       CALL histwrite(nid_mth,"thv"//str1,itau_w,zx_tmp_3d,
    149      .                                   iim*jjmp1*klev,ndex3d)
    150        enddo      ! fin boucle nspecv
    151 
    152        do k=3,nspecv,5
     161      z4_tmp_3d=zx_tmp_3d
     162      CALL histwrite(nid_day,"thv"//str1,itau_w,z4_tmp_3d,
     163     .                                   iim*jjmp1*klev,ndex3d)
     164       enddo      ! fin boucle NSPECV
     165
     166       do k=7,NSPECV,10
    153167         do i=1,klon
    154168         do l=1,klev
    … …  
    166180      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
    167181      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    168       CALL histwrite(nid_mth,"khv"//str1,itau_w,zx_tmp_3d,
    169      .                                   iim*jjmp1*klev,ndex3d)
    170        enddo      ! fin boucle nspecv
     182      z4_tmp_3d=zx_tmp_3d
     183      CALL histwrite(nid_day,"khv"//str1,itau_w,z4_tmp_3d,
     184     .                                   iim*jjmp1*klev,ndex3d)
     185       enddo      ! fin boucle NSPECV
     186
     187       do k=7,NSPECV,10
     188         do i=1,klon
     189         do l=1,klev
     190           t_tauhvd(i,l)=TAUGVD(i,klev-l+1,k)
     191         enddo
     192         enddo
     193         write(str1,'(i2.2)') k
     194      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
     195      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     196      z4_tmp_3d=zx_tmp_3d
     197      CALL histwrite(nid_day,"tgv"//str1,itau_w,z4_tmp_3d,
     198     .                                   iim*jjmp1*klev,ndex3d)
     199       enddo      ! fin boucle NSPECV
     200
     201       do k=7,NSPECV,10
     202         do i=1,klon
     203         do l=1,klev
     204         if(l.ne.klev)
     205     s    t_khvd(i,l)=TAUGVD(i,klev-l+1,k)
     206     s    -TAUGVD(i,klev-l+1-1,k)
     207
     208          if(l.eq.klev)
     209     s    t_khvd(i,l)=TAUGVD(i,klev-l+1,k)
     210
     211         t_khvd(i,l)=t_khvd(i,l)/(zzlev(i,l+1)-zzlev(i,l))
     212         enddo
     213         enddo
     214         write(str1,'(i2.2)') k
     215      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
     216      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     217      z4_tmp_3d=zx_tmp_3d
     218      CALL histwrite(nid_day,"kgv"//str1,itau_w,z4_tmp_3d,
     219     .                                   iim*jjmp1*klev,ndex3d)
     220       enddo      ! fin boucle NSPECV
     221
     222       do k=8,NSPECI,10
     223         do i=1,klon
     224         do l=1,klev
     225           t_tauhvd(i,l)=TAUHID(i,klev-l+1,k)
     226         enddo
     227         enddo
     228         write(str1,'(i2.2)') k
     229      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
     230      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     231      z4_tmp_3d=zx_tmp_3d
     232      CALL histwrite(nid_day,"thi"//str1,itau_w,z4_tmp_3d,
     233     .                                   iim*jjmp1*klev,ndex3d)
     234       enddo      ! fin boucle NSPECI
     235
     236       do k=8,NSPECI,10
     237         do i=1,klon
     238         do l=1,klev
     239         if(l.ne.klev)
     240     s    t_khvd(i,l)=TAUHID(i,klev-l+1,k)
     241     s    -TAUHID(i,klev-l+1-1,k)
     242
     243          if(l.eq.klev)
     244     s    t_khvd(i,l)=TAUHID(i,klev-l+1,k)
     245
     246         t_khvd(i,l)=t_khvd(i,l)/(zzlev(i,l+1)-zzlev(i,l))
     247         enddo
     248         enddo
     249         write(str1,'(i2.2)') k
     250      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
     251      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     252      z4_tmp_3d=zx_tmp_3d
     253      CALL histwrite(nid_day,"khi"//str1,itau_w,z4_tmp_3d,
     254     .                                   iim*jjmp1*klev,ndex3d)
     255       enddo      ! fin boucle NSPECI
     256
     257       do k=8,NSPECI,10
     258         do i=1,klon
     259         do l=1,klev
     260           t_tauhvd(i,l)=TAUGID(i,klev-l+1,k)
     261         enddo
     262         enddo
     263         write(str1,'(i2.2)') k
     264      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
     265      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     266      z4_tmp_3d=zx_tmp_3d
     267      CALL histwrite(nid_day,"tgi"//str1,itau_w,z4_tmp_3d,
     268     .                                   iim*jjmp1*klev,ndex3d)
     269       enddo      ! fin boucle NSPECI
     270
     271       do k=8,NSPECI,10
     272         do i=1,klon
     273         do l=1,klev
     274         if(l.ne.klev)
     275     s    t_khvd(i,l)=TAUGID(i,klev-l+1,k)
     276     s    -TAUGID(i,klev-l+1-1,k)
     277
     278          if(l.eq.klev)
     279     s    t_khvd(i,l)=TAUGID(i,klev-l+1,k)
     280
     281         t_khvd(i,l)=t_khvd(i,l)/(zzlev(i,l+1)-zzlev(i,l))
     282         enddo
     283         enddo
     284         write(str1,'(i2.2)') k
     285      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
     286      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
     287      z4_tmp_3d=zx_tmp_3d
     288      CALL histwrite(nid_day,"kgi"//str1,itau_w,z4_tmp_3d,
     289     .                                   iim*jjmp1*klev,ndex3d)
     290       enddo      ! fin boucle NSPECI
    171291
    172292      ENDIF !lev_histmth.GE.3
Note: See TracChangeset for help on using the changeset viewer.