Ignore:
Timestamp:
Jun 4, 2007, 4:34:47 PM (17 years ago)
Author:
Laurent Fairhead
Message:

Merge entre la version V3_conv et le HEAD
YM, JG, LF

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/trunk/libf/phylmd/initphysto.F

    r541 r766  
    77     .  (infile,
    88     .  rlon, rlat, tstep,t_ops,t_wrt,nq,fileid)
    9 
     9       
     10       USE dimphy
    1011       USE IOIPSL
    11 
     12       USE iophy
    1213      implicit none
    1314
    … …  
    5152#include "indicesol.h"
    5253#include "control.h"
    53 #include "dimphy.h"
     54cym#include "dimphy.h"
    5455
    5556C   Arguments
    … …  
    9192      tau0 = 0
    9293       
    93         CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlon,zx_lon)
    94          DO i = 1, iim
    95             zx_lon(i,1) = rlon(i+1)
    96             zx_lon(i,jjm+1) = rlon(i+1)
    97          ENDDO
    98          CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlat,zx_lat)
    99 
    100 
    101       call histbeg(infile, iim, zx_lon(:,1), jjm+1, zx_lat(1,:),
    102      .             1, iim, 1, jjm+1,
    103      .             tau0, zjulian, tstep, nhoriid, fileid)
    104 
     94cym     CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlon,zx_lon)
     95cym         DO i = 1, iim
     96cym            zx_lon(i,1) = rlon(i+1)
     97cym            zx_lon(i,jjm+1) = rlon(i+1)
     98cym         ENDDO
     99cym         CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlat,zx_lat)
     100
     101
     102      call histbeg_phy(infile,tau0, zjulian, tstep,
     103     .                 nhoriid, fileid)
     104
     105c$OMP MASTER   
    105106C  Appel a histvert pour la grille verticale
    106107C
    … …  
    108109            nivsigs(l)=float(l)
    109110         ENDDO
    110        
     111
    111112        write(*,*) 'avant histvert ds initphysto'
    112113
    … …  
    120121
    121122       CALL histdef(fileid, "phis", "Surface geop. height", "-",
    122      .                iim,jjm+1,nhoriid, 1,1,1, -99, 32,
     123     .                iim,jjphy_nb,nhoriid, 1,1,1, -99, 32,
    123124     .                "once", t_ops, t_wrt)
    124125c
    … …  
    126127
    127128         CALL histdef(fileid, "aire", "Grid area", "-",
    128      .                iim,jjm+1,nhoriid, 1,1,1, -99, 32,
     129     .                iim,jjphy_nb,nhoriid, 1,1,1, -99, 32,
    129130     .                "once", t_ops, t_wrt)
    130131         write(*,*) 'apres aire ds initphysto'
    131132
     133cym     Attention dtime et istphy ne sont pas �rit ---> a �iminer ?
    132134         CALL histdef(fileid, "dtime", "tps phys ", "s",
    133135     .                1,1,nhoriid, 1,1,1, -99, 32,
    … …  
    141143C
    142144      call histdef(fileid, 't', 'Temperature', 'K',
    143      .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
     145     .             iim, jjphy_nb, nhoriid, llm, 1, llm, zvertiid,
    144146     .             32, 'inst(X)', t_ops, t_wrt)
    145147        write(*,*) 'apres t ds initphysto'
    … …  
    147149C
    148150      call histdef(fileid, 'mfu', 'flx m. pan. mt', 'kg m/s',
    149      .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
     151     .             iim, jjphy_nb, nhoriid, llm, 1, llm, zvertiid,
    150152     .             32, 'inst(X)', t_ops, t_wrt)
    151153        write(*,*) 'apres mfu ds initphysto'
    … …  
    154156C
    155157      call histdef(fileid, 'mfd', 'flx m. pan. des', 'kg m/s',
    156      .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
     158     .             iim, jjphy_nb, nhoriid, llm, 1, llm, zvertiid,
    157159     .             32, 'inst(X)', t_ops, t_wrt)
    158160
    … …  
    161163C
    162164      call histdef(fileid, 'en_u', 'flx ent pan mt', 'kg m/s',
    163      .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
     165     .             iim, jjphy_nb, nhoriid, llm, 1, llm, zvertiid,
    164166     .             32, 'inst(X)', t_ops, t_wrt)
    165167               write(*,*) 'apres en_u ds initphysto'
    … …  
    168170C
    169171      call histdef(fileid, 'de_u', 'flx det pan mt', 'kg m/s',
    170      .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
     172     .             iim, jjphy_nb, nhoriid, llm, 1, llm, zvertiid,
    171173     .             32, 'inst(X)', t_ops, t_wrt)
    172174
    … …  
    175177C
    176178      call histdef(fileid, 'en_d', 'flx ent pan dt', 'kg m/s',
    177      .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
     179     .             iim, jjphy_nb, nhoriid, llm, 1, llm, zvertiid,
    178180     .             32, 'inst(X)', t_ops, t_wrt)
    179181C
    … …  
    183185C
    184186      call histdef(fileid, 'de_d', 'flx det pan dt', 'kg m/s',
    185      .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
     187     .             iim, jjphy_nb, nhoriid, llm, 1, llm, zvertiid,
    186188     .             32, 'inst(X)', t_ops, t_wrt)
    187189
    … …  
    189191       
    190192        call histdef(fileid, "coefh", " ", " ",
    191      .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
     193     .             iim, jjphy_nb, nhoriid, llm, 1, llm, zvertiid,
    192194     .             32, "inst(X)", t_ops, t_wrt)
    193195
    194196c abderrahmane le 16 09 02
    195197        call histdef(fileid, "fm_th", " ", " ",
    196      .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
     198     .             iim, jjphy_nb, nhoriid, llm, 1, llm, zvertiid,
    197199     .             32, "inst(X)", t_ops, t_wrt)
    198200
    199201        call histdef(fileid, "en_th", " ", " ",
    200      .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
     202     .             iim, jjphy_nb, nhoriid, llm, 1, llm, zvertiid,
    201203     .             32, "inst(X)", t_ops, t_wrt)
    202204c fin aj
    … …  
    205207
    206208        call histdef(fileid, 'frac_impa', ' ', ' ',
    207      .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
     209     .             iim, jjphy_nb, nhoriid, llm, 1, llm, zvertiid,
    208210     .             32, 'inst(X)', t_ops, t_wrt)
    209211       
    210212        call histdef(fileid, 'frac_nucl', ' ', ' ',
    211      .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
     213     .             iim, jjphy_nb, nhoriid, llm, 1, llm, zvertiid,
    212214     .             32, 'inst(X)', t_ops, t_wrt)
    213215
    … …  
    216218c
    217219      CALL histdef(fileid, "pyu1", " ", " ",
    218      .                iim,jjm+1,nhoriid, 1,1,1, -99, 32,
     220     .                iim,jjphy_nb,nhoriid, 1,1,1, -99, 32,
    219221     .                "inst(X)", t_ops, t_wrt)
    220222
    … …  
    223225c
    224226        CALL histdef(fileid, "pyv1", " ", " ",
    225      .                iim,jjm+1,nhoriid, 1,1,1, -99, 32,
     227     .                iim,jjphy_nb,nhoriid, 1,1,1, -99, 32,
    226228     .                "inst(X)", t_ops, t_wrt)
    227229       
    … …  
    231233c
    232234        call histdef(fileid, "ftsol1", " ", " ",
    233      .             iim, jjm+1, nhoriid, 1, 1,1, -99,32,
     235     .             iim, jjphy_nb, nhoriid, 1, 1,1, -99,32,
    234236     .             "inst(X)", t_ops, t_wrt)
    235237
    … …  
    238240c
    239241        call histdef(fileid, "ftsol2", " ", " ",
    240      .             iim, jjm+1, nhoriid, 1, 1,1, -99,32,
     242     .             iim, jjphy_nb, nhoriid, 1, 1,1, -99,32,
    241243     .             "inst(X)", t_ops, t_wrt)
    242244
    … …  
    245247c
    246248        call histdef(fileid, "ftsol3", " ", " ",
    247      .             iim, jjm+1, nhoriid, 1, 1,1, -99,
     249     .             iim, jjphy_nb, nhoriid, 1, 1,1, -99,
    248250     .             32, "inst(X)", t_ops, t_wrt)
    249251
    … …  
    252254c
    253255        call histdef(fileid, "ftsol4", " ", " ",
    254      .             iim, jjm+1, nhoriid, 1, 1,1, -99,
     256     .             iim, jjphy_nb, nhoriid, 1, 1,1, -99,
    255257     .             32, "inst(X)", t_ops, t_wrt)
    256258       
    … …  
    259261c
    260262        call histdef(fileid, "rain", " ", " ",
    261      .             iim, jjm+1, nhoriid, 1, 1,1, -99,
     263     .             iim, jjphy_nb, nhoriid, 1, 1,1, -99,
    262264     .             32, "inst(X)", t_ops, t_wrt)
    263265
    … …  
    266268c
    267269        call histdef(fileid, "psrf1", " ", " ",
    268      .             iim, jjm+1, nhoriid, 1, 1, 1, -99,
     270     .             iim, jjphy_nb, nhoriid, 1, 1, 1, -99,
    269271     .             32, "inst(X)", t_ops, t_wrt)
    270272       
    … …  
    273275c
    274276        call histdef(fileid, "psrf2", " ", " ",
    275      .             iim, jjm+1, nhoriid, 1, 1, 1, -99,
     277     .             iim, jjphy_nb, nhoriid, 1, 1, 1, -99,
    276278     .             32, "inst(X)", t_ops, t_wrt)
    277279
    … …  
    280282c
    281283        call histdef(fileid, "psrf3", " ", " ",
    282      .             iim, jjm+1, nhoriid, 1, 1, 1, -99,
     284     .             iim, jjphy_nb, nhoriid, 1, 1, 1, -99,
    283285     .             32, "inst(X)", t_ops, t_wrt)
    284286
    … …  
    287289c
    288290        call histdef(fileid, "psrf4", " ", " ",
    289      .             iim, jjm+1, nhoriid, 1, 1, 1, -99,
     291     .             iim, jjphy_nb, nhoriid, 1, 1, 1, -99,
    290292     .             32, "inst(X)", t_ops, t_wrt)
    291293       
    … …  
    295297c     if (ok_sync) call histsync(fileid)
    296298      if (ok_sync) call histsync
    297 
     299c$OMP END MASTER
    298300       
    299301
Note: See TracChangeset for help on using the changeset viewer.