/[lmdze]/trunk/libf/phylmd/phytrac.f90
ViewVC logotype

Diff of /trunk/libf/phylmd/phytrac.f90

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 14 by guez, Fri Jul 25 19:59:34 2008 UTC revision 15 by guez, Fri Aug 1 15:24:12 2008 UTC
# Line 509  contains Line 509  contains
509        itau_w = itau_phy + itap        itau_w = itau_phy + itap
510    
511        CALL gr_fi_ecrit(1, klon, iim, jjm+1, pphis, zx_tmp_2d)        CALL gr_fi_ecrit(1, klon, iim, jjm+1, pphis, zx_tmp_2d)
512        CALL histwrite(nid_tra, "phis", itau_w, zx_tmp_2d, iim*(jjm+1), ndex2d)        CALL histwrite(nid_tra, "phis", itau_w, zx_tmp_2d)
513    
514        CALL gr_fi_ecrit(1, klon, iim, jjm+1, airephy, zx_tmp_2d)              CALL gr_fi_ecrit(1, klon, iim, jjm+1, airephy, zx_tmp_2d)      
515        CALL histwrite(nid_tra, "aire", itau_w, zx_tmp_2d, iim*(jjm+1), ndex2d)        CALL histwrite(nid_tra, "aire", itau_w, zx_tmp_2d)
516    
517        CALL gr_fi_ecrit(llm, klon, iim, jjm+1, zmasse, zx_tmp_3d)              CALL gr_fi_ecrit(llm, klon, iim, jjm+1, zmasse, zx_tmp_3d)      
518        CALL histwrite(nid_tra, "zmasse", itau_w, zx_tmp_3d, iim*(jjm+1)*llm, &        CALL histwrite(nid_tra, "zmasse", itau_w, zx_tmp_3d)
            ndex3d)  
519    
520        DO it=1, nqmax        DO it=1, nqmax
521           CALL gr_fi_ecrit(llm, klon, iim, jjm+1, tr_seri(1, 1, it), zx_tmp_3d)           CALL gr_fi_ecrit(llm, klon, iim, jjm+1, tr_seri(1, 1, it), zx_tmp_3d)
522           CALL histwrite(nid_tra, tnom(it+2), itau_w, zx_tmp_3d, &           CALL histwrite(nid_tra, tnom(it+2), itau_w, zx_tmp_3d)
               iim*(jjm+1)*llm, ndex3d)  
523           if (lessivage) THEN           if (lessivage) THEN
524              CALL gr_fi_ecrit(llm, klon, iim, jjm+1, flestottr(1, 1, it), &              CALL gr_fi_ecrit(llm, klon, iim, jjm+1, flestottr(1, 1, it), &
525                   zx_tmp_3d)                   zx_tmp_3d)
526              CALL histwrite(nid_tra, "fl"//tnom(it+2), itau_w, zx_tmp_3d, &              CALL histwrite(nid_tra, "fl"//tnom(it+2), itau_w, zx_tmp_3d)
                  iim*(jjm+1)*llm, ndex3d)  
527           endif           endif
528    
529           CALL gr_fi_ecrit(llm, klon, iim, jjm+1, d_tr_th(1, 1, it), zx_tmp_3d)           CALL gr_fi_ecrit(llm, klon, iim, jjm+1, d_tr_th(1, 1, it), zx_tmp_3d)
530           CALL histwrite(nid_tra, "d_tr_th_"//tnom(it+2), itau_w, zx_tmp_3d, &           CALL histwrite(nid_tra, "d_tr_th_"//tnom(it+2), itau_w, zx_tmp_3d)
               iim*(jjm+1)*llm, ndex3d)  
531           CALL gr_fi_ecrit(llm, klon, iim, jjm+1, d_tr_cv(1, 1, it), zx_tmp_3d)           CALL gr_fi_ecrit(llm, klon, iim, jjm+1, d_tr_cv(1, 1, it), zx_tmp_3d)
532           CALL histwrite(nid_tra, "d_tr_cv_"//tnom(it+2), itau_w, zx_tmp_3d, &           CALL histwrite(nid_tra, "d_tr_cv_"//tnom(it+2), itau_w, zx_tmp_3d)
               iim*(jjm+1)*llm, ndex3d)  
533           CALL gr_fi_ecrit(llm, klon, iim, jjm+1, d_tr_cl(1, 1, it), zx_tmp_3d)           CALL gr_fi_ecrit(llm, klon, iim, jjm+1, d_tr_cl(1, 1, it), zx_tmp_3d)
534           CALL histwrite(nid_tra, "d_tr_cl_"//tnom(it+2), itau_w, zx_tmp_3d, &           CALL histwrite(nid_tra, "d_tr_cl_"//tnom(it+2), itau_w, zx_tmp_3d)
               iim*(jjm+1)*llm, ndex3d)  
535        ENDDO        ENDDO
536    
       CALL gr_fi_ecrit(1, klon, iim, jjm+1, yu1, zx_tmp_2d)  
       CALL histwrite(nid_tra, "pyu1", itau_w, zx_tmp_2d, &  
            iim*(jjm+1), ndex2d)  
   
       CALL gr_fi_ecrit(1, klon, iim, jjm+1, yv1, zx_tmp_2d)  
       CALL histwrite(nid_tra, "pyv1", itau_w, zx_tmp_2d, &  
            iim*(jjm+1), ndex2d)  
   
       CALL gr_fi_ecrit(1, klon, iim, jjm+1, pftsol1, zx_tmp_2d)  
       CALL histwrite(nid_tra, "ftsol1", itau_w, zx_tmp_2d, &  
            iim*(jjm+1), ndex2d)  
   
       CALL gr_fi_ecrit(1, klon, iim, jjm+1, pftsol2, zx_tmp_2d)  
       CALL histwrite(nid_tra, "ftsol2", itau_w, zx_tmp_2d, &  
            iim*(jjm+1), ndex2d)  
   
       CALL gr_fi_ecrit(1, klon, iim, jjm+1, pftsol3, zx_tmp_2d)  
       CALL histwrite(nid_tra, "ftsol3", itau_w, zx_tmp_2d, &  
            iim*(jjm+1), ndex2d)  
   
       CALL gr_fi_ecrit(1, klon, iim, jjm+1, pftsol4, zx_tmp_2d)  
       CALL histwrite(nid_tra, "ftsol4", itau_w, zx_tmp_2d, &  
            iim*(jjm+1), ndex2d)  
   
       CALL gr_fi_ecrit(1, klon, iim, jjm+1, ppsrf1, zx_tmp_2d)  
       CALL histwrite(nid_tra, "psrf1", itau_w, zx_tmp_2d, &  
            iim*(jjm+1), ndex2d)  
   
       CALL gr_fi_ecrit(1, klon, iim, jjm+1, ppsrf2, zx_tmp_2d)  
       CALL histwrite(nid_tra, "psrf2", itau_w, zx_tmp_2d, &  
            iim*(jjm+1), ndex2d)  
   
       CALL gr_fi_ecrit(1, klon, iim, jjm+1, ppsrf3, zx_tmp_2d)  
       CALL histwrite(nid_tra, "psrf3", itau_w, zx_tmp_2d, &  
            iim*(jjm+1), ndex2d)  
   
       CALL gr_fi_ecrit(1, klon, iim, jjm+1, ppsrf4, zx_tmp_2d)  
       CALL histwrite(nid_tra, "psrf4", itau_w, zx_tmp_2d, &  
            iim*(jjm+1), ndex2d)  
537        CALL gr_fi_ecrit(llm, klon, iim, jjm+1, pplay, zx_tmp_3d)        CALL gr_fi_ecrit(llm, klon, iim, jjm+1, pplay, zx_tmp_3d)
538        CALL histwrite(nid_tra, "pplay", itau_w, zx_tmp_3d, &        CALL histwrite(nid_tra, "pplay", itau_w, zx_tmp_3d)
            iim*(jjm+1)*llm, ndex3d)  
539    
540        CALL gr_fi_ecrit(llm, klon, iim, jjm+1, t_seri, zx_tmp_3d)        CALL gr_fi_ecrit(llm, klon, iim, jjm+1, t_seri, zx_tmp_3d)
541        CALL histwrite(nid_tra, "t", itau_w, zx_tmp_3d, &        CALL histwrite(nid_tra, "t", itau_w, zx_tmp_3d)
            iim*(jjm+1)*llm, ndex3d)  
       CALL gr_fi_ecrit(llm, klon, iim, jjm+1, pmfu, zx_tmp_3d)  
       CALL histwrite(nid_tra, "mfu", itau_w, zx_tmp_3d, &  
            iim*(jjm+1)*llm, ndex3d)  
       CALL gr_fi_ecrit(llm, klon, iim, jjm+1, pmfd, zx_tmp_3d)  
       CALL histwrite(nid_tra, "mfd", itau_w, zx_tmp_3d, &  
            iim*(jjm+1)*llm, ndex3d)  
       CALL gr_fi_ecrit(llm, klon, iim, jjm+1, pen_u, zx_tmp_3d)  
       CALL histwrite(nid_tra, "en_u", itau_w, zx_tmp_3d, &  
            iim*(jjm+1)*llm, ndex3d)  
       CALL gr_fi_ecrit(llm, klon, iim, jjm+1, pen_d, zx_tmp_3d)  
       CALL histwrite(nid_tra, "en_d", itau_w, zx_tmp_3d, &  
            iim*(jjm+1)*llm, ndex3d)  
       CALL gr_fi_ecrit(llm, klon, iim, jjm+1, pde_d, zx_tmp_3d)  
       CALL histwrite(nid_tra, "de_d", itau_w, zx_tmp_3d, &  
            iim*(jjm+1)*llm, ndex3d)  
       CALL gr_fi_ecrit(llm, klon, iim, jjm+1, pde_u, zx_tmp_3d)  
       CALL histwrite(nid_tra, "de_u", itau_w, zx_tmp_3d, &  
            iim*(jjm+1)*llm, ndex3d)  
       CALL gr_fi_ecrit(llm, klon, iim, jjm+1, coefh, zx_tmp_3d)  
       CALL histwrite(nid_tra, "coefh", itau_w, zx_tmp_3d, &  
            iim*(jjm+1)*llm, ndex3d)  
   
       ! abder  
542    
543        if (ok_sync) then        if (ok_sync) then
544           call histsync(nid_tra)           call histsync(nid_tra)
# Line 630  contains Line 560  contains
560      use advtrac_m, only: niadv, tnom, ttext      use advtrac_m, only: niadv, tnom, ttext
561      use dimphy, only: klon      use dimphy, only: klon
562      use clesphys, only: ecrit_tra      use clesphys, only: ecrit_tra
563      use grid_change, only: gr_phy_write      use grid_change, only: gr_phy_write_2d
564      use phyetat0_m, only: rlon, rlat      use phyetat0_m, only: rlon, rlat
565    
566      INTEGER, intent(out):: nid_tra      INTEGER, intent(out):: nid_tra
# Line 653  contains Line 583  contains
583      !---------------------------------------------------------      !---------------------------------------------------------
584    
585      CALL ymds2ju(annee_ref, month=1, day=day_ref, sec=0.0, julian=zjulian)      CALL ymds2ju(annee_ref, month=1, day=day_ref, sec=0.0, julian=zjulian)
586      zx_lat(:, :) = gr_phy_write(rlat)      zx_lat(:, :) = gr_phy_write_2d(rlat)
587      CALL histbeg_totreg("histrac", iim, rlon(2:iim+1), jjm+1, zx_lat(1, :), &      CALL histbeg_totreg("histrac", rlon(2:iim+1), zx_lat(1, :), &
588           1, iim, 1, jjm+1, itau_phy, zjulian, pdtphys, nhori, nid_tra)           1, iim, 1, jjm+1, itau_phy, zjulian, pdtphys, nhori, nid_tra)
589      CALL histvert(nid_tra, "presnivs", "Vertical levels", "mb", llm, &      CALL histvert(nid_tra, "presnivs", "Vertical levels", "mb", llm, &
590           presnivs, nvert)           presnivs, nvert)
# Line 663  contains Line 593  contains
593      zout = pdtphys * REAL(ecrit_tra)      zout = pdtphys * REAL(ecrit_tra)
594    
595      CALL histdef(nid_tra, "phis", "Surface geop. height", "-", &      CALL histdef(nid_tra, "phis", "Surface geop. height", "-", &
596           iim, jjm+1, nhori, 1, 1, 1, -99, 32, &           iim, jjm+1, nhori, 1, 1, 1, -99, &
597           "once",  zsto, zout)           "once",  zsto, zout)
598      CALL histdef(nid_tra, "aire", "Grid area", "-", &      CALL histdef(nid_tra, "aire", "Grid area", "-", &
599           iim, jjm+1, nhori, 1, 1, 1, -99, 32, &           iim, jjm+1, nhori, 1, 1, 1, -99, &
600           "once",  zsto, zout)           "once",  zsto, zout)
601      CALL histdef(nid_tra, "zmasse", "column density of air in cell", &      CALL histdef(nid_tra, "zmasse", "column density of air in cell", &
602           "kg m-2", iim, jjm + 1, nhori, llm, 1, llm, nvert, 32, "ave(X)", &           "kg m-2", iim, jjm + 1, nhori, llm, 1, llm, nvert, "ave(X)", &
603           zsto, zout)           zsto, zout)
604    
605      DO it=1, nqmax      DO it = 1, nqmax
606         ! champ 2D         ! champ 2D
607         iq=it+2         iq=it+2
608         iiq=niadv(iq)         iiq=niadv(iq)
609         CALL histdef(nid_tra, tnom(iq), ttext(iiq), "U/kga", iim, jjm+1, &         CALL histdef(nid_tra, tnom(iq), ttext(iiq), "U/kga", iim, jjm+1, &
610              nhori, llm, 1, llm, nvert, 32, "ave(X)", zsto, zout)              nhori, llm, 1, llm, nvert, "ave(X)", zsto, zout)
611         if (lessivage) THEN         if (lessivage) THEN
612            CALL histdef(nid_tra, "fl"//tnom(iq), "Flux "//ttext(iiq), &            CALL histdef(nid_tra, "fl"//tnom(iq), "Flux "//ttext(iiq), &
613                 "U/m2/s", iim, jjm+1, nhori, llm, 1, llm, nvert, 32, &                 "U/m2/s", iim, jjm+1, nhori, llm, 1, llm, nvert, &
614                 "ave(X)", zsto, zout)                 "ave(X)", zsto, zout)
615         endif         endif
616    
617         !---Ajout Olivia         !---Ajout Olivia
618         CALL histdef(nid_tra, "d_tr_th_"//tnom(iq), &         CALL histdef(nid_tra, "d_tr_th_"//tnom(iq), &
619              "tendance thermique"// ttext(iiq), "?", &              "tendance thermique"// ttext(iiq), "?", &
620              iim, jjm+1, nhori, llm, 1, llm, nvert, 32, &              iim, jjm+1, nhori, llm, 1, llm, nvert, &
621              "ave(X)", zsto, zout)              "ave(X)", zsto, zout)
622         CALL histdef(nid_tra, "d_tr_cv_"//tnom(iq), &         CALL histdef(nid_tra, "d_tr_cv_"//tnom(iq), &
623              "tendance convection"// ttext(iiq), "?", &              "tendance convection"// ttext(iiq), "?", &
624              iim, jjm+1, nhori, llm, 1, llm, nvert, 32, &              iim, jjm+1, nhori, llm, 1, llm, nvert, &
625              "ave(X)", zsto, zout)              "ave(X)", zsto, zout)
626         CALL histdef(nid_tra, "d_tr_cl_"//tnom(iq), &         CALL histdef(nid_tra, "d_tr_cl_"//tnom(iq), &
627              "tendance couche limite"// ttext(iiq), "?", &              "tendance couche limite"// ttext(iiq), "?", &
628              iim, jjm+1, nhori, llm, 1, llm, nvert, 32, &              iim, jjm+1, nhori, llm, 1, llm, nvert, &
629              "ave(X)", zsto, zout)              "ave(X)", zsto, zout)
630         !---fin Olivia             !---fin Olivia    
631    
632      ENDDO      ENDDO
633    
634      CALL histdef(nid_tra, "pyu1", "Vent niv 1", "-", &      CALL histdef(nid_tra, "pplay", "", "-", &
635           iim, jjm+1, nhori, 1, 1, 1, -99, 32, &           iim, jjm+1, nhori, llm, 1, llm, nvert, &
          "inst(X)", zout, zout)  
   
     CALL histdef(nid_tra, "pyv1", "Vent niv 1", "-", &  
          iim, jjm+1, nhori, 1, 1, 1, -99, 32, &  
          "inst(X)",  zout, zout)  
     CALL histdef(nid_tra, "psrf1", "nature sol", "-", &  
          iim, jjm+1, nhori, 1, 1, 1, -99, 32, &  
          "inst(X)",  zout, zout)  
     CALL histdef(nid_tra, "psrf2", "nature sol", "-", &  
          iim, jjm+1, nhori, 1, 1, 1, -99, 32, &  
          "inst(X)",  zout, zout)  
     CALL histdef(nid_tra, "psrf3", "nature sol", "-", &  
          iim, jjm+1, nhori, 1, 1, 1, -99, 32, &  
          "inst(X)",  zout, zout)  
     CALL histdef(nid_tra, "psrf4", "nature sol", "-", &  
          iim, jjm+1, nhori, 1, 1, 1, -99, 32, &  
          "inst(X)",  zout, zout)  
     CALL histdef(nid_tra, "ftsol1", "temper sol", "-", &  
          iim, jjm+1, nhori, 1, 1, 1, -99, 32, &  
          "inst(X)",  zout, zout)  
     CALL histdef(nid_tra, "ftsol2", "temper sol", "-", &  
          iim, jjm+1, nhori, 1, 1, 1, -99, 32, &  
          "inst(X)",  zout, zout)  
     CALL histdef(nid_tra, "ftsol3", "temper sol", "-", &  
          iim, jjm+1, nhori, 1, 1, 1, -99, 32, &  
          "inst(X)",  zout, zout)  
     CALL histdef(nid_tra, "ftsol4", "temper sol", "-", &  
          iim, jjm+1, nhori, 1, 1, 1, -99, 32, &  
          "inst(X)",  zout, zout)  
     CALL histdef(nid_tra, "pplay", "flux u mont", "-", &  
          iim, jjm+1, nhori, llm, 1, llm, nvert, 32, &  
636           "inst(X)", zout, zout)           "inst(X)", zout, zout)
637      CALL histdef(nid_tra, "t", "flux u mont", "-", &      CALL histdef(nid_tra, "t", "", "-", &
638           iim, jjm+1, nhori, llm, 1, llm, nvert, 32, &           iim, jjm+1, nhori, llm, 1, llm, nvert, &
639           "inst(X)", zout, zout)           "inst(X)", zout, zout)
     CALL histdef(nid_tra, "mfu", "flux u mont", "-", &  
          iim, jjm+1, nhori, llm, 1, llm, nvert, 32, &  
          "ave(X)", zsto, zout)  
     CALL histdef(nid_tra, "mfd", "flux u decen", "-", &  
          iim, jjm+1, nhori, llm, 1, llm, nvert, 32, &  
          "ave(X)", zsto, zout)  
     CALL histdef(nid_tra, "en_u", "flux u mont", "-", &  
          iim, jjm+1, nhori, llm, 1, llm, nvert, 32, &  
          "ave(X)", zsto, zout)  
     CALL histdef(nid_tra, "en_d", "flux u mont", "-", &  
          iim, jjm+1, nhori, llm, 1, llm, nvert, 32, &  
          "ave(X)", zsto, zout)  
     CALL histdef(nid_tra, "de_d", "flux u mont", "-", &  
          iim, jjm+1, nhori, llm, 1, llm, nvert, 32, &  
          "ave(X)", zsto, zout)  
     CALL histdef(nid_tra, "de_u", "flux u decen", "-", &  
          iim, jjm+1, nhori, llm, 1, llm, nvert, 32, &  
          "ave(X)", zsto, zout)  
     CALL histdef(nid_tra, "coefh", "turbulent coef", "-", &  
          iim, jjm+1, nhori, llm, 1, llm, nvert, 32, &  
          "ave(X)", zsto, zout)  
640    
641      CALL histend(nid_tra)      CALL histend(nid_tra)
642    
Legend:
Removed from v.14  
changed lines
  Added in v.15
  ViewVC Help
Powered by ViewVC 1.1.21