Changeset 902 for trunk/LMDZ.VENUS/libf


Ignore:
Timestamp:
Mar 13, 2013, 11:06:21 AM (12 years ago)
Author:
slebonnois
Message:

SL: small modifs on Venus outputs

Location:
trunk/LMDZ.VENUS/libf/phyvenus
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.VENUS/libf/phyvenus/ini_histday.h

    r892 r902  
    66         zsto = dtime
    77         zout = dtime * REAL(ecrit_day)
    8          zsto1= dtime * REAL(ecrit_day)
    98
    109         idayref = day_ref
     
    1211
    1312         call histbeg_phy("histday.nc",itau_phy,
    14      .                    zjulian,zsto,nhori,nid_day)
     13     .                    zjulian,dtime,nhori,nid_day)
    1514
    1615!$OMP MASTER
     
    3938     .                "ave(X)", zsto,zout)
    4039c
    41          CALL histdef(nid_day, "ue", "Zonal energy transport", "-",
    42      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    43      .                "ave(X)", zsto,zout)
    44 c
    45          CALL histdef(nid_day, "ve", "Merid energy transport", "-",
    46      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    47      .                "ave(X)", zsto,zout)
     40c        CALL histdef(nid_day, "ue", "Zonal energy transport", "-",
     41c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     42c    .                "ave(X)", zsto,zout)
     43c
     44c        CALL histdef(nid_day, "ve", "Merid energy transport", "-",
     45c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     46c    .                "ave(X)", zsto,zout)
     47c
     48c        CALL histdef(nid_day, "cdragh", "Drag coef on T", "-",
     49c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     50c    .                "ave(X)", zsto,zout)
     51c
     52c        CALL histdef(nid_day, "cdragm", "Drag coef on U", "-",
     53c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     54c    .                "ave(X)", zsto,zout)
    4855c
    4956      ENDIF !lev_histday.GE.1
     
    7885     .                "ave(X)", zsto,zout)
    7986c
     87         CALL histdef(nid_day, "dudyn", "Dynamics dU", "m/s2",
     88     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     89     .                "ave(X)", zsto,zout)
     90c
     91         CALL histdef(nid_day, "duvdf", "Boundary-layer dU", "m/s2",
     92     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     93     .                "ave(X)", zsto,zout)
     94c
     95c        CALL histdef(nid_day, "mang", "Angular momentum", "kg m2/s",
     96c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     97c    .                "ave(X)", zsto,zout)
     98c
     99c        CALL histdef(nid_day, "Kz", "vertical diffusion coef", "m2/s",
     100c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     101c    .                "ave(X)", zsto,zout)
     102c
     103c plusieurs traceurs
     104          if (iflag_trac.eq.1) THEN
     105            DO iq=1,nqmax
     106             IF (iq.LE.99) THEN
     107          WRITE(str2,'(i2.2)') iq
     108          CALL histdef(nid_day, tname(iq), ttext(iq), "ppm",
     109     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     110     .                "ave(X)", zsto,zout)
     111             ELSE
     112          PRINT*, "Trop de traceurs"
     113          CALL abort
     114             ENDIF
     115            ENDDO
     116          endif
     117c
     118         CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
     119     .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     120     .                "ave(X)", zsto,zout)
     121c
    80122      ENDIF !lev_histday.GE.2
    81123c
     
    87129c 2D
    88130c
    89          CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
    90      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    91      .                "ave(X)", zsto1,zout)
    92 c
    93131         CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
    94132     .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    95      .                "ave(X)", zsto1,zout)
     133     .                "ave(X)", zsto,zout)
    96134c
    97135         CALL histdef(nid_day, "sols", "Solar rad. at surf.", "W/m2",
    98136     .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    99      .                "ave(X)", zsto1,zout)
     137     .                "ave(X)", zsto,zout)
    100138c
    101139         CALL histdef(nid_day, "soll", "IR rad. at surface", "W/m2",
    102140     .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    103      .                "ave(X)", zsto1,zout)
     141     .                "ave(X)", zsto,zout)
    104142c
    105143c 3D
     
    107145         CALL histdef(nid_day, "SWnet", "Net SW flux","W/m2",
    108146     .                iim,jj_nb,nhori, klev,1,klev,nvert,
    109      .                32, "ave(X)", zsto1,zout)
     147     .                32, "ave(X)", zsto,zout)
    110148c
    111149         CALL histdef(nid_day, "LWnet", "Net LW flux","W/m2",
    112150     .                iim,jj_nb,nhori, klev,1,klev,nvert,
    113      .                32, "ave(X)", zsto1,zout)
     151     .                32, "ave(X)", zsto,zout)
     152c
     153         CALL histdef(nid_day, "fluxvdf", "PBL net flux","W/m2",
     154     .                iim,jj_nb,nhori, klev,1,klev,nvert,
     155     .                32, "ave(X)", zsto,zout)
     156c
     157         CALL histdef(nid_day, "fluxdyn", "Dyn. net flux","W/m2",
     158     .                iim,jj_nb,nhori, klev,1,klev,nvert,
     159     .                32, "ave(X)", zsto,zout)
     160c
     161         CALL histdef(nid_day, "fluxajs", "Dry adj. net flux","W/m2",
     162     .                iim,jj_nb,nhori, klev,1,klev,nvert,
     163     .                32, "ave(X)", zsto,zout)
     164c
     165c        CALL histdef(nid_day, "fluxec", "Cin. net flux","W/m2",
     166c    .                iim,jj_nb,nhori, klev,1,klev,nvert,
     167c    .                32, "ave(X)", zsto,zout)
    114168c
    115169      ENDIF !lev_histday.GE.3
     
    122176     .                "ave(X)", zsto,zout)
    123177c
    124 c         CALL histdef(nid_day, "dtphy", "Physics dT", "K/s",
    125 c     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    126 c     .                "ave(X)", zsto,zout)
     178c        CALL histdef(nid_day, "dtphy", "Physics dT", "K/s",
     179c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     180c    .                "ave(X)", zsto,zout)
    127181c
    128182         CALL histdef(nid_day, "dtvdf", "Boundary-layer dT", "K/s",
     
    142196     .                "ave(X)", zsto,zout)
    143197c
    144 c         CALL histdef(nid_day, "dtec", "Cinetic dissip dT", "K/s",
    145 c     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    146 c     .                "ave(X)", zsto,zout)
    147 c
    148          CALL histdef(nid_mth, "dudyn", "Dynamics dU", "m/s2",
    149      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    150      .                "ave(X)", zsto,zout)
    151 c
    152          CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
    153      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    154      .                "ave(X)", zsto,zout)
     198c        CALL histdef(nid_day, "dtec", "Cinetic dissip dT", "K/s",
     199c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     200c    .                "ave(X)", zsto,zout)
    155201c
    156202         CALL histdef(nid_mth, "duajs", "Dry convection dU", "m/s2",
     
    174220c-------------------------------------------------------
    175221      IF(lev_histday.GE.5) THEN
    176 c
    177          if (nqmax.GE.3) THEN
    178            DO iq=1,nqmax-2
    179             IF (iq.LE.99) THEN
    180          WRITE(str2,'(i2.2)') iq
    181          CALL histdef(nid_day, "trac"//str2, "Tracer No."//str2, "-",
    182      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    183      .                "ave(X)", zsto,zout)
    184             ELSE
    185          PRINT*, "Trop de traceurs"
    186          CALL abort
    187             ENDIF
    188            ENDDO
    189          endif
    190222c
    191223c        call histdef(nid_day, "taux",
  • trunk/LMDZ.VENUS/libf/phyvenus/ini_histins.h

    r892 r902  
    44      IF (ok_instan) THEN
    55
    6           zsto1= dtime
    7           zsto = dtime * ecrit_ins
    8           zout = dtime * ecrit_ins
     6         zsto = dtime * ecrit_ins
     7         zout = dtime * ecrit_ins
    98
    109         idayref = day_ref
     
    1211
    1312         call histbeg_phy("histins.nc",itau_phy,
    14      .                    zjulian,zsto,nhori,nid_ins)
     13     .                    zjulian,dtime,nhori,nid_ins)
    1514
    1615!$OMP MASTER
     
    4746c    .                "ins(X)", zsto,zout)
    4847c
    49          CALL histdef(nid_ins, "cdragh", "Drag coef on T", "-",
    50      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    51      .                "ins(X)", zsto,zout)
    52 c
    53          CALL histdef(nid_ins, "cdragm", "Drag coef on U", "-",
    54      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    55      .                "ins(X)", zsto,zout)
     48c        CALL histdef(nid_ins, "cdragh", "Drag coef on T", "-",
     49c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     50c    .                "ins(X)", zsto,zout)
     51c
     52c        CALL histdef(nid_ins, "cdragm", "Drag coef on U", "-",
     53c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     54c    .                "ins(X)", zsto,zout)
    5655c
    5756      ENDIF !lev_histday.GE.1
     
    9897c    .                "ins(X)", zsto,zout)
    9998c
    100          CALL histdef(nid_ins, "Kz", "vertical diffusion coef", "m2/s",
    101      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    102      .                "ins(X)", zsto,zout)
     99c        CALL histdef(nid_ins, "Kz", "vertical diffusion coef", "m2/s",
     100c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     101c    .                "ins(X)", zsto,zout)
    103102c
    104103c plusieurs traceurs
  • trunk/LMDZ.VENUS/libf/phyvenus/ini_histmth.h

    r892 r902  
    77         zout = dtime * ecrit_mth
    88
    9 c zsto1: pour des valeurs "instantanees" mensuelles
    10          zsto1 = dtime * ecrit_mth
    11 c zsto2: pour des flux radiatifs calcules avec un autre pas de temps
    12          zsto2 = dtime * radpas
    13          PRINT*,' zsto,zsto1,zsto2,zout=',zsto, zsto1, zsto2,zout
    14 
    159         idayref = day_ref
    1610         CALL ymds2ju(annee_ref, 1, idayref, zero, zjulian)
    1711
    1812         call histbeg_phy("histmth.nc",itau_phy,
    19      .                    zjulian,zsto,nhori,nid_mth)
     13     .                    zjulian,dtime,nhori,nid_mth)
    2014
    2115!$OMP MASTER
     
    5246c    .                "ave(X)", zsto,zout)
    5347c
    54          CALL histdef(nid_mth, "cdragh", "Drag coef on T", "-",
    55      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    56      .                "ave(X)", zsto,zout)
    57 c
    58          CALL histdef(nid_mth, "cdragm", "Drag coef on U", "-",
    59      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    60      .                "ave(X)", zsto,zout)
     48c        CALL histdef(nid_mth, "cdragh", "Drag coef on T", "-",
     49c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     50c    .                "ave(X)", zsto,zout)
     51c
     52c        CALL histdef(nid_mth, "cdragm", "Drag coef on U", "-",
     53c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     54c    .                "ave(X)", zsto,zout)
    6155c
    6256      ENDIF !lev_histmth.GE.1
     
    10397c    .                "ave(X)", zsto,zout)
    10498c
    105          CALL histdef(nid_mth, "Kz", "vertical diffusion coef", "m2/s",
    106      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    107      .                "ave(X)", zsto,zout)
     99c        CALL histdef(nid_mth, "Kz", "vertical diffusion coef", "m2/s",
     100c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     101c    .                "ave(X)", zsto,zout)
    108102c
    109103c plusieurs traceurs
     
    124118         CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
    125119     .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    126      .                "ave(X)", zsto2,zout)
     120     .                "ave(X)", zsto,zout)
    127121c
    128122      ENDIF !lev_histmth.GE.2
     
    137131         CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
    138132     .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    139      .                "ave(X)", zsto2,zout)
     133     .                "ave(X)", zsto,zout)
    140134c
    141135         CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
    142136     .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    143      .                "ave(X)", zsto2,zout)
     137     .                "ave(X)", zsto,zout)
    144138c
    145139         CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
    146140     .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    147      .                "ave(X)", zsto2,zout)
     141     .                "ave(X)", zsto,zout)
    148142c
    149143c 3D
     
    151145         CALL histdef(nid_mth, "SWnet", "Net SW flux","W/m2",
    152146     .                iim,jj_nb,nhori, klev,1,klev,nvert,
    153      .                32, "ave(X)", zsto2,zout)
     147     .                32, "ave(X)", zsto,zout)
    154148c
    155149         CALL histdef(nid_mth, "LWnet", "Net LW flux","W/m2",
    156150     .                iim,jj_nb,nhori, klev,1,klev,nvert,
    157      .                32, "ave(X)", zsto2,zout)
     151     .                32, "ave(X)", zsto,zout)
    158152c
    159153         CALL histdef(nid_mth, "fluxvdf", "PBL net flux","W/m2",
  • trunk/LMDZ.VENUS/libf/phyvenus/write_histday.h

    r892 r902  
    44      IF (ok_journe) THEN
    55
    6          zsto = dtime
    7          zout = dtime * REAL(ecrit_day)
    86         itau_w = itau_phy + itap
    97
     
    5452
    5553      ENDIF !lev_histday.GE.2
    56 c
     54
    5755c-------------------------------------------------------
    5856      IF(lev_histday.GE.3) THEN
    59 c
     57
    6058cccccccccccccccccc  Radiative transfer
    6159
     
    7674
    7775      ENDIF !lev_histday.GE.3
    78 c
     76
    7977c-------------------------------------------------------
    8078      IF(lev_histday.GE.4) THEN
     
    105103
    106104      ENDIF !lev_histday.GE.4
    107 c
     105
    108106c-------------------------------------------------------
    109107      IF(lev_histday.GE.5) THEN
     
    116114      ENDIF !lev_histday.GE.5
    117115c-------------------------------------------------------
    118 c
     116
    119117      if (ok_sync) then
    120118        call histsync(nid_day)
  • trunk/LMDZ.VENUS/libf/phyvenus/write_histins.h

    r892 r902  
    44      IF (ok_instan) THEN
    55
    6          zsto = dtime * ecrit_ins
    7          zout = dtime * ecrit_ins
    86         itau_w = itau_phy + itap
    97
     
    107105
    108106c-------------------------------------------------------
    109        IF(lev_histday.GE.5) THEN
     107      IF(lev_histday.GE.5) THEN
    110108
    111109c     call histwrite_phy(nid_ins,.false.,"taux_",itau_w,fluxu)
  • trunk/LMDZ.VENUS/libf/phyvenus/write_histmth.h

    r892 r902  
    44      IF (ok_mensuel) THEN
    55
    6          zsto = dtime
    7          zout = dtime * ecrit_mth
    86         itau_w = itau_phy + itap
    97
     
    105103
    106104      ENDIF !lev_histmth.GE.4
    107 c
     105
    108106c-------------------------------------------------------
    109        IF(lev_histmth.GE.5) THEN
     107      IF(lev_histmth.GE.5) THEN
    110108
    111109c     call histwrite_phy(nid_mth,.false.,"taux_",itau_w,fluxu)
     
    116114      ENDIF !lev_histmth.GE.5
    117115c-------------------------------------------------------
    118 c
     116
    119117      if (ok_sync) then
    120118        call histsync(nid_mth)
Note: See TracChangeset for help on using the changeset viewer.