Ignore:
Timestamp:
Jun 4, 2007, 4:34:47 PM (18 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/ini_histISCCP.h

    r684 r766  
    44      IF (ok_isccp) THEN
    55c
     6c$OMP MASTER
    67      ndex2d = 0
    78      ndex3d = 0
    … …  
    5354c definition coordonnees lon,lat en globale
    5455c
    55         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
    56         DO i = 1, iim
    57           zx_lon(i,1) = rlon(i+1)
    58           zx_lon(i,jjmp1) = rlon(i+1)
    59         ENDDO
    60 
    61         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
     56cym        CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
     57cym        DO i = 1, iim
     58cym          zx_lon(i,1) = rlon(i+1)
     59cym          zx_lon(i,jjmp1) = rlon(i+1)
     60cym        ENDDO
     61
     62cym        CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
    6263c
    6364cIM BEG region
     65cym Desole dans un premier temps le mode region ne marchera pas
     66cym Il faudra voir dans un second temps pour l'implementer
     67cym Mais cela posera des problemes au niveau de la reconstruction
     68
    6469          imin_ins=1
    6570          imax_ins=iim
    6671          jmin_ins=1
    6772          jmax_ins=jjmp1
    68           do i=1,iim-1
    69              if(zx_lon(i,1).lt.lonmin_ins) imin_ins=i
    70              if(zx_lon(i,1).le.lonmax_ins) imax_ins=i+1
    71           enddo
    72           do j=1,jjmp1
    73              if(zx_lat(1,j).ge.latmin_ins) jmax_ins=j
    74              if(zx_lat(1,j).gt.latmax_ins) jmin_ins=j
    75           enddo
     73cym          do i=1,iim-1
     74cym             if(zx_lon(i,1).lt.lonmin_ins) imin_ins=i
     75cym             if(zx_lon(i,1).le.lonmax_ins) imax_ins=i+1
     76cym          enddo
     77cym          do j=1,jjmp1
     78cym             if(zx_lat(1,j).ge.latmin_ins) jmax_ins=j
     79cym             if(zx_lat(1,j).gt.latmax_ins) jmin_ins=j
     80cym          enddo
    7681c
    7782          print*,'On stoke le fichier histISCCP instantanne sur ',
    7883     s   imin_ins,imax_ins,jmin_ins,jmax_ins
    79           print*,'On stoke le fichier histISCCP instantanne sur ',
    80      s   zx_lon(imin_ins,1),zx_lon(imax_ins,1),
    81      s   zx_lat(1,jmin_ins),zx_lat(1,jmax_ins)
     84cym          print*,'On stoke le fichier histISCCP instantanne sur ',
     85cym     s   zx_lon(imin_ins,1),zx_lon(imax_ins,1),
     86cym     s   zx_lat(1,jmin_ins),zx_lat(1,jmax_ins)
    8287cIM END region
    8388c
    8489        IF(1.EQ.0) THEN
    85          CALL histbeg("histISCCP.nc", iim,zx_lon(:,1),jjmp1,zx_lat(1,:),
    86      .                 1, iim, 1, jjmp1,
    87      .                 itau_phy, zjulian, dtime,
     90cym         CALL histbeg("histISCCP.nc", iim,zx_lon(:,1),jjmp1,zx_lat(1,:),
     91cym     .                 1, iim, 1, jjmp1,
     92cym     .                 itau_phy, zjulian, dtime,
     93cym     .                 nhori, nid_isccp)
     94         CALL histbeg_phy("histISCCP.nc", itau_phy, zjulian, dtime,
    8895     .                 nhori, nid_isccp)
     96
    8997        ENDIF !(1.EQ.0) THEN
    9098c
    91          CALL histbeg("histISCCP.nc", iim,zx_lon(:,1),
    92      .                 jjmp1,zx_lat(1,:),
    93      .                 imin_ins,imax_ins-imin_ins+1,
    94      .                 jmin_ins,jmax_ins-jmin_ins+1,
    95      .                 itau_phy, zjulian, dtime,
     99cym         CALL histbeg("histISCCP.nc", iim,zx_lon(:,1),
     100cym     .                 jjmp1,zx_lat(1,:),
     101cym     .                 imin_ins,imax_ins-imin_ins+1,
     102cym     .                 jmin_ins,jmax_ins-jmin_ins+1,
     103cym     .                 itau_phy, zjulian, dtime,
     104cym     .                 nhori, nid_isccp)
     105
     106         CALL histbeg_phy("histISCCP.nc", itau_phy, zjulian, dtime,
    96107     .                 nhori, nid_isccp)
     108
    97109c
    98110        IF(type_run.EQ."ENSP".OR.type_run.EQ."CLIM") THEN
    … …  
    118130          CALL histdef(nid_isccp, "cldISCCP_"//taulev(k)//typinout(n),
    119131     .                "LMDZ ISCCP cld", "%",
    120      .                iim, jjmp1,nhori,lmaxm1,1,lmaxm1,nvert,32,
     132     .                iim, jjphy_nb,nhori,lmaxm1,1,lmaxm1,nvert,32,
    121133     .                "ave(X)", zcals(n),zout_isccp(n))
    122134         ENDDO
    … …  
    124136         CALL histdef(nid_isccp, "nsunlit"//typinout(n),
    125137     .                "Nb of calls with sunlit ", "%",
    126      .                iim, jjmp1,nhori,1,1,1,-99,32,
     138     .                iim, jjphy_nb,nhori,1,1,1,-99,32,
    127139     .                "ave(X)", zcals(n),zout_isccp(n))
    128140c
    … …  
    140152           CALL histdef(nid_isccp, pclev(l)//taulev(k)//typinout(n),
    141153     .                "LMDZ ISCCP cld "//cnameisccp(l,k), "%",
    142      .                iim, jjmp1,nhori,1,1,1,-99,32,
     154     .                iim, jjphy_nb,nhori,1,1,1,-99,32,
    143155     .                "ave(X)", zcals(n),zout_isccp(n))
    144156c
    … …  
    149161          CALL histdef(nid_isccp, "nsunlit"//typinout(n),
    150162     .                "Nb of calls with sunlit ", "%",
    151      .                iim, jjmp1,nhori,1,1,1,-99,32,
     163     .                iim, jjphy_nb,nhori,1,1,1,-99,32,
    152164     .                "ave(X)", zcals(n),zout_isccp(n))
    153165c
    … …  
    155167          CALL histdef(nid_isccp, "cirr",
    156168     .                "Cirrus lk ISCCP-D2", "%",
    157      .                iim, jjmp1,nhori,1,1,1,-99,32,
     169     .                iim, jjphy_nb,nhori,1,1,1,-99,32,
    158170     .                "ave(X)", zcals(n),zout_isccp(n))
    159171c
    160172          CALL histdef(nid_isccp, "cist",
    161173     .                "CiSt lk ISCCP-D2", "%",
    162      .                iim, jjmp1,nhori,1,1,1,-99,32,
     174     .                iim, jjphy_nb,nhori,1,1,1,-99,32,
    163175     .                "ave(X)", zcals(n),zout_isccp(n))
    164176c
    165177          CALL histdef(nid_isccp, "deep",
    166178     .                "Deep lk ISCCP-D2", "%",
    167      .                iim, jjmp1,nhori,1,1,1,-99,32,
     179     .                iim, jjphy_nb,nhori,1,1,1,-99,32,
    168180     .                "ave(X)", zcals(n),zout_isccp(n))
    169181c
    170182          CALL histdef(nid_isccp, "alcu",
    171183     .                "AlCu lk ISCCP-D2", "%",
    172      .                iim, jjmp1,nhori,1,1,1,-99,32,
     184     .                iim, jjphy_nb,nhori,1,1,1,-99,32,
    173185     .                "ave(X)", zcals(n),zout_isccp(n))
    174186c
    175187          CALL histdef(nid_isccp, "alst",
    176188     .                "AlSt lk ISCCP-D2", "%",
    177      .                iim, jjmp1,nhori,1,1,1,-99,32,
     189     .                iim, jjphy_nb,nhori,1,1,1,-99,32,
    178190     .                "ave(X)", zcals(n),zout_isccp(n))
    179191c
    180192          CALL histdef(nid_isccp, "nist",
    181193     .                "NiSt lk ISCCP-D2", "%",
    182      .                iim, jjmp1,nhori,1,1,1,-99,32,
     194     .                iim, jjphy_nb,nhori,1,1,1,-99,32,
    183195     .                "ave(X)", zcals(n),zout_isccp(n))
    184196c
    185197          CALL histdef(nid_isccp, "cumu",
    186198     .                "Cumu lk ISCCP-D2", "%",
    187      .                iim, jjmp1,nhori,1,1,1,-99,32,
     199     .                iim, jjphy_nb,nhori,1,1,1,-99,32,
    188200     .                "ave(X)", zcals(n),zout_isccp(n))
    189201c
    190202          CALL histdef(nid_isccp, "stcu",
    191203     .                "StCu lk ISCCP-D2", "%",
    192      .                iim, jjmp1,nhori,1,1,1,-99,32,
     204     .                iim, jjphy_nb,nhori,1,1,1,-99,32,
    193205     .                "ave(X)", zcals(n),zout_isccp(n))
    194206c
    195207          CALL histdef(nid_isccp, "stra",
    196208     .                "Stra lk ISCCP-D2", "%",
    197      .                iim, jjmp1,nhori,1,1,1,-99,32,
     209     .                iim, jjphy_nb,nhori,1,1,1,-99,32,
    198210     .                "ave(X)", zcals(n),zout_isccp(n))
    199211c
    … …  
    201213          CALL histdef(nid_isccp, "thin",
    202214     .                "Opt. thin ISCCP-D2 like clouds", "%",
    203      .                iim, jjmp1,nhori,lmax3,1,lmax3,nvert3,32,
     215     .                iim, jjphy_nb,nhori,lmax3,1,lmax3,nvert3,32,
    204216     .                "ave(X)", zcals(n),zout_isccp(n))
    205217c
    206218          CALL histdef(nid_isccp, "mid",
    207219     .                "Opt. intermediate ISCCP-D2 like clouds", "%",
    208      .                iim, jjmp1,nhori,lmax3,1,lmax3,nvert3,32,
     220     .                iim, jjphy_nb,nhori,lmax3,1,lmax3,nvert3,32,
    209221     .                "ave(X)", zcals(n),zout_isccp(n))
    210222c
    211223          CALL histdef(nid_isccp, "thick",
    212224     .                "Opt. thick ISCCP-D2 like clouds", "%",
    213      .                iim, jjmp1,nhori,lmax3,1,lmax3,nvert3,32,
     225     .                iim, jjphy_nb,nhori,lmax3,1,lmax3,nvert3,32,
    214226     .                "ave(X)", zcals(n),zout_isccp(n))
    215227c
    … …  
    224236          CALL histdef(nid_isccp,"boxptop"//verticaxe(n),
    225237     .                "Boxptop axe"//verticaxe(n), "mb",
    226      .                iim, jjmp1,nhori,
     238     .                iim, jjphy_nb,nhori,
    227239     .                ncol(n),1,ncol(n),nvlev(n),32,
    228240cIM  .                ncolmx,1,ncolmx,nvlev,32,
    … …  
    237249          CALL histdef(nid_isccp, "seed"//verticaxe(n),
    238250     .                "seed axe"//verticaxe(n), "-",
    239      .                iim, jjmp1,nhori,1,1,1,-99,32,
     251     .                iim, jjphy_nb,nhori,1,1,1,-99,32,
    240252cIM  .                "inst(X)", dtime,dtime)
    241253     .                "ave(X)", zcals(n),zout_isccp(n))
    … …  
    245257        CALL histend(nid_isccp)
    246258c
     259c$OMP END MASTER
    247260      ENDIF ! ok_isccp
Note: See TracChangeset for help on using the changeset viewer.