Ignore:
Timestamp:
Dec 15, 2011, 3:07:31 AM (13 years ago)
Author:
slebonnois
Message:

Update of Titan physics for clouds.

Location:
trunk/LMDZ.TITAN/libf/phytitan
Files:
11 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.TITAN/libf/phytitan/brume3D.F

    r175 r474  
    113113           do j=1,nrad
    114114             v1=vitesse(i,j,0)
     115c             ho que c'est moche ! -> taused = RT/(Mn2*g)*1/vaer = H/v
    115116             taused(i,j)=(8.314*t(i)/28.e-3/1.35)/v1
    116117           enddo
    … …  
    126127         call sedif(dice3)
    127128c        En theorie, dice3 est NEGATIF (en sedimentant on ne fait que perdre des aerosols)
    128 c        Le precipitations sont comptées positivement.
    129          precip(ihor,5)=-dice3/rhol    ! m3/m2=m
     129c        Les precipitations sont comptees positivement. (ET ON NE PREND QUE DES VALEURS POSITIVES)
     130         precip(ihor,5)=AMAX1(-dice3/rhol,0.)    ! m3/m2=m
    130131
    131132         li=3-li

  • trunk/LMDZ.TITAN/libf/phytitan/conf_phys.F90

    r175 r474  
    391391  clouds = 1
    392392  call getin('clouds',clouds)
    393   if (microfi.ne.1) clouds = 0      ! On  ne fait pas de nuages sans microphysique !
     393  if (microfi.lt.1) clouds = 0      ! On  ne fait pas de nuages sans microphysique !
    394394
    395395!

  • trunk/LMDZ.TITAN/libf/phytitan/diagmuphy.h

    r175 r474  
    1313       REAL tau_aer(ngd,llm,nrad)
    1414! ---- Production de glace (negatif si disparition)
    15        REAL solesp(ngd,llm+1,3)
     15       REAL solesp(ngd,llm,3)
     16! ---- Evaporation CH4
     17       REAL evapch4(ngd)
    1618! ---- occurences des nuages
    1719       REAL occcld_m(ngd,llm,12)
    … …  
    2527       REAL satc2h2(ngd,llm)
    2628! ---- precipitations (CH4, C2H6, C2H2, noyaux, aerosols)
    27        REAL prec(ngd,5)
     29       REAL precip(ngd,5)
    2830! ---- rayon moyen des gouttes
    2931       REAL rmcloud(ngd,llm)
    3032   
    3133       common/diagmuphy/flxesp_i,tau_drop,tau_aer,                      &
    32      &                  solesp,occcld_m,occcld,                         &
     34     &                  solesp,evapch4,occcld_m,occcld,                 &
    3335     &                  satch4,satc2h6,satc2h2,                         &
    34      &                  prec,rmcloud
     36     &                  precip,rmcloud
    3537

  • trunk/LMDZ.TITAN/libf/phytitan/ini_histday.h

    r306 r474  
    114114     .                    "ave(X)", zsto,zout)
    115115
    116              if (clouds.eq.1) then
    117 
     116            if (clouds.eq.1) then
    118117             CALL histdef(nid_day, "qnoy","nb tot noy" , "n/m2",
    119118     .                    iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    … …  
    142141c ----- RESERVOIR DE SURFACE
    143142               CALL histdef(nid_day, "reserv", "Reservoir surface","m",
     143     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     144     .                "ave(X)", zsto,zout)
     145c --------------
     146c ----- ECHANGE GAZ SURF/ATM (evaporation)
     147               CALL histdef(nid_day, "evapch4", "Evaporation CH4","m",
    144148     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    145149     .                "ave(X)", zsto,zout)
    … …  
    171175     .                "m",iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    172176     .                "ave(X)", zsto,zout)
    173              endif
     177            endif
    174178          endif
    175179c --------------

  • trunk/LMDZ.TITAN/libf/phytitan/ini_histmth.h

    r306 r474  
    118118     .                    "ave(X)", zsto,zout)
    119119
    120              if (clouds.eq.1) then
    121 
     120            if (clouds.eq.1) then
    122121             CALL histdef(nid_mth, "qnoy","nb tot noy" , "n/m2",
    123122     .                    iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    … …  
    149148     .                "ave(X)", zsto,zout)
    150149c --------------
    151 c ----- PRECIPITATIONS (precipitations cumulatives)
    152                CALL histdef(nid_mth,"prech4","Precip CH4","m",
     150c ----- ECHANGE GAZ SURF/ATM (evaporation)
     151               CALL histdef(nid_mth, "evapch4", "Evaporation CH4","m",
    153152     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    154153     .                "ave(X)", zsto,zout)
     154c --------------
     155c ----- PRECIPITATIONS (precipitations moyennes)
     156               CALL histdef(nid_mth,"prech4","Precip CH4","um/s",
     157     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     158     .                "ave(X)", zsto,zout)
    155159               CALL histdef(nid_mth,"prec2h6","Precip C2H6",
    156      .                "m",iim,jjmp1,nhori, 1,1,1, -99, 32,
     160     .                "um/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
    157161     .                "ave(X)", zsto,zout)
    158162               CALL histdef(nid_mth,"prec2h2","Precip C2H2",
    159      .                "m",iim,jjmp1,nhori, 1,1,1, -99, 32,
     163     .                "um/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
    160164     .                "ave(X)", zsto,zout)
    161165               CALL histdef(nid_mth,"prenoy","Precip NOY",
    162      .                "m",iim,jjmp1,nhori, 1,1,1, -99, 32,
     166     .                "um/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
    163167     .                "ave(X)", zsto,zout)
    164168               CALL histdef(nid_mth,"preaer","Precip AER",
    165      .                "m",iim,jjmp1,nhori, 1,1,1, -99, 32,
     169     .                "um/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
    166170     .                "ave(X)", zsto,zout)
    167171c --------------
    … …  
    177181     .              "ave(X)", zsto,zout)
    178182c --------------
     183c ----- Source/puits GLACE
     184               CALL histdef(nid_mth,"solch4", "dQ gl CH4",
     185     .              "m3/m3",iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     186     .              "ave(X)", zsto,zout)
     187               CALL histdef(nid_mth,"solc2h6", "dQ gl C2H6",
     188     .              "m3/m3",iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     189     .              "ave(X)", zsto,zout)
     190               CALL histdef(nid_mth,"solc2h2", "dQ gl C2H2",
     191     .              "m3/m3",iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     192     .              "ave(X)", zsto,zout)
     193c --------------
    179194c ----- RAYON DES GOUTTES
    180195               CALL histdef(nid_mth,"rcldbar", "rayon moyen goutte",
    181196     .                "m",iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    182197     .                "ave(X)", zsto,zout)
    183              endif
     198            endif
    184199          endif
    185200c --------------

  • trunk/LMDZ.TITAN/libf/phytitan/muphys3D.F

    r176 r474  
    107107         REAL pdq(ngrid,klev,nmicro)
    108108         REAL flxesp_i(ngrid,klev,3)    ! flx esp GLACE
    109          REAL solesp(ngrid,klev+1,3)    ! tx prod glace (puit/source)  ! ind klev+1 = surface ?
     109         REAL solesp(ngrid,klev,3)      ! tx prod glace (puit/source)
    110110         REAL tau_drop(ngrid,klev)
    111111         REAL tau_aer(ngrid,klev,nrad)
    … …  
    134134
    135135
    136          real ttq(ngrid,klev,nmicro,2)
    137          real tttq(ngrid,klev,nmicro,2)
     136c         real ttq(ngrid,klev,nmicro,2)
     137c         real tttq(ngrid,klev,nmicro,2)
    138138
    139139
    … …  
    351351           do j=1,klev
    352352             do i=1,nrad
     353!              solesp en m3/m3 pour passer en m3/m2 il faut faire :
     354!              (c1i(j,i)*dzb(j) -q(IHOR,j,i+2*nrad))
    353355               solesp(ihor,klev+1-j,1)=solesp(ihor,klev+1-j,1) +
    354356     &         (c1i(j,i)-q(IHOR,j,i+2*nrad)/dzb(j))

  • trunk/LMDZ.TITAN/libf/phytitan/phytrac.F

    r306 r474  
    7070c Local variables
    7171      REAL qaer0(klon,klev,nqmax)
     72      REAL prec(klon,5)
    7273
    7374c  ASTUCE POUR EVITER klon... EN ATTENDANT MIEUX
    … …  
    108109      integer      i,j,k,l,iq,ig0
    109110     
    110       REAL zprec(jjm+1,5),zsolesp(jjm+1,klev+1,3),
     111      REAL zprec(jjm+1,5),zsolesp(jjm+1,klev,3),
    111112     &     zflxesp_i(jjm+1,klev,3)
    112113      REAL ztau_drop(jjm+1,klev),ztau_aer(jjm+1,klev,nrad)
    … …  
    163164          tau_aer(:,:,:)  = 0.
    164165          solesp(:,:,:)   = 0.
    165           prec(:,:)       = 0.
     166          precip(:,:)     = 0.   ! c'est uniquement une sortie en um/s
     167c
     168          prec(:,:)       = 0.   ! c'est la variable temporaire des precipitions de la microfi
     169                                 ! prec est en m (metre precipitable)
    166170
    167171c-----------------------------------------------------------------------
    … …  
    507511         IF (clouds.eq.1) THEN
    508512           DO iq=1,3
    509              DO l=1,llm+1
    510                if (l.le.llm) flxesp_i(1,l,iq) = zflxesp_i(1,l,iq)
     513             DO l=1,llm
     514               flxesp_i(1,l,iq) = zflxesp_i(1,l,iq)
    511515               solesp(1,l,iq) = zsolesp(1,l,iq)
    512516               ig0 = 2
    513517               DO j=2,jjm
    514518                 DO i = 1, iim
    515                    if (l.le.llm) flxesp_i(ig0,l,iq)=zflxesp_i(1,l,iq)
    516                    solesp(ig0,l,iq) = zsolesp(1,l,iq)   
     519                   flxesp_i(ig0,l,iq)=zflxesp_i(j,l,iq)
     520                   solesp(ig0,l,iq) = zsolesp(j,l,iq)   
    517521                   ig0 = ig0 + 1
    518522                 ENDDO
    519523               ENDDO
    520                if (l.le.llm) flxesp_i(ig0,l,iq)=zflxesp_i(jjm+1,l,iq)
     524               flxesp_i(ig0,l,iq)=zflxesp_i(jjm+1,l,iq)
    521525               solesp(ig0,l,iq) = zsolesp(jjm+1,l,iq)
    522526             ENDDO
    … …  
    568572     &                pu,pv,pplev,pzlay,pzlev,
    569573     &                gaz1,gaz2,gaz3,
    570      &                ftsol,solesp,reservoir)
     574     &                ftsol,evapch4,reservoir)
    571575 
    572576c        ch4b=0.
    … …  
    778782c  CONDENSATION VIA MICROFI
    779783c----------------------
    780 c La microphysique avec nuages doit se faire obligatoirement en 3D.
     784c La microphysique avec nuages doit se faire obligatoirement en 3D.  (FAUX ACTUELLEMENT)
    781785c Rien n'empeche de faire la chimie en 2D. Cependant pour prendre en compte la
    782786c condensation due a la microfi (en 3D) on recalcule la tendance finale pour
    … …  
    836840
    837841c--------------------------------------------------
     842c  CONVERSION PRECIPITATION :
     843c  en microns/secondes
     844c--------------------------------------------------
     845        precip = prec * 1.e6 / ptimestep
     846
     847
     848c--------------------------------------------------
    838849c CALCUL DU FLUX DE CHALEUR LATENTE D'EVAPORATION
    839850c DU METHANE
    … …  
    848859     &     (7.08*ftm**0.354+10.95*1.1e-2*ftm**0.456)
    849860     &     /mch4
    850            ! solcxhy en m3/m2 {ok}
     861           ! evapch4 en m3/m2 {ok}
    851862           ! 425 en kg/m3
    852863           ! Lv en J/kg       {ok}
    853864           ! ptimestep en s   {ok}
    854            fclat(i)=(solesp(i,klev+1,1)*Lvch4*rhoi_ch4)   ! en J/m2/s
     865           fclat(i)=(evapch4(i)*Lvch4*rhoi_ch4)   ! en J/m2/s
    855866         ENDDO
    856867       ENDIF
    … …  
    873884     &           ((qaer(i,j,iq+2*nrad)/qaer(i,j,iq+nrad)+
    874885     &           qaer(i,j,iq+3*nrad)/qaer(i,j,iq+nrad) +
    875      &           v_e(iq))*0.75/3.1415926353)**(1./3.)
     886     &           v_e(iq))*0.75/RPI)**(1./3.)
    876887               ENDIF
    877888             ENDDO
    … …  
    908919     &         print*, 'PB AVEC XFRAC:', i,j,xgsn,vcl
    909920               rmcloud(i,j)=          ! rayon moyen des gouttes
    910      &         (vcl/nuc*0.75/3.1415926353)**(1./3.)
     921     &         (vcl/nuc*0.75/RPI)**(1./3.)
    911922               xfrac(i,j,1)=xgsn/vcl         ! fraction volumique noyau/goutte
    912923               xfrac(i,j,2)=xmsn/vcl         ! fraction volumique CH4/goutte
    … …  
    918929               ncount(i,j) = ncount(i,j)+1
    919930             ENDIF
    920              socccld=socccld+3.1415926*(rmcloud(i,j)**2.)*nuc
     931             socccld=socccld+RPI*(rmcloud(i,j)**2.)*nuc
    921932             occcld(i,j)=socccld
    922933           ENDDO
    … …  
    940951c         12   31.622778       100000.00
    941952c
     953c        mise a zero de occld_m
     954         occcld_m=0.
    942955         DO i=1,klon
    943956           DO j=1,klev

  • trunk/LMDZ.TITAN/libf/phytitan/snuages3D.F

    r175 r474  
    240240 
    241241c        En theorie, les diceX sont NEGATIF (en sedimentant on ne fait que perdre de la glace)
    242 c        Le precipitations sont comptées positivement.
    243          precip(ihor,1)=-dice1/rhoi_ch4     ! m3/m2 = m :)
    244          precip(ihor,2)=-dice2/rhoi_c2h6
    245          precip(ihor,3)=-dice3/rhoi_c2h2
    246          precip(ihor,4)=-dice4/rhol
     242c        Les precipitations sont comptees positivement. (ET ON NE PREND QUE DES VALEURS POSITIVES !)
     243
     244         precip(ihor,1)=AMAX1(-dice1/rhoi_ch4,0)     ! m3/m2 = m :)
     245         precip(ihor,2)=AMAX1(-dice2/rhoi_c2h6,0)
     246         precip(ihor,3)=AMAX1(-dice3/rhoi_c2h2,0)
     247         precip(ihor,4)=AMAX1(-dice4/rhol,0)
    247248
    248249         do i=1,nz

  • trunk/LMDZ.TITAN/libf/phytitan/sources.F

    r175 r474  
    33     $                      pplev,pzlay,pzlev,
    44     $                      gaz1,gaz2,gaz3,
    5      $                      ptsrf,solesp,reserv)
     5     $                      ptsrf,evapch4,reserv)
    66
    77c=======================================================================
    … …  
    3737         REAL gaz1(ngrid,nlay),gaz2(ngrid,nlay),gaz3(ngrid,nlay)
    3838         REAL ptsrf(ngrid)
    39          REAL solesp(ngrid,klev+1,3)
     39         REAL evapch4(ngrid)
    4040c
    4141c   local:
    … …  
    5151         REAL prodc2h6,prodc2h2
    5252         real reserv(ngrid),restemp,drestemp
    53          REAL evapch4
     53         REAL zevapch4
    5454
    5555         real umin
    5656         data umin/1.e-12/
    5757         save umin
    58 
     58c
    5959c
    6060c-----------------------------------------------------------------------
    … …  
    8484         zgz3 = gaz3
    8585
     86         evapch4 = 0.
     87
    8688c-----------------------------------------------------------------------
    8789c     2. calcul de  cd :
    … …  
    103105*   Conditions CH4
    104106         DO ig=1,ngrid
    105            evapch4=0.
     107           zevapch4=0.
    106108           restemp=0.
    107109           gg=RG*RA**2/(RA+pzlay(ig,1))**2.
    … …  
    110112           ws=sqrt(pu(ig)**2.+pv(ig)**2.)*(10./pzlay(ig,1))**0.2
    111113           ch=1.5*sqrt(zcdv(ig))
    112 c           write(100,'(I4,3(ES24.17,1X))')
    113 c     &     ig,sqrt(pu(ig)**2.+pv(ig)**2.),
    114 c     &     (10./pzlay(ig,1))**0.2,ws
    115114           call ch4sat(ptsrf(ig),pplev(ig,1),qch4) ! qch4=kg/kg
    116115           qch4=qch4*0.50 ! ici on impose 50% d'humidité au sol
    … …  
    123122             flux=flux*0.1 ! fraction occupée par les lacs
    124123           endif
    125 c           write(101,'(I4,5(ES24.17,1X))')
    126 c     &     ig,reserv(ig),ws,ch,zrho,(pzlev(ig,1+1)-pzlev(ig,1))
    127 c           flux=flux/zrho/(pzlev(ig,1+1)-pzlev(ig,1)) ! dx/ds
    128124
    129125           zmem=zgz1(ig,1)
    … …  
    146142           IF (restemp.ge.0.) THEN
    147143             reserv(ig) = reserv(ig) + drestemp
    148              evapch4    = evapch4    + drestemp
     144             zevapch4   = zevapch4   + drestemp
    149145           ELSE
    150146*          Il n'y a pas suffisamment de méthane; on re-évalue le flux d'évaporation
    151147*          Quelle nouvelle concentration zgz1(ig,1) atteint-on en évaporant tout ?
    152              zgz1(ig,1)=reserv(ig)/(xmair*(pzlev(ig,1+1)-pzlev(ig,1))
     148             zgz1(ig,1)= reserv(ig)/(xmair*(pzlev(ig,1+1)-pzlev(ig,1))
    153149     &                 *16./xmuair/425.)+zmem
    154              evapch4=evapch4-reserv(ig)
     150             zevapch4  = zevapch4-reserv(ig)
    155151
    156152             if(reserv(ig).eq.0.)
    … …  
    160156           ENDIF
    161157c         
    162            solesp(ig,klev+1,1)=solesp(ig,klev+1,1)+evapch4  ! <0 si volume évaporé (m3/m2)
    163 
    164 c         write(498,'(I3,4(ES24.17,1X))') ig,gaz1(ig,1),zgz1(ig,1),flux
     158           evapch4(ig)=evapch4(ig)+zevapch4  ! evapch4 doit etre < 0
     159
    165160         ENDDO
    166 c         write(498,*) ""
    167 c         write(101,*) ""
     161
    168162*   Conditions C2H6
    169163

  • trunk/LMDZ.TITAN/libf/phytitan/write_histday.h

    r306 r474  
    114114
    115115             if (clouds.eq.1) then
    116 
    117116c    -------   NB NOY TOT
    118117               do i=1,klon
    … …  
    170169     .                        iim*jjmp1,ndex2d)
    171170c --------------
     171c ----- ECHANGE GAZ SURF/ATM (evaporation)
     172       CALL gr_fi_ecrit(1, klon,iim,jjmp1,evapch4,zx_tmp_2d)
     173       CALL histwrite(nid_day,"evapch4",itau_w,zx_tmp_2d,
     174     .                        iim*jjmp1,ndex2d)
     175c --------------
    172176c ----- PRECIPITATIONS
    173177c       -----  CH4
    174        CALL gr_fi_ecrit(1, klon,iim,jjmp1,prec(:,1),zx_tmp_2d)
     178       CALL gr_fi_ecrit(1, klon,iim,jjmp1,precip(1:klon,1),zx_tmp_2d)
    175179       CALL histwrite(nid_day,"prech4",itau_w,zx_tmp_2d,
    176180     .                        iim*jjmp1,ndex2d)
    177181c       -----  C2H6
    178        CALL gr_fi_ecrit(1, klon,iim,jjmp1,prec(:,2),zx_tmp_2d)
     182       CALL gr_fi_ecrit(1, klon,iim,jjmp1,precip(1:klon,2),zx_tmp_2d)
    179183       CALL histwrite(nid_day,"prec2h6",itau_w,zx_tmp_2d,
    180184     .                        iim*jjmp1,ndex2d)
    181185c       -----  C2H2
    182        CALL gr_fi_ecrit(1, klon,iim,jjmp1,prec(:,3),zx_tmp_2d)
     186       CALL gr_fi_ecrit(1, klon,iim,jjmp1,precip(1:klon,3),zx_tmp_2d)
    183187       CALL histwrite(nid_day,"prec2h2",itau_w,zx_tmp_2d,
    184188     .                        iim*jjmp1,ndex2d)
    … …  
    186190c --------------
    187191c ----- FLUX GLACE
    188        CALL gr_fi_ecrit(klev,klon,iim,jjmp1,flxesp_i(1,1,1),zx_tmp_3d)
     192       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,
     193     .              flxesp_i(1:klon,1:klev,1),zx_tmp_3d)
    189194       CALL histwrite(nid_day,"flxgl1", itau_w, zx_tmp_3d,
    190195     .                                   iim*jjmp1*klev,ndex3d)
    191        CALL gr_fi_ecrit(klev,klon,iim,jjmp1,flxesp_i(1,1,2),zx_tmp_3d)
     196       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,
     197     .              flxesp_i(1:klon,1:klev,2),zx_tmp_3d)
    192198       CALL histwrite(nid_day,"flxgl2", itau_w, zx_tmp_3d,
    193199     .                                   iim*jjmp1*klev,ndex3d)
    194        CALL gr_fi_ecrit(klev,klon,iim,jjmp1,flxesp_i(1,1,3),zx_tmp_3d)
     200       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,
     201     .              flxesp_i(1:klon,1:klev,3),zx_tmp_3d)
    195202       CALL histwrite(nid_day,"flxgl3", itau_w, zx_tmp_3d,
    196203     .                                   iim*jjmp1*klev,ndex3d)

  • trunk/LMDZ.TITAN/libf/phytitan/write_histmth.h

    r306 r474  
    9797          if (microfi.ge.1) then
    9898c          DO iq=1,nmicro
    99 c      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qaer(1,1,iq), zx_tmp_3d)
     99c      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,
     100c    .                  qaer(1:klon,1:klev,iq), zx_tmp_3d)
    100101c      CALL histwrite(nid_mth,tname(iq),itau_w,zx_tmp_3d,
    101102c    .                                   iim*jjmp1*klev,ndex3d)
    … …  
    112113
    113114             if (clouds.eq.1) then
    114 
    115115c    -------   NB NOY TOT
    116116               do i=1,klon
    … …  
    168168     .                        iim*jjmp1,ndex2d)
    169169c --------------
     170c ----- ECHANGE GAZ SURF/ATM (evaporation)
     171       CALL gr_fi_ecrit(1, klon,iim,jjmp1,evapch4,zx_tmp_2d)
     172       CALL histwrite(nid_mth,"evapch4",itau_w,zx_tmp_2d,
     173     .                        iim*jjmp1,ndex2d)
     174c --------------
    170175c ----- PRECIPITATIONS
    171176c       -----  CH4
    172        CALL gr_fi_ecrit(1, klon,iim,jjmp1,prec(:,1),zx_tmp_2d)
     177       CALL gr_fi_ecrit(1, klon,iim,jjmp1,precip(1:klon,1),zx_tmp_2d)
    173178       CALL histwrite(nid_mth,"prech4",itau_w,zx_tmp_2d,
    174179     .                        iim*jjmp1,ndex2d)
    175180c       -----  C2H6
    176        CALL gr_fi_ecrit(1, klon,iim,jjmp1,prec(:,2),zx_tmp_2d)
     181       CALL gr_fi_ecrit(1, klon,iim,jjmp1,precip(1:klon,2),zx_tmp_2d)
    177182       CALL histwrite(nid_mth,"prec2h6",itau_w,zx_tmp_2d,
    178183     .                        iim*jjmp1,ndex2d)
    179184c       -----  C2H2
    180        CALL gr_fi_ecrit(1, klon,iim,jjmp1,prec(:,3),zx_tmp_2d)
     185       CALL gr_fi_ecrit(1, klon,iim,jjmp1,precip(1:klon,3),zx_tmp_2d)
    181186       CALL histwrite(nid_mth,"prec2h2",itau_w,zx_tmp_2d,
    182187     .                        iim*jjmp1,ndex2d)
    183188c
     189c       -----  NOY
     190       CALL gr_fi_ecrit(1, klon,iim,jjmp1,precip(1:klon,4),zx_tmp_2d)
     191       CALL histwrite(nid_mth,"prenoy",itau_w,zx_tmp_2d,
     192     .                        iim*jjmp1,ndex2d)
     193c       -----  AER
     194       CALL gr_fi_ecrit(1, klon,iim,jjmp1,precip(1:klon,5),zx_tmp_2d)
     195       CALL histwrite(nid_mth,"preaer",itau_w,zx_tmp_2d,
     196     .                        iim*jjmp1,ndex2d)
    184197c --------------
    185198c ----- FLUX GLACE
    186        CALL gr_fi_ecrit(klev,klon,iim,jjmp1,flxesp_i(1,1,1),zx_tmp_3d)
     199c       -----  CH4
     200       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,
     201     .                  flxesp_i(1:klon,1:klev,1),zx_tmp_3d)
    187202       CALL histwrite(nid_mth,"flxgl1", itau_w, zx_tmp_3d,
    188203     .                                   iim*jjmp1*klev,ndex3d)
    189        CALL gr_fi_ecrit(klev,klon,iim,jjmp1,flxesp_i(1,1,2),zx_tmp_3d)
     204c       -----  C2H6
     205       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,
     206     .                  flxesp_i(1:klon,1:klev,2),zx_tmp_3d)
    190207       CALL histwrite(nid_mth,"flxgl2", itau_w, zx_tmp_3d,
    191208     .                                   iim*jjmp1*klev,ndex3d)
    192        CALL gr_fi_ecrit(klev,klon,iim,jjmp1,flxesp_i(1,1,3),zx_tmp_3d)
     209c       -----  C2H2
     210       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,
     211     .                  flxesp_i(1:klon,1:klev,3),zx_tmp_3d)
    193212       CALL histwrite(nid_mth,"flxgl3", itau_w, zx_tmp_3d,
     213     .                                   iim*jjmp1*klev,ndex3d)
     214c --------------
     215c ----- Source/puits GLACE
     216c       -----  CH4
     217       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,
     218     .                          solesp(1:klon,1:klev,1),zx_tmp_3d)
     219       CALL histwrite(nid_mth,"solch4", itau_w, zx_tmp_3d,
     220     .                                   iim*jjmp1*klev,ndex3d)
     221c       -----  C2H6
     222       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,
     223     .                          solesp(1:klon,1:klev,2),zx_tmp_3d)
     224       CALL histwrite(nid_mth,"solc2h6", itau_w, zx_tmp_3d,
     225     .                                   iim*jjmp1*klev,ndex3d)
     226c       -----  C2H2
     227       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,
     228     .                          solesp(1:klon,1:klev,3),zx_tmp_3d)
     229       CALL histwrite(nid_mth,"solc2h2", itau_w, zx_tmp_3d,
    194230     .                                   iim*jjmp1*klev,ndex3d)
    195231c
    … …  
    414450             write(str1,'(i2.2)') k
    415451             zx_tmp_fi3d(1:klon,1:klev)=occcld_m(1:klon,1:klev,k)
     452             CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    416453             CALL histwrite(nid_mth,"occcld"//str1,itau_w,zx_tmp_3d,
    417454     .                                   iim*jjmp1*klev,ndex3d)
Note: See TracChangeset for help on using the changeset viewer.