Changeset 110


Ignore:
Timestamp:
Apr 15, 2011, 3:44:44 PM (14 years ago)
Author:
slebonnois
Message:

SL: corrections de bugs suite a compilations venus et titan de la version 109.

Location:
trunk
Files:
14 edited

Legend:

Unmodified
Added
Removed

  • trunk/chantiers/commit_importants.log

    r109 r110  
    802802  Ajout dissip_horiz.tex (et .pdf) dans la documentation.
    803803
    804 *********************
    805 **** commit_v109 ****
    806 *********************
     804****************************
     805**** commit_v109 (+110) ****
     806****************************
    807807
    808808** Discretisation verticale

  • trunk/libf/dyn3d/disvert_noterre.F

    r109 r110  
    199199c     Specifique.  Ce choix est spécifié ici ET dans exner_hyb.F
    200200
    201       DO l = 1, llm
     201      DO l = 1, llm-1
    202202       aps(l) =  0.5 *( ap(l) +ap(l+1))
    203203       bps(l) =  0.5 *( bp(l) +bp(l+1))
    204204      ENDDO
    205205     
    206       if (hybrid) then
     206c     if (hybrid) then
    207207         aps(llm) = aps(llm-1)**2 / aps(llm-2)
    208208         bps(llm) = 0.5*(bp(llm) + bp(llm+1))
    209       else
    210          bps(llm) = bps(llm-1)**2 / bps(llm-2)
    211          aps(llm) = 0.
    212       end if
     209c     else
     210c        bps(llm) = bps(llm-1)**2 / bps(llm-2)
     211c        aps(llm) = 0.
     212c     end if
    213213
    214214      PRINT *,' BPs '

  • trunk/libf/dyn3d/disvert_terre.F90

    r109 r110  
    140140  write(lunout, *) presnivs
    141141
    142 END SUBROUTINE disvert
     142END SUBROUTINE disvert_terre

  • trunk/libf/dyn3d/inidissip.F

    r108 r110  
    185185
    186186         write(*,*) 'Dissipation : '
    187          write(*,*) 'Multiplication de la dissipation en altitude :',
     187         write(*,*) 'Multiplication de la dissipation en altitude :'
    188188         write(*,*) '  dissip_fac_mid =', dissip_fac_mid
    189189

  • trunk/libf/dyn3d/leapfrog.F

    r109 r110  
    470470c      -------------------
    471471         IF (ok_strato) THEN
    472            CALL top_bound( vcov,ucov,teta,masse,dutop,dvtop,dtetatop)
     472           CALL top_bound( vcov,ucov,teta,phi,masse,
     473     $                     dutop,dvtop,dtetatop)
    473474c dqtop=0, dptop=0
    474475           CALL addfi( dtphys, leapf, forward   ,
    … …  
    516517        ! Sponge layer (if any)
    517518        IF (ok_strato) THEN
    518           CALL top_bound(vcov,ucov,teta,masse,dutop,dvtop,dtetatop)
     519          CALL top_bound(vcov,ucov,teta,phi,
     520     $                   masse,dutop,dvtop,dtetatop)
    519521c dqtop=0, dptop=0
    520522          CALL addfi( dtvr, leapf, forward   ,

  • trunk/libf/dyn3d/top_bound.F

    r108 r110  
    1       SUBROUTINE top_bound( vcov,ucov,teta,masse, du,dv,dh )
     1      SUBROUTINE top_bound( vcov,ucov,teta,phi,masse, du,dv,dh )
    22      IMPLICIT NONE
    33c
    … …  
    3535
    3636      REAL ucov(iip1,jjp1,llm),vcov(iip1,jjm,llm),teta(iip1,jjp1,llm)
     37      REAL phi(iip1,jjp1,llm)                  ! geopotentiel
    3738      REAL masse(iip1,jjp1,llm)
    3839      REAL dv(iip1,jjm,llm),du(iip1,jjp1,llm),dh(iip1,jjp1,llm)
    … …  
    7778         print*,'TOP_BOUND mode',mode_top_bound
    7879         print*,'Coeffs pour la couche eponge a l equateur'
    79          print*,'p (Pa)  z(km)  tau (s)   dt*rdamp'
     80         print*,'p (Pa)  z(km)  tau (s)'
    8081         do l=1,llm
    8182           if (rdamp(l).ne.0.) then
    8283            zkm        = phi(iip1/2,jjp1/2,l)/(1000*g)
    83           print*,presnivs(l),zkm,
    84      .          1./rdamp(l),
    85      .          dt*rdamp(l)
     84          print*,presnivs(l),zkm,1./rdamp(l)
    8685           endif
    8786         enddo

  • trunk/libf/dyn3dpar/disvert_noterre.F

    r109 r110  
    199199c     Specifique.  Ce choix est spécifié ici ET dans exner_hyb.F
    200200
    201       DO l = 1, llm
     201      DO l = 1, llm-1
    202202       aps(l) =  0.5 *( ap(l) +ap(l+1))
    203203       bps(l) =  0.5 *( bp(l) +bp(l+1))
    204204      ENDDO
    205205     
    206       if (hybrid) then
     206c     if (hybrid) then
    207207         aps(llm) = aps(llm-1)**2 / aps(llm-2)
    208208         bps(llm) = 0.5*(bp(llm) + bp(llm+1))
    209       else
    210          bps(llm) = bps(llm-1)**2 / bps(llm-2)
    211          aps(llm) = 0.
    212       end if
     209c     else
     210c        bps(llm) = bps(llm-1)**2 / bps(llm-2)
     211c        aps(llm) = 0.
     212c     end if
    213213
    214214      PRINT *,' BPs '

  • trunk/libf/dyn3dpar/disvert_terre.F90

    r109 r110  
    140140  write(lunout, *) presnivs
    141141
    142 END SUBROUTINE disvert
     142END SUBROUTINE disvert_terre

  • trunk/libf/dyn3dpar/inidissip.F

    r109 r110  
    185185
    186186         write(*,*) 'Dissipation : '
    187          write(*,*) 'Multiplication de la dissipation en altitude :',
     187         write(*,*) 'Multiplication de la dissipation en altitude :'
    188188         write(*,*) '  dissip_fac_mid =', dissip_fac_mid
    189189

  • trunk/libf/dyn3dpar/leapfrog_p.F

    r109 r110  
    947947c      -------------------
    948948         IF (ok_strato) THEN
    949            CALL top_bound_p( vcov,ucov,teta,masse,dutop,dvtop,dtetatop)
     949           CALL top_bound_p( vcov,ucov,teta,phi,masse,
     950     $                       dutop,dvtop,dtetatop)
    950951           CALL addfi_p( dtphys, leapf, forward   ,
    951952     $                  ucov, vcov, teta , q   ,ps ,
    … …  
    10691070        ! Sponge layer (if any)
    10701071        IF (ok_strato) THEN
    1071           ! set dufi,dvfi,... to zero
     1072          ! set dutop,dvtop,... to zero
    10721073          ijb=ij_begin
    10731074          ije=ij_end
    10741075!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    10751076          do l=1,llm
    1076             dufi(ijb:ije,l)=0
    1077             dtetafi(ijb:ije,l)=0
    1078             dqfi(ijb:ije,l,1:nqtot)=0
     1077            dutop(ijb:ije,l)=0
     1078            dtetatop(ijb:ije,l)=0
     1079            dqtop(ijb:ije,l,1:nqtot)=0
    10791080          enddo
    10801081!$OMP END DO
    10811082!$OMP SINGLE
    1082           dpfi(ijb:ije)=0
     1083          dptop(ijb:ije)=0
    10831084!$OMP END SINGLE
    10841085          ijb=ij_begin
    … …  
    10871088!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    10881089          do l=1,llm
    1089             dvfi(ijb:ije,l)=0
     1090            dvtop(ijb:ije,l)=0
    10901091          enddo
    10911092!$OMP END DO
    10921093
    1093           CALL top_bound_p(vcov,ucov,teta,masse,dufi,dvfi,dtetafi)
     1094          CALL top_bound_p(vcov,ucov,teta,phi,masse,
     1095     $                     dutop,dvtop,dtetatop)
    10941096          CALL addfi_p( dtvr, leapf, forward   ,
    10951097     $                  ucov, vcov, teta , q   ,ps ,
    1096      $                 dufi, dvfi, dtetafi , dqfi ,dpfi  )
     1098     $                 dutop, dvtop, dtetatop , dqtop ,dptop  )
    10971099!$OMP BARRIER
    10981100        ENDIF ! of IF (ok_strato)

  • trunk/libf/dyn3dpar/top_bound_p.F

    r108 r110  
    1       SUBROUTINE top_bound_p( vcov,ucov,teta,masse, du,dv,dh )
     1      SUBROUTINE top_bound_p( vcov,ucov,teta,phi,masse, du,dv,dh )
    22      USE parallel
    33      IMPLICIT NONE
    … …  
    3535
    3636      REAL ucov(iip1,jjp1,llm),vcov(iip1,jjm,llm),teta(iip1,jjp1,llm)
     37      REAL phi(iip1,jjp1,llm)                  ! geopotentiel
    3738      REAL masse(iip1,jjp1,llm)
    3839      REAL dv(iip1,jjm,llm),du(iip1,jjp1,llm),dh(iip1,jjp1,llm)
    … …  
    7273         print*,'TOP_BOUND mode',mode_top_bound
    7374         print*,'Coeffs pour la couche eponge a l equateur'
    74          print*,'p (Pa)  z(km)  tau (s)   dt*rdamp'
     75         print*,'p (Pa)  z(km)  tau (s)'
    7576         do l=1,llm
    7677           if (rdamp(l).ne.0.) then
    7778            zkm        = phi(iip1/2,jjp1/2,l)/(1000*g)
    78           print*,presnivs(l),zkm,
    79      .          1./rdamp(l),
    80      .          dt*rdamp(l)
     79          print*,presnivs(l),zkm,1./rdamp(l)
    8180           endif
    8281         enddo

  • trunk/libf/phytitan/write_histday.h

    r106 r110  
    105105           DO iq=1,nmicro
    106106       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qaer(1,1,iq), zx_tmp_3d)
    107        z4_tmp_3d=zx_tmp_3d
    108        CALL histwrite(nid_day,tname(iq),itau_w,z4_tmp_3d,
     107       CALL histwrite(nid_day,tname(iq),itau_w,zx_tmp_3d,
    109108     .                                   iim*jjmp1*klev,ndex3d)
    110109           ENDDO
    … …  
    113112           DO iq=nmicro+1,nqmax
    114113       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,tr_seri(1,1,iq),zx_tmp_3d)
    115        z4_tmp_3d=zx_tmp_3d
    116        CALL histwrite(nid_day,tname(iq),itau_w,z4_tmp_3d,
     114       CALL histwrite(nid_day,tname(iq),itau_w,zx_tmp_3d,
    117115     .                                   iim*jjmp1*klev,ndex3d)
    118116           ENDDO
    … …  
    156154      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
    157155      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,
     156      CALL histwrite(nid_day,"thv"//str1,itau_w,zx_tmp_3d,
    160157     .                                   iim*jjmp1*klev,ndex3d)
    161158       enddo      ! fin boucle NSPECV
    … …  
    177174      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
    178175      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,
     176      CALL histwrite(nid_day,"khv"//str1,itau_w,zx_tmp_3d,
    181177     .                                   iim*jjmp1*klev,ndex3d)
    182178       enddo      ! fin boucle NSPECV
    … …  
    191187      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
    192188      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,
     189      CALL histwrite(nid_day,"tgv"//str1,itau_w,zx_tmp_3d,
    195190     .                                   iim*jjmp1*klev,ndex3d)
    196191       enddo      ! fin boucle NSPECV
    … …  
    212207      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
    213208      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,
     209      CALL histwrite(nid_day,"kgv"//str1,itau_w,zx_tmp_3d,
    216210     .                                   iim*jjmp1*klev,ndex3d)
    217211       enddo      ! fin boucle NSPECV
    … …  
    226220      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
    227221      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,
     222      CALL histwrite(nid_day,"thi"//str1,itau_w,zx_tmp_3d,
    230223     .                                   iim*jjmp1*klev,ndex3d)
    231224       enddo      ! fin boucle NSPECI
    … …  
    247240      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
    248241      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,
     242      CALL histwrite(nid_day,"khi"//str1,itau_w,zx_tmp_3d,
    251243     .                                   iim*jjmp1*klev,ndex3d)
    252244       enddo      ! fin boucle NSPECI
    … …  
    261253      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
    262254      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,
     255      CALL histwrite(nid_day,"tgi"//str1,itau_w,zx_tmp_3d,
    265256     .                                   iim*jjmp1*klev,ndex3d)
    266257       enddo      ! fin boucle NSPECI
    … …  
    282273      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
    283274      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,
     275      CALL histwrite(nid_day,"kgi"//str1,itau_w,zx_tmp_3d,
    286276     .                                   iim*jjmp1*klev,ndex3d)
    287277       enddo      ! fin boucle NSPECI

  • trunk/libf/phytitan/write_histins.h

    r106 r110  
    154154      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
    155155      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,
     156      CALL histwrite(nid_ins,"thv"//str1,itau_w,zx_tmp_3d,
    158157     .                                   iim*jjmp1*klev,ndex3d)
    159158       enddo      ! fin boucle NSPECV
    … …  
    175174      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
    176175      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,
     176      CALL histwrite(nid_ins,"khv"//str1,itau_w,zx_tmp_3d,
    179177     .                                   iim*jjmp1*klev,ndex3d)
    180178       enddo      ! fin boucle NSPECV
    … …  
    189187      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
    190188      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,
     189      CALL histwrite(nid_ins,"tgv"//str1,itau_w,zx_tmp_3d,
    193190     .                                   iim*jjmp1*klev,ndex3d)
    194191       enddo      ! fin boucle NSPECV
    … …  
    210207      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
    211208      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,
     209      CALL histwrite(nid_ins,"kgv"//str1,itau_w,zx_tmp_3d,
    214210     .                                   iim*jjmp1*klev,ndex3d)
    215211       enddo      ! fin boucle NSPECV
    … …  
    224220      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
    225221      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,
     222      CALL histwrite(nid_ins,"thi"//str1,itau_w,zx_tmp_3d,
    228223     .                                   iim*jjmp1*klev,ndex3d)
    229224       enddo      ! fin boucle NSPECI
    … …  
    245240      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
    246241      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,
     242      CALL histwrite(nid_ins,"khi"//str1,itau_w,zx_tmp_3d,
    249243     .                                   iim*jjmp1*klev,ndex3d)
    250244       enddo      ! fin boucle NSPECI
    … …  
    259253      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
    260254      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,
     255      CALL histwrite(nid_ins,"tgi"//str1,itau_w,zx_tmp_3d,
    263256     .                                   iim*jjmp1*klev,ndex3d)
    264257       enddo      ! fin boucle NSPECI
    … …  
    280273      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
    281274      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,
     275      CALL histwrite(nid_ins,"kgi"//str1,itau_w,zx_tmp_3d,
    284276     .                                   iim*jjmp1*klev,ndex3d)
    285277       enddo      ! fin boucle NSPECI

  • trunk/libf/phytitan/write_histmth.h

    r106 r110  
    4949c
    5050c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
    51 c     z4_tmp_2d=zx_tmp_2d
    52 c     CALL histwrite(nid_mth,"cdragh",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     51c     CALL histwrite(nid_mth,"cdragh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    5352c
    5453c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
    55 c     z4_tmp_2d=zx_tmp_2d
    56 c     CALL histwrite(nid_mth,"cdragm",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
     54c     CALL histwrite(nid_mth,"cdragm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    5755c
    5856      ENDIF !lev_histmth.GE.1
    … …  
    8886c
    8987c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ycoefh, zx_tmp_3d)
    90 c     z4_tmp_3d=zx_tmp_3d
    91 c     CALL histwrite(nid_mth,"Kz",itau_w,z4_tmp_3d,
     88c     CALL histwrite(nid_mth,"Kz",itau_w,zx_tmp_3d,
    9289c    .                                   iim*jjmp1*klev,ndex3d)
    9390c
    … …  
    159156      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
    160157      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    161       z4_tmp_3d=zx_tmp_3d
    162       CALL histwrite(nid_day,"thv"//str1,itau_w,z4_tmp_3d,
     158      CALL histwrite(nid_day,"thv"//str1,itau_w,zx_tmp_3d,
    163159     .                                   iim*jjmp1*klev,ndex3d)
    164160       enddo      ! fin boucle NSPECV
    … …  
    180176      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
    181177      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    182       z4_tmp_3d=zx_tmp_3d
    183       CALL histwrite(nid_day,"khv"//str1,itau_w,z4_tmp_3d,
     178      CALL histwrite(nid_day,"khv"//str1,itau_w,zx_tmp_3d,
    184179     .                                   iim*jjmp1*klev,ndex3d)
    185180       enddo      ! fin boucle NSPECV
    … …  
    194189      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
    195190      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,
     191      CALL histwrite(nid_day,"tgv"//str1,itau_w,zx_tmp_3d,
    198192     .                                   iim*jjmp1*klev,ndex3d)
    199193       enddo      ! fin boucle NSPECV
    … …  
    215209      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
    216210      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,
     211      CALL histwrite(nid_day,"kgv"//str1,itau_w,zx_tmp_3d,
    219212     .                                   iim*jjmp1*klev,ndex3d)
    220213       enddo      ! fin boucle NSPECV
    … …  
    229222      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
    230223      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,
     224      CALL histwrite(nid_day,"thi"//str1,itau_w,zx_tmp_3d,
    233225     .                                   iim*jjmp1*klev,ndex3d)
    234226       enddo      ! fin boucle NSPECI
    … …  
    250242      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
    251243      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,
     244      CALL histwrite(nid_day,"khi"//str1,itau_w,zx_tmp_3d,
    254245     .                                   iim*jjmp1*klev,ndex3d)
    255246       enddo      ! fin boucle NSPECI
    … …  
    264255      zx_tmp_fi3d(1:klon,1:klev)=t_tauhvd(1:klon,1:klev)
    265256      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,
     257      CALL histwrite(nid_day,"tgi"//str1,itau_w,zx_tmp_3d,
    268258     .                                   iim*jjmp1*klev,ndex3d)
    269259       enddo      ! fin boucle NSPECI
    … …  
    285275      zx_tmp_fi3d(1:klon,1:klev)=t_khvd(1:klon,1:klev)
    286276      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,
     277      CALL histwrite(nid_day,"kgi"//str1,itau_w,zx_tmp_3d,
    289278     .                                   iim*jjmp1*klev,ndex3d)
    290279       enddo      ! fin boucle NSPECI
Note: See TracChangeset for help on using the changeset viewer.