Ignore:
Timestamp:
Oct 7, 2013, 6:42:03 PM (11 years ago)
Author:
slebonnois
Message:

SL: Titan runs ! see DOC/chantiers/commit_importants.log

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r902 r1056  
    11      IF (ok_journe) THEN
    2 c
     2
    33         zsto = dtime
    4          zout = dtime * FLOAT(ecrit_day)
     4         zout = dtime * REAL(ecrit_day)
    55c zsto1: pour des flux radiatifs calcules tous les radpas appels physiq
    6          zsto1= dtime * FLOAT(radpas)
    7 c
     6         zsto1= dtime * REAL(radpas)
     7
    88         idayref = day_ref
    99         CALL ymds2ju(annee_ref, 1, idayref, zero, zjulian)
    10 c
    11          CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlond,zx_lon)
    12          DO i = 1, iim
    13             zx_lon(i,1) = rlond(i+jjmp1-jjm)
    14             zx_lon(i,jjmp1) = rlond(i+jjmp1-jjm)
    15          ENDDO
    16          CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlatd,zx_lat)
    17          CALL histbeg("histday", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
    18      .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
     10
     11         CALL histbeg_phy("histday.nc", itau_phy, zjulian, dtime,
    1912     .                 nhori, nid_day)
    20          write(*,*)'Journee ', itau_phy, zjulian
    21          CALL histvert(nid_day, "presnivs", "Vertical levels", "mb",
     13
     14!$OMP MASTER
     15         CALL histvert(nid_day, "presnivs", "Vertical levels", "Pa",
    2216     .                 klev, presnivs, nvert)
    2317
    2418c-------------------------------------------------------
    2519      IF(lev_histday.GE.1) THEN
    26 c
     20
    2721ccccccccccccc 2D fields, invariables
    28 c
     22
    2923         CALL histdef(nid_day, "phis", "Surface geop. height", "-",
    30      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     24     .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    3125     .                "once",  zsto,zout)
    32 c
     26
    3327         CALL histdef(nid_day, "aire", "Grid area", "-",
    34      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     28     .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    3529     .                "once",  zsto,zout)
    36 c
     30
    3731ccccccc axe Ls
    3832         CALL histdef(nid_day, "ls", "Solar longitude", "degrees",
    39      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    40      .                "ave(X)", zsto,zout)
    41 c
     33     .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     34     .                "ave(X)", zsto,zout)
     35
    4236ccccccccccccc 2D fields, variables
    43 c
     37
    4438         CALL histdef(nid_day, "tsol", "Surface Temperature", "K",
    45      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    46      .                "ave(X)", zsto,zout)
    47 c
     39     .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     40     .                "ave(X)", zsto,zout)
     41
    4842         CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",
    49      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    50      .                "ave(X)", zsto,zout)
    51 c
     43     .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     44     .                "ave(X)", zsto,zout)
     45
    5246c        CALL histdef(nid_day, "ue", "Zonal energy transport", "-",
    53 c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     47c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    5448c    .                "ave(X)", zsto,zout)
    55 c
     49
    5650c        CALL histdef(nid_day, "ve", "Merid energy transport", "-",
    57 c     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     51c     .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    5852c     .                "ave(X)", zsto,zout)
    59 c
     53
    6054      ENDIF !lev_histday.GE.1
    61 c
     55
    6256c-------------------------------------------------------
    6357      IF(lev_histday.GE.2) THEN
    64 c
     58
    6559ccccccccccccc 3D fields, basics
    66 c
     60
    6761         CALL histdef(nid_day, "temp", "Air temperature", "K",
    68      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    69      .                "ave(X)", zsto,zout)
    70 c
     62     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     63     .                "ave(X)", zsto,zout)
     64
    7165         CALL histdef(nid_day, "pres", "Air pressure", "Pa",
    72      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    73      .                "ave(X)", zsto,zout)
    74 c
     66     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     67     .                "ave(X)", zsto,zout)
     68
    7569         CALL histdef(nid_day, "geop", "Geopotential height", "m",
    76      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    77      .                "ave(X)", zsto,zout)
    78 c
     70     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     71     .                "ave(X)", zsto,zout)
     72
    7973         CALL histdef(nid_day, "vitu", "Zonal wind", "m/s",
    80      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    81      .                "ave(X)", zsto,zout)
    82 c
     74     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     75     .                "ave(X)", zsto,zout)
     76
    8377         CALL histdef(nid_day, "vitv", "Meridional wind", "m/s",
    84      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    85      .                "ave(X)", zsto,zout)
    86 c
     78     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     79     .                "ave(X)", zsto,zout)
     80
    8781         CALL histdef(nid_day, "vitw", "Vertical wind", "Pa/s",
    88      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    89      .                "ave(X)", zsto,zout)
    90 c
     82     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     83     .                "ave(X)", zsto,zout)
     84
    9185         CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
    92      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     86     .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    9387     .                "ave(X)", zsto1,zout)
    94 c
     88
     89         CALL histdef(nid_day, "duvdf", "Boundary-layer dU", "m/s2",
     90     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     91     .                "ave(X)", zsto,zout)
     92
    9593         CALL histdef(nid_day, "dudyn", "Dynamics dU", "m/s2",
    96      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    97      .                "ave(X)", zsto,zout)
    98 c
    99          CALL histdef(nid_day, "duvdf", "Boundary-layer dU", "m/s2",
    100      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    101      .                "ave(X)", zsto,zout)
    102 c
     94     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     95     .                "ave(X)", zsto,zout)
     96
    10397cccccccccccccccccc  Tracers
    104 c
     98
    10599         if (iflag_trac.eq.1) THEN
    106100          if (microfi.ge.1) then
    107101c           DO iq=1,nmicro
    108102c             CALL histdef(nid_day, tname(iq), ttext(iq), "n/m2",
    109 c     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     103c     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    110104c     .                "ave(X)", zsto,zout)
    111105c           ENDDO
    112106             CALL histdef(nid_day, "qaer","nb tot aer" , "n/m2",
    113      .                    iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     107     .                    iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    114108     .                    "ave(X)", zsto,zout)
    115109
    116110            if (clouds.eq.1) then
    117111             CALL histdef(nid_day, "qnoy","nb tot noy" , "n/m2",
    118      .                    iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     112     .                    iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    119113     .                    "ave(X)", zsto,zout)
    120114             CALL histdef(nid_day, "qgl1","V tot gl1" , "m3/m2",
    121      .                    iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     115     .                    iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    122116     .                    "ave(X)", zsto,zout)
    123117             CALL histdef(nid_day, "qgl2","V tot gl2" , "m3/m2",
    124      .                    iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     118     .                    iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    125119     .                    "ave(X)", zsto,zout)
    126120             CALL histdef(nid_day, "qgl3","V tot gl3" , "m3/m2",
    127      .                    iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     121     .                    iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    128122     .                    "ave(X)", zsto,zout)
    129123c--------------
    130124c ----- SATURATION ESP NUAGES
    131125               CALL histdef(nid_day,"ch4sat", "saturation CH4", "--",
    132      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     126     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    133127     .                "ave(X)", zsto,zout)
    134128               CALL histdef(nid_day,"c2h6sat", "saturation C2H6", "--",
    135      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     129     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    136130     .                "ave(X)", zsto,zout)
    137131               CALL histdef(nid_day,"c2h2sat", "saturation C2H2", "--",
    138      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     132     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    139133     .                "ave(X)", zsto,zout)
    140134c --------------
    141135c ----- RESERVOIR DE SURFACE
    142136               CALL histdef(nid_day, "reserv", "Reservoir surface","m",
    143      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     137     .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    144138     .                "ave(X)", zsto,zout)
    145139c --------------
    146140c ----- ECHANGE GAZ SURF/ATM (evaporation)
    147141               CALL histdef(nid_day, "evapch4", "Evaporation CH4","m",
    148      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     142     .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    149143     .                "ave(X)", zsto,zout)
    150144c --------------
    151145c ----- PRECIPITATIONS (precipitations cumulatives)
    152146               CALL histdef(nid_day,"prech4","Precip CH4","m",
    153      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     147     .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    154148     .                "ave(X)", zsto,zout)
    155149               CALL histdef(nid_day,"prec2h6","Precip C2H6",
    156      .                "m",iim,jjmp1,nhori, 1,1,1, -99, 32,
     150     .                "m",iim,jj_nb,nhori, 1,1,1, nvert, 32,
    157151     .                "ave(X)", zsto,zout)
    158152               CALL histdef(nid_day,"prec2h2","Precip C2H2",
    159      .                "m",iim,jjmp1,nhori, 1,1,1, -99, 32,
     153     .                "m",iim,jj_nb,nhori, 1,1,1, nvert, 32,
     154     .                "ave(X)", zsto,zout)
     155               CALL histdef(nid_day,"prenoy","Precip NOY",
     156     .                "um/s",iim,jj_nb,nhori, 1,1,1, nvert, 32,
     157     .                "ave(X)", zsto,zout)
     158               CALL histdef(nid_day,"preaer","Precip AER",
     159     .                "um/s",iim,jj_nb,nhori, 1,1,1, nvert, 32,
    160160     .                "ave(X)", zsto,zout)
    161161c --------------
    162162c ----- FLUX GLACE
    163163               CALL histdef(nid_day,"flxgl1", "flux gl CH4",
    164      .              "kg/m2/s",iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     164     .              "kg/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    165165     .              "ave(X)", zsto,zout)
    166166               CALL histdef(nid_day,"flxgl2", "flux gl C2H6",
    167      .              "kg/m2/s",iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     167     .              "kg/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    168168     .              "ave(X)", zsto,zout)
    169169               CALL histdef(nid_day,"flxgl3", "flux gl C2H2",
    170      .              "kg/m2/s",iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     170     .              "kg/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    171171     .              "ave(X)", zsto,zout)
    172172c --------------
    173173c ----- RAYON DES GOUTTES
    174174               CALL histdef(nid_day,"rcldbar", "rayon moyen goutte",
    175      .                "m",iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     175     .                "m",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    176176     .                "ave(X)", zsto,zout)
    177177            endif
    … …  
    182182           DO iq=nmicro+1,nqmax
    183183         CALL histdef(nid_day, tname(iq), ttext(iq), "ppm",
    184      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     184     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    185185     .                "ave(X)", zsto,zout)
    186186           ENDDO
    187187          endif
    188188         endif
    189 c
     189
    190190      ENDIF !lev_histday.GE.2
    191 c
     191
    192192c-------------------------------------------------------
    193193      IF(lev_histday.GE.3) THEN
    194 c
     194
    195195cccccccccccccccccc  Radiative transfer
    196 c
     196
    197197c 2D
    198 c
     198
    199199         CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
    200      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     200     .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    201201     .                "ave(X)", zsto1,zout)
    202 c
     202
    203203         CALL histdef(nid_day, "sols", "Solar rad. at surf.", "W/m2",
    204      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     204     .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    205205     .                "ave(X)", zsto1,zout)
    206 c
     206
    207207         CALL histdef(nid_day, "soll", "IR rad. at surface", "W/m2",
    208      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     208     .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    209209     .                "ave(X)", zsto1,zout)
    210 c
     210
    211211c 3D
    212 c
     212
    213213         CALL histdef(nid_day, "SWnet", "Net SW flux","W/m2",
    214      .                iim,jjmp1,nhori, klev,1,klev,nvert,
     214     .                iim,jj_nb,nhori, klev,1,klev,nvert,
    215215     .                32, "ave(X)", zsto1,zout)
    216 c
     216
    217217         CALL histdef(nid_day, "LWnet", "Net LW flux","W/m2",
    218      .                iim,jjmp1,nhori, klev,1,klev,nvert,
     218     .                iim,jj_nb,nhori, klev,1,klev,nvert,
    219219     .                32, "ave(X)", zsto1,zout)
    220 c
     220
    221221c --------------
    222222c ----- OPACITE BRUME
    223223         DO k=7,NSPECV,10
    224            write(str1,'(i2.2)') k
    225          CALL histdef(nid_day,"thv"//str1,"Haze Opa Vis",
    226      .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
    227      .                "ave(X)",zsto1,zout)
    228          ENDDO
    229 c
     224           write(str2,'(i2.2)') k
     225         CALL histdef(nid_day,"thv"//str2,"Haze Opa Vis",
     226     .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     227     .                "ave(X)",zsto1,zout)
     228         ENDDO
     229
    230230         DO k=8,NSPECI,10
    231            write(str1,'(i2.2)') k
    232          CALL histdef(nid_day,"thi"//str1,"Haze Opa IR",
    233      .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
    234      .                "ave(X)",zsto1,zout)
    235          ENDDO
    236 c
     231           write(str2,'(i2.2)') k
     232         CALL histdef(nid_day,"thi"//str2,"Haze Opa IR",
     233     .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     234     .                "ave(X)",zsto1,zout)
     235         ENDDO
     236
    237237c --------------
    238238c ----- EXTINCTION BRUME
    239239         DO k=7,NSPECV,10
    240            write(str1,'(i2.2)') k
    241          CALL histdef(nid_day,"khv"//str1,"Haze ext Vis ",
    242      .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
    243      .                "ave(X)",zsto1,zout)
    244          ENDDO
    245 c
     240           write(str2,'(i2.2)') k
     241         CALL histdef(nid_day,"khv"//str2,"Haze ext Vis ",
     242     .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     243     .                "ave(X)",zsto1,zout)
     244         ENDDO
     245
    246246         DO k=8,NSPECI,10
    247            write(str1,'(i2.2)') k
    248          CALL histdef(nid_day,"khi"//str1,"Haze ext IR ",
    249      .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
    250      .                "ave(X)",zsto1,zout)
    251          ENDDO
    252 c
     247           write(str2,'(i2.2)') k
     248         CALL histdef(nid_day,"khi"//str2,"Haze ext IR ",
     249     .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     250     .                "ave(X)",zsto1,zout)
     251         ENDDO
     252
    253253c --------------
    254254c ----- OPACITE GAZ
    255255         DO k=7,NSPECV,10
    256            write(str1,'(i2.2)') k
    257          CALL histdef(nid_day,"tgv"//str1,"Gas Opa Vis",
    258      .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
    259      .                "ave(X)",zsto1,zout)
    260          ENDDO
    261 c
     256           write(str2,'(i2.2)') k
     257         CALL histdef(nid_day,"tgv"//str2,"Gas Opa Vis",
     258     .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     259     .                "ave(X)",zsto1,zout)
     260         ENDDO
     261
    262262         DO k=8,NSPECI,10
    263            write(str1,'(i2.2)') k
    264          CALL histdef(nid_day,"tgi"//str1,"Gas Opa IR",
    265      .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
    266      .                "ave(X)",zsto1,zout)
    267          ENDDO
    268 c
     263           write(str2,'(i2.2)') k
     264         CALL histdef(nid_day,"tgi"//str2,"Gas Opa IR",
     265     .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     266     .                "ave(X)",zsto1,zout)
     267         ENDDO
     268
    269269c --------------
    270270c ----- EXTINCTION GAZ
    271271         DO k=7,NSPECV,10
    272            write(str1,'(i2.2)') k
    273          CALL histdef(nid_day,"kgv"//str1,"Gas ext Vis ",
    274      .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
    275      .                "ave(X)",zsto1,zout)
    276          ENDDO
    277 c
     272           write(str2,'(i2.2)') k
     273         CALL histdef(nid_day,"kgv"//str2,"Gas ext Vis ",
     274     .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     275     .                "ave(X)",zsto1,zout)
     276         ENDDO
     277
    278278         DO k=8,NSPECI,10
    279            write(str1,'(i2.2)') k
    280          CALL histdef(nid_day,"kgi"//str1,"Gas ext IR ",
    281      .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
    282      .                "ave(X)",zsto1,zout)
    283          ENDDO
    284 c
     279           write(str2,'(i2.2)') k
     280         CALL histdef(nid_day,"kgi"//str2,"Gas ext IR ",
     281     .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     282     .                "ave(X)",zsto1,zout)
     283         ENDDO
     284
    285285c --------------
    286286c ----- OPACITE NUAGES
    287287         if (clouds.eq.1) then
    288288           CALL histdef(nid_day,"tcld","Cld Opa proxy",
    289      .                "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     289     .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
    290290     .                "ave(X)",zsto,zout)
    291 c
     291
    292292c --------------
    293293c ----- EXTINCTION NUAGES
    294294           CALL histdef(nid_day,"kcld","Cld Ext proxy",
    295      .                "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
     295     .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
    296296     .                "ave(X)",zsto,zout)
    297297         endif
    298 c
     298
    299299      ENDIF !lev_histday.GE.3
    300 c
     300
    301301c-------------------------------------------------------
    302302      IF(lev_histday.GE.4) THEN
    303 c
     303
    304304         CALL histdef(nid_day, "dtdyn", "Dynamics dT", "K/s",
    305      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    306      .                "ave(X)", zsto,zout)
    307 c
     305     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     306     .                "ave(X)", zsto,zout)
     307
    308308         CALL histdef(nid_day, "dtphy", "Physics dT", "K/s",
    309      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    310      .                "ave(X)", zsto,zout)
    311 c
     309     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     310     .                "ave(X)", zsto,zout)
     311
    312312         CALL histdef(nid_day, "dtvdf", "Boundary-layer dT", "K/s",
    313      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    314      .                "ave(X)", zsto,zout)
    315 c
     313     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     314     .                "ave(X)", zsto,zout)
     315
    316316         CALL histdef(nid_day, "dtajs", "Dry adjust. dT", "K/s",
    317      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    318      .                "ave(X)", zsto,zout)
    319 c
     317     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     318     .                "ave(X)", zsto,zout)
     319
    320320         CALL histdef(nid_day, "dtswr", "SW radiation dT", "K/s",
    321      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    322      .                "ave(X)", zsto,zout)
    323 c
     321     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     322     .                "ave(X)", zsto,zout)
     323
    324324         CALL histdef(nid_day, "dtlwr", "LW radiation dT", "K/s",
    325      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    326      .                "ave(X)", zsto,zout)
    327 c
     325     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     326     .                "ave(X)", zsto,zout)
     327
    328328c        CALL histdef(nid_day, "dtec", "Cinetic dissip dT", "K/s",
    329 c    .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     329c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    330330c    .                "ave(X)", zsto,zout)
    331 c
     331
    332332      ENDIF !lev_histday.GE.4
    333 c
     333
    334334c-------------------------------------------------------
    335335      IF(lev_histday.GE.5) THEN
    336 c
    337 c
     336
     337
    338338c        call histdef(nid_day, "taux",
    339339c    $         "Zonal wind stress", "Pa", 
    340 c    $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     340c    $         iim,jj_nb,nhori, 1,1,1, nvert, 32,
    341341c    $         "ave(X)", zsto,zout)
    342 c
     342
    343343c        call histdef(nid_day, "tauy",
    344344c    $         "Meridional xind stress", "Pa", 
    345 c    $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     345c    $         iim,jj_nb,nhori, 1,1,1, nvert, 32,
    346346c    $         "ave(X)", zsto,zout)
    347 c
     347
    348348c        CALL histdef(nid_day, "cdrm", "Momentum drag coef.", "-",
    349 c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     349c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    350350c    .                "ave(X)", zsto,zout)
    351 c
     351
    352352c        CALL histdef(nid_day, "cdrh", "Heat drag coef.", "-",
    353 c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     353c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
    354354c    .                "ave(X)", zsto,zout)
    355 c
     355
    356356      ENDIF !lev_histday.GE.5
    357357c-------------------------------------------------------
    358 c
     358
    359359         CALL histend(nid_day)
    360 c
    361          ndex2d = 0
    362          ndex3d = 0
    363 c
     360
    364361      ENDIF ! fin de test sur ok_journe
Note: See TracChangeset for help on using the changeset viewer.