Changeset 590


Ignore:
Timestamp:
Feb 7, 2005, 5:41:35 PM (19 years ago)
Author:
Laurent Fairhead
Message:

Synchro avec LOOP PC
LF

Location:
LMDZ4/trunk/libf/phylmd
Files:
6 edited

Legend:

Unmodified
Added
Removed

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

    r541 r590  
    241241c
    242242      REAL yt2m(klon), yq2m(klon), yu10m(klon)
     243c -- LOOP
     244       REAL yu10mx(klon)
     245       REAL yu10my(klon)
     246       REAL ywindsp(klon)
     247c -- LOOP
    243248c
    244249      REAL uzon(klon), vmer(klon)
    … …  
    351356      ytsoil = 999999.
    352357      yrugoro = 0.
    353 
     358c -- LOOP
     359      yu10mx = 0.0
     360      yu10my = 0.0
     361      ywindsp = 0.0
     362c -- LOOP
    354363      DO nsrf = 1, nbsrf
    355364      DO i = 1, klon
    … …  
    458467        ypaprs(j,klev+1) = paprs(i,klev+1)
    459468        y_run_off_lic_0(j) = run_off_lic_0(i)
     469c -- LOOP
     470       yu10mx(j) = u10m(i,nsrf)
     471       yu10my(j) = v10m(i,nsrf)
     472       ywindsp(j) = SQRT(yu10mx(j)*yu10mx(j) + yu10my(j)*yu10my(j) )
     473c -- LOOP
    460474      END DO
    461475C
    … …  
    667681     e          yt,yq,yts,ypaprs,ypplay,
    668682     e          ydelp,yrads,yalb, yalblw, ysnow, yqsurf,
    669      e          yrain_f, ysnow_f, yfder, ytaux, ytauy,
     683     e          yrain_f, ysnow_f, yfder, ytaux, ytauy,
     684c -- LOOP
     685     e          ywindsp,
     686c -- LOOP
    670687c$$$     e          ysollw, ysolsw,
    671688     e          ysollw, ysollwdown, ysolsw,yfluxlat,
    … …  
    900917     e                delp,radsol,albedo,alblw,snow,qsurf,
    901918     e                precip_rain, precip_snow, fder, taux, tauy,
     919c -- LOOP
     920     e                ywindsp,
     921c -- LOOP
    902922     $                sollw, sollwdown, swnet,fluxlat,
    903923     s                pctsrf_new, agesno,
    … …  
    958978      character*6 ocean
    959979      integer npas, nexca
     980c -- LOOP
     981       REAL yu10mx(klon)
     982       REAL yu10my(klon)
     983       REAL ywindsp(klon)
     984c -- LOOP
     985
    960986
    961987c
    … …  
    11841210     e tq_cdrag, petAcoef, peqAcoef, petBcoef, peqBcoef,
    11851211     e precip_rain, precip_snow, sollw, sollwdown, swnet, swdown,
    1186      e fder, taux, tauy, rugos, rugoro,
     1212     e fder, taux, tauy,
     1213c -- LOOP
     1214     e ywindsp,
     1215c -- LOOP
     1216     e rugos, rugoro,
    11871217     e albedo, snow, qsurf,
    11881218     e ts, p1lay, psref, radsol,

  • LMDZ4/trunk/libf/phylmd/ini_histhf.h

    r524 r590  
    113113     .                "ave(X)", zsto1,zout)
    114114c
    115 c
     115c -- LOOP
    116116         CALL histdef(nid_hf, "SWdownOR",
    117      .                "Sfce incident SW radiation OR", "W/m^2",
    118      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    119      .                "ave(X)", zsto1,zout)
     117     .                "Sfce incident SW down radiation OR", "W/m^2",
     118     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     119     .                "ave(X)", zsto1,zout)
     120c
     121         CALL histdef(nid_hf, "SWnetOR",
     122     .                "Sfce incident SW net radiation OR", "W/m^2",
     123     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     124     .                "ave(X)", zsto1,zout)
     125c -- LOOP
    120126c
    121127         CALL histdef(nid_hf, "LWdownOR",

  • LMDZ4/trunk/libf/phylmd/interface_surf.F90

    r560 r590  
    6464      & tq_cdrag, petAcoef, peqAcoef, petBcoef, peqBcoef, &
    6565      & precip_rain, precip_snow, sollw, sollwdown, swnet, swdown, &
    66       & fder, taux, tauy, rugos, rugoro, &
     66      & fder, taux, tauy, &
     67! -- LOOP
     68      & windsp, &
     69! -- LOOP
     70      & rugos, rugoro, &
    6771      & albedo, snow, qsurf, &
    6872      & tsurf, p1lay, ps, radsol, &
    … …  
    127131!   fder         derivee des flux (pour le couplage)
    128132!   taux, tauy   tension de vents
     133! -- LOOP
     134!   windsp       module du vent a 10m
     135! -- LOOP
    129136!   rugos        rugosite
    130137!   zmasq        masque terre/ocean
    … …  
    175182  real, dimension(klon), intent(IN) :: zmasq
    176183  real, dimension(klon), intent(IN) :: taux, tauy, rugos, rugoro
     184! -- LOOP
     185  real, dimension(klon), intent(IN) :: windsp
     186! -- LOOP
    177187  character (len = 6)  :: ocean
    178188  integer              :: npas, nexca ! nombre et pas de temps couplage
    … …  
    444454      & ocean, npas, nexca, debut, lafin, &
    445455      & swdown, sollw, precip_rain, precip_snow, evap, tsurf, &
    446       & fluxlat, fluxsens, fder, albedo, taux, tauy, zmasq, &
    447       & tsurf_new, alb_new, pctsrf_new)
     456      & fluxlat, fluxsens, fder, albedo, taux, tauy, &
     457! -- LOOP
     458      & windsp, &
     459! -- LOOP
     460      & zmasq, &
     461      & tsurf_new, alb_new, &
     462      & pctsrf_new)
    448463
    449464!    else if (ocean == 'slab  ') then
    … …  
    498513      & ocean, npas, nexca, debut, lafin, &
    499514      & swdown, sollw, precip_rain, precip_snow, evap, tsurf, &
    500       & fluxlat, fluxsens, fder, albedo, taux, tauy, zmasq, &
    501       & tsurf_new, alb_new, pctsrf_new)
     515      & fluxlat, fluxsens, fder, albedo, taux, tauy, &
     516! -- LOOP
     517      & windsp, &
     518! -- LOOP
     519      & zmasq, &
     520      & tsurf_new, alb_new, &
     521      & pctsrf_new)
    502522
    503523!    else if (ocean == 'slab  ') then
    … …  
    549569      & ocean, npas, nexca, debut, lafin, &
    550570      & swdown, sollw, precip_rain, precip_snow, evap, tsurf, &
    551       & fluxlat, fluxsens, fder, albedo, taux, tauy, zmasq, &
    552       & tsurf_new, alb_new, pctsrf_new)
     571      & fluxlat, fluxsens, fder, albedo, taux, tauy, &
     572! -- LOOP
     573      & windsp, &
     574! -- LOOP
     575      & zmasq, &
     576      & tsurf_new, alb_new, &
     577      & pctsrf_new)
    553578
    554579      tsurf_temp = tsurf_new
    … …  
    644669      & ocean, npas, nexca, debut, lafin, &
    645670      & swdown, sollw, precip_rain, precip_snow, evap, tsurf, &
    646       & fluxlat, fluxsens, fder, albedo, taux, tauy, zmasq, &
    647       & tsurf_new, alb_new, pctsrf_new)
     671      & fluxlat, fluxsens, fder, albedo, taux, tauy, &
     672! -- LOOP
     673      & windsp, &
     674! -- LOOP
     675      & zmasq, &
     676      & tsurf_new, alb_new, &
     677      & pctsrf_new)
    648678
    649679!    else if (ocean == 'slab  ') then
    … …  
    11981228      & ocean, npas, nexca, debut, lafin, &
    11991229      & swdown, lwdown, precip_rain, precip_snow, evap, tsurf, &
    1200       & fluxlat, fluxsens, fder, albsol, taux, tauy, zmasq, &
    1201       & tsurf_new, alb_new, pctsrf_new)
     1230      & fluxlat, fluxsens, fder, albsol, taux, tauy, &
     1231! -- LOOP
     1232      & windsp, &
     1233! -- LOOP
     1234      & zmasq, &
     1235      & tsurf_new, alb_new, &
     1236      & pctsrf_new)
    12021237
    12031238! Cette routine sert d'interface entre le modele atmospherique et un
    … …  
    12421277!   taux         tension de vent en x
    12431278!   tauy         tension de vent en y
     1279! -- LOOP
     1280!    windsp       module du vent a 10m
     1281! -- LOOP
    12441282!   nexca        frequence de couplage
    12451283!   zmasq        masque terre/ocean
    … …  
    12691307  real, dimension(klon), intent(IN) :: precip_rain, precip_snow
    12701308  real, dimension(klon), intent(IN) :: tsurf, fder, albsol, taux, tauy
     1309! -- LOOP
     1310   real, dimension(klon), intent(IN) :: windsp
     1311! -- LOOP
    12711312  INTEGER              :: nexca, npas, kstep
    12721313  real, dimension(klon), intent(IN) :: zmasq
    … …  
    12811322! Variables locales
    12821323  integer                    :: j, error, sum_error, ig, cpl_index,i
     1324! -- LOOP
     1325  INTEGER :: nsrf
     1326! -- LOOP
    12831327  character (len = 20) :: modname = 'interfoce_cpl'
    12841328  character (len = 80) :: abort_message
    … …  
    12881332  real, allocatable, dimension(:,:),save :: cpl_snow, cpl_evap, cpl_tsol
    12891333  real, allocatable, dimension(:,:),save :: cpl_fder, cpl_albe, cpl_taux
     1334! -- LOOP
     1335  real, allocatable, dimension(:,:),save :: cpl_windsp
     1336! -- LOOP
    12901337  real, allocatable, dimension(:,:),save :: cpl_tauy
    12911338  REAL, ALLOCATABLE, DIMENSION(:,:),SAVE :: cpl_rriv, cpl_rcoa, cpl_rlic
    … …  
    12951342  real, allocatable, dimension(:,:,:),save :: tmp_snow, tmp_evap, tmp_tsol
    12961343  real, allocatable, dimension(:,:,:),save :: tmp_fder, tmp_albe, tmp_taux
     1344! -- LOOP
     1345  real, allocatable, dimension(:,:,:),save :: tmp_windsp
     1346! -- LOOP
    12971347!!$  real, allocatable, dimension(:,:,:),save :: tmp_tauy, tmp_rriv, tmp_rcoa
    12981348  REAL, ALLOCATABLE, DIMENSION(:,:,:),SAVE :: tmp_tauy
    … …  
    13021352  REAL, DIMENSION(iim, jjm+1) :: wri_evap_sea, wri_rcoa, wri_rriv
    13031353  REAL, DIMENSION(iim, jjm+1) :: wri_rain, wri_snow, wri_taux, wri_tauy
     1354! -- LOOP
     1355   REAL, DIMENSION(iim, jjm+1) :: wri_windsp
     1356! -- LOOP
    13041357  REAL, DIMENSION(iim, jjm+1) :: wri_calv
    13051358  REAL, DIMENSION(iim, jjm+1) :: wri_tauxx, wri_tauyy, wri_tauzz
    … …  
    13341387!med  integer :: itau_w
    13351388  integer,save :: itau_w
     1389! -- LOOP
     1390  integer :: nb_interf_cpl
     1391! -- LOOP
    13361392#include "param_cou.h"
    13371393#include "inc_cpl.h"
    … …  
    13691425    allocate(cpl_albe(klon,2), stat = error); sum_error = sum_error + error
    13701426    allocate(cpl_taux(klon,2), stat = error); sum_error = sum_error + error
     1427! -- LOOP
     1428     allocate(cpl_windsp(klon,2), stat = error); sum_error = sum_error + error
     1429! -- LOOP
    13711430    allocate(cpl_tauy(klon,2), stat = error); sum_error = sum_error + error
    13721431    ALLOCATE(cpl_rriv(iim,jjm+1), stat=error); sum_error = sum_error + error
    … …  
    13861445    cpl_evap = 0.; cpl_tsol = 0.; cpl_fder = 0.; cpl_albe = 0.
    13871446    cpl_taux = 0.; cpl_tauy = 0.; cpl_rriv = 0.; cpl_rcoa = 0.; cpl_rlic = 0.
     1447! -- LOOP
     1448     cpl_windsp = 0.
     1449! -- LOOP
    13881450
    13891451    sum_error = 0
    … …  
    14651527
    14661528! calcul des fluxs a passer
    1467 
     1529! -- LOOP
     1530  nb_interf_cpl = nb_interf_cpl + 1 
     1531  if (check) write(lunout,*)'passage dans interface_surf.F90 :  ',nb_interf_cpl
     1532! -- LOOP
    14681533  cpl_index = 1
    14691534  if (nisurf == is_sic) cpl_index = 2
    14701535  if (cumul) then
     1536! -- LOOP
     1537      if (check) write(lunout,*)'passage dans cumul '
     1538      if (check) write(lunout,*)'valeur de cpl_index ', cpl_index
     1539! -- LOOP 
    14711540    if (check) write(*,*) modname, 'cumul des champs'
    14721541    do ig = 1, knon
    … …  
    14921561      cpl_tauy(ig,cpl_index) = cpl_tauy(ig,cpl_index) &
    14931562       &                          + tauy(ig)        / FLOAT(nexca)
     1563! -- LOOP
     1564      IF (cpl_index .EQ. 1) THEN
     1565      cpl_windsp(ig,cpl_index) = cpl_windsp(ig,cpl_index) &
     1566       &                          + windsp(ig)      / FLOAT(nexca)
     1567      ENDIF
     1568! -- LOOP
    14941569    enddo
    14951570    IF (cpl_index .EQ. 1) THEN
    … …  
    16031678      allocate(tmp_taux(iim,jjm+1,2), stat=error); sum_error = sum_error + error
    16041679      allocate(tmp_tauy(iim,jjm+1,2), stat=error); sum_error = sum_error + error
     1680! -- LOOP
     1681       allocate(tmp_windsp(iim,jjm+1,2), stat=error); sum_error = sum_error + error
     1682! -- LOOP
    16051683!!$      allocate(tmp_rriv(iim,jjm+1,2), stat=error); sum_error = sum_error + error
    16061684!!$      allocate(tmp_rcoa(iim,jjm+1,2), stat=error); sum_error = sum_error + error
    … …  
    16251703    call gath2cpl(cpl_albe(1,cpl_index), tmp_albe(1,1,cpl_index), klon, knon,iim,jjm,                  knindex)
    16261704    call gath2cpl(cpl_taux(1,cpl_index), tmp_taux(1,1,cpl_index), klon, knon,iim,jjm,                  knindex)
     1705! -- LOOP
     1706     call gath2cpl(cpl_windsp(1,cpl_index), tmp_windsp(1,1,cpl_index), klon, knon,iim,jjm,             knindex)
     1707! -- LOOP
    16271708    call gath2cpl(cpl_tauy(1,cpl_index), tmp_tauy(1,1,cpl_index), klon, knon,iim,jjm,                  knindex)
    16281709
    … …  
    16331714      wri_rain = 0.; wri_snow = 0.; wri_rcoa = 0.; wri_rriv = 0.
    16341715      wri_taux = 0.; wri_tauy = 0.
     1716! -- LOOP
     1717       wri_windsp = 0.
     1718! -- LOOP     
    16351719      call gath2cpl(pctsrf(1,is_oce), tamp_srf(1,1,1), klon, klon, iim, jjm, tamp_ind)
    16361720      call gath2cpl(pctsrf(1,is_sic), tamp_srf(1,1,2), klon, klon, iim, jjm, tamp_ind)
    … …  
    16431727      wri_evap_ice = tmp_evap(:,:,2)
    16441728      wri_evap_sea = tmp_evap(:,:,1)
     1729! -- LOOP
     1730       wri_windsp = tmp_windsp(:,:,1)
     1731! -- LOOP
     1732
    16451733!!$PB
    16461734      wri_rriv = cpl_rriv(:,:)
    … …  
    16961784! envoi au coupleur
    16971785!
    1698       CALL histwrite(nidct,cl_writ(1),itau_w,wri_sol_ice,iim*(jjm+1),ndexct)
    1699       CALL histwrite(nidct,cl_writ(2),itau_w,wri_sol_sea,iim*(jjm+1),ndexct)
    1700       CALL histwrite(nidct,cl_writ(3),itau_w,wri_nsol_ice,iim*(jjm+1),ndexct)
    1701       CALL histwrite(nidct,cl_writ(4),itau_w,wri_nsol_sea,iim*(jjm+1),ndexct)
    1702       CALL histwrite(nidct,cl_writ(5),itau_w,wri_fder_ice,iim*(jjm+1),ndexct)
    1703       CALL histwrite(nidct,cl_writ(6),itau_w,wri_evap_ice,iim*(jjm+1),ndexct)
    1704       CALL histwrite(nidct,cl_writ(7),itau_w,wri_evap_sea,iim*(jjm+1),ndexct)
    1705       CALL histwrite(nidct,cl_writ(8),itau_w,wri_rain,iim*(jjm+1),ndexct)
    1706       CALL histwrite(nidct,cl_writ(9),itau_w,wri_snow,iim*(jjm+1),ndexct)
    1707       CALL histwrite(nidct,cl_writ(10),itau_w,wri_rcoa,iim*(jjm+1),ndexct)
    1708       CALL histwrite(nidct,cl_writ(11),itau_w,wri_rriv,iim*(jjm+1),ndexct)
    1709       CALL histwrite(nidct,cl_writ(12),itau_w,wri_calv,iim*(jjm+1),ndexct)
    1710       CALL histwrite(nidct,cl_writ(13),itau_w,wri_tauxx,iim*(jjm+1),ndexct)
    1711       CALL histwrite(nidct,cl_writ(14),itau_w,wri_tauyy,iim*(jjm+1),ndexct)
    1712       CALL histwrite(nidct,cl_writ(15),itau_w,wri_tauzz,iim*(jjm+1),ndexct)
    1713       CALL histwrite(nidct,cl_writ(16),itau_w,wri_tauxx,iim*(jjm+1),ndexct)
    1714       CALL histwrite(nidct,cl_writ(17),itau_w,wri_tauyy,iim*(jjm+1),ndexct)
    1715       CALL histwrite(nidct,cl_writ(18),itau_w,wri_tauzz,iim*(jjm+1),ndexct)
     1786      CALL histwrite(nidct,cl_writ(8),itau_w,wri_sol_ice,iim*(jjm+1),ndexct)
     1787      CALL histwrite(nidct,cl_writ(9),itau_w,wri_sol_sea,iim*(jjm+1),ndexct)
     1788      CALL histwrite(nidct,cl_writ(10),itau_w,wri_nsol_ice,iim*(jjm+1),ndexct)
     1789      CALL histwrite(nidct,cl_writ(11),itau_w,wri_nsol_sea,iim*(jjm+1),ndexct)
     1790      CALL histwrite(nidct,cl_writ(12),itau_w,wri_fder_ice,iim*(jjm+1),ndexct)
     1791      CALL histwrite(nidct,cl_writ(13),itau_w,wri_evap_ice,iim*(jjm+1),ndexct)
     1792      CALL histwrite(nidct,cl_writ(14),itau_w,wri_evap_sea,iim*(jjm+1),ndexct)
     1793      CALL histwrite(nidct,cl_writ(15),itau_w,wri_rain,iim*(jjm+1),ndexct)
     1794      CALL histwrite(nidct,cl_writ(16),itau_w,wri_snow,iim*(jjm+1),ndexct)
     1795      CALL histwrite(nidct,cl_writ(17),itau_w,wri_rcoa,iim*(jjm+1),ndexct)
     1796      CALL histwrite(nidct,cl_writ(18),itau_w,wri_rriv,iim*(jjm+1),ndexct)
     1797      CALL histwrite(nidct,cl_writ(19),itau_w,wri_calv,iim*(jjm+1),ndexct)
     1798      CALL histwrite(nidct,cl_writ(1),itau_w,wri_tauxx,iim*(jjm+1),ndexct)
     1799      CALL histwrite(nidct,cl_writ(2),itau_w,wri_tauyy,iim*(jjm+1),ndexct)
     1800      CALL histwrite(nidct,cl_writ(3),itau_w,wri_tauzz,iim*(jjm+1),ndexct)
     1801      CALL histwrite(nidct,cl_writ(4),itau_w,wri_tauxx,iim*(jjm+1),ndexct)
     1802      CALL histwrite(nidct,cl_writ(5),itau_w,wri_tauyy,iim*(jjm+1),ndexct)
     1803      CALL histwrite(nidct,cl_writ(6),itau_w,wri_tauzz,iim*(jjm+1),ndexct)
     1804! -- LOOP
     1805      CALL histwrite(nidct,cl_writ(7),itau_w,wri_windsp,iim*(jjm+1),ndexct)
     1806! -- LOOP
    17161807      CALL histsync(nidct)
    17171808! pas utile      IF (lafin) CALL histclo(nidct)
    … …  
    17231814      & wri_nsol_sea, wri_fder_ice, wri_evap_ice, wri_evap_sea, wri_rain, &
    17241815      & wri_snow, wri_rcoa, wri_rriv, wri_calv, wri_tauxx, wri_tauyy,     &
    1725       & wri_tauzz, wri_tauxx, wri_tauyy, wri_tauzz,lafin )
     1816      & wri_tauzz, wri_tauxx, wri_tauyy, wri_tauzz,                       &
     1817! -- LOOP
     1818      & wri_windsp,lafin)
     1819! -- LOOP
    17261820#else
    17271821      call intocpl(itime, (jjm+1)*iim, wri_sol_ice, wri_sol_sea, wri_nsol_ice,&
    17281822      & wri_nsol_sea, wri_fder_ice, wri_evap_ice, wri_evap_sea, wri_rain, &
    17291823      & wri_snow, wri_rcoa, wri_rriv, wri_calv, wri_tauxx, wri_tauyy,     &
    1730       & wri_tauzz, wri_tauxx, wri_tauyy, wri_tauzz,lafin )
     1824      & wri_tauzz, wri_tauxx, wri_tauyy, wri_tauzz,                       &
     1825! -- LOOP
     1826      & wri_windsp,lafin)
     1827! -- LOOP
    17311828#endif
    17321829#endif
    … …  
    17351832      cpl_evap = 0.; cpl_tsol = 0.; cpl_fder = 0.; cpl_albe = 0.
    17361833      cpl_taux = 0.; cpl_tauy = 0.; cpl_rriv = 0.; cpl_rcoa = 0.; cpl_rlic = 0.
     1834! -- LOOP
     1835      cpl_windsp = 0.
     1836! -- LOOP
    17371837!
    17381838! deallocation memoire variables temporaires
    … …  
    17491849      deallocate(tmp_taux, stat=error); sum_error = sum_error + error
    17501850      deallocate(tmp_tauy, stat=error); sum_error = sum_error + error
     1851! -- LOOP
     1852      deallocate(tmp_windsp, stat=error); sum_error = sum_error + error
     1853! -- LOOP
    17511854!!$PB
    17521855!!$      deallocate(tmp_rriv, stat=error); sum_error = sum_error + error

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

    r541 r590  
    22! $Header$
    33!
     4C $Id$
    45C****
    56C

  • LMDZ4/trunk/libf/phylmd/oasis.dummy

    r524 r590  
    2020c     INCLUDE "param.h"
    2121c
    22       INTEGER kastp, kexch, kstep,imjm
     22      INTEGER kastp, kexch, kstep,imjm,klon
    2323      INTEGER iparal(3)
    2424      INTEGER ifcpl, idt, info, imxtag, istep, jf
     25
     26c -- LOOP
     27c
     28#include "dimensions.h"
     29      INTEGER jjmp1
     30      PARAMETER (jjmp1=jjm+1-1/jjm)
     31#include "dimphy.h"
     32       REAL zwindsp(klon)
     33c
     34c -- LOOP
    2535c
    2636#include "param_cou.h"
    … …  
    7282c         must be the same as (1) of the field  definition in namcouple:
    7383c
    74       cl_writ(1)='COSHFICE'
    75       cl_writ(2)='COSHFOCE'
    76       cl_writ(3)='CONSFICE'
    77       cl_writ(4)='CONSFOCE'
    78       cl_writ(5)='CODFLXDT'
     84      cl_writ(8)='COSHFICE'
     85      cl_writ(9)='COSHFOCE'
     86      cl_writ(10)='CONSFICE'
     87      cl_writ(11)='CONSFOCE'
     88      cl_writ(12)='CODFLXDT'
    7989c      cl_writ(6)='COICTEMP'
    80       cl_writ(6)='COTFSICE'
    81       cl_writ(7)='COTFSOCE'
    82       cl_writ(8)='COTOLPSU'
    83       cl_writ(9)='COTOSPSU'
    84       cl_writ(10)='CORUNCOA'
    85       cl_writ(11)='CORIVFLU'
    86       cl_writ(12)='COCALVIN'
     90      cl_writ(13)='COTFSICE'
     91      cl_writ(14)='COTFSOCE'
     92      cl_writ(15)='COTOLPSU'
     93      cl_writ(16)='COTOSPSU'
     94      cl_writ(17)='CORUNCOA'
     95      cl_writ(18)='CORIVFLU'
     96      cl_writ(19)='COCALVIN'
    8797c$$$      cl_writ(13)='COZOTAUX'
    8898c$$$      cl_writ(14)='COZOTAUV'
    8999c$$$      cl_writ(15)='COMETAUY'
    90100c$$$      cl_writ(16)='COMETAUU'
    91       cl_writ(13)='COTAUXXU'
    92       cl_writ(14)='COTAUYYU'
    93       cl_writ(15)='COTAUZZU'
    94       cl_writ(16)='COTAUXXV'
    95       cl_writ(17)='COTAUYYV'
    96       cl_writ(18)='COTAUZZV'
     101      cl_writ(1)='COTAUXXU'
     102      cl_writ(2)='COTAUYYU'
     103      cl_writ(3)='COTAUZZU'
     104      cl_writ(4)='COTAUXXV'
     105      cl_writ(5)='COTAUYYV'
     106      cl_writ(6)='COTAUZZV'
     107c -- LOOP
     108      cl_writ(7)='COWINDSP'
     109c -- LOOP
    97110c
    98111c     Define files name for fields exchanged from atmos to coupler,
    … …  
    117130      cl_f_writ(17)='flxatmos'
    118131      cl_f_writ(18)='flxatmos'
     132c -- LOOP
     133      cl_f_writ(19)='flxatmos'
     134c -- LOOP
    119135
    120136c
    … …  
    310326
    311327c $Id$
     328c -- LOOP
    312329      SUBROUTINE intocpl(kt, imjm, fsolice, fsolwat, fnsolice, fnsolwat,
    313330     $    fnsicedt, evice, evwat, lpre, spre, dirunoff, rivrunoff,
    314      $    calving, tauxx_u, tauyy_u, tauzz_u, tauxx_v, tauyy_v, tauzz_v
    315      $    , last)
     331     $    calving, tauxx_u, tauyy_u, tauzz_u, tauxx_v, tauyy_v, tauzz_v,
     332     $    windsp, last)
     333c -- LOOP
    316334c ======================================================================
    317335c S. Valcke (02/99) adapted From L.Z.X Li: this subroutine provides the
    … …  
    321339c ======================================================================
    322340      IMPLICIT NONE
     341c -- LOOP
     342c
     343#include "dimensions.h"
     344      INTEGER jjmp1
     345      PARAMETER (jjmp1=jjm+1-1/jjm)
     346#include "dimphy.h"
     347c      REAL zu10m(klon), zv10m(klon)
     348       REAL zwindsp(klon)
     349c
     350c -- LOOP
     351c
     352
    323353      INTEGER kt, imjm
    324354c
    … …  
    345375      REAL tauzz_u(imjm)
    346376      REAL tauzz_v(imjm)
     377c -- LOOP
     378       REAL windsp(imjm)
     379c -- LOOP
    347380      LOGICAL last
    348381c
    … …  
    424457c         WRITE fields to files         
    425458          DO jf=1, jpflda2o1 + jpflda2o2
    426             IF (jf.eq.1)
     459            IF (jf.eq.8)
    427460     $          CALL locwrite(cl_writ(jf),fsolice, imjm,
    428461     $          file_unit_field(jf), ierror)
    429             IF (jf.eq.2)
     462            IF (jf.eq.9)
    430463     $          CALL locwrite(cl_writ(jf),fsolwat, imjm,
    431464     $          file_unit_field(jf), ierror)
    432             IF (jf.eq.3)
     465            IF (jf.eq.10)
    433466     $          CALL locwrite(cl_writ(jf),fnsolice, imjm,
    434467     $          file_unit_field(jf), ierror)
    435             IF (jf.eq.4)
     468            IF (jf.eq.11)
    436469     $          CALL locwrite(cl_writ(jf),fnsolwat, imjm,
    437470     $          file_unit_field(jf), ierror)
    438             IF (jf.eq.5)
     471            IF (jf.eq.12)
    439472     $          CALL locwrite(cl_writ(jf),fnsicedt, imjm,
    440473     $          file_unit_field(jf), ierror)
    441 c            IF (jf.eq.6)
     474c            IF (jf.eq.13)
    442475c     $          CALL locwrite(cl_writ(jf),ictemp, imjm,
    443476c     $          file_unit_field(jf), ierror)
    444             IF (jf.eq.6)
     477            IF (jf.eq.13)
    445478     $          CALL locwrite(cl_writ(jf),evice, imjm,
    446479     $          file_unit_field(jf), ierror)
    447             IF (jf.eq.7)
     480            IF (jf.eq.14)
    448481     $          CALL locwrite(cl_writ(jf),evwat, imjm,
    449482     $          file_unit_field(jf), ierror)
    450             IF (jf.eq.8)
     483            IF (jf.eq.15)
    451484     $          CALL locwrite(cl_writ(jf),lpre, imjm,
    452485     $          file_unit_field(jf), ierror)
    453             IF (jf.eq.9)
     486            IF (jf.eq.16)
    454487     $          CALL locwrite(cl_writ(jf),spre, imjm,
    455488     $          file_unit_field(jf), ierror)
    456             IF (jf.eq.10)
     489            IF (jf.eq.17)
    457490     $          CALL locwrite(cl_writ(jf),dirunoff, imjm,
    458491     $          file_unit_field(jf), ierror)
    459             IF (jf.eq.11)
     492            IF (jf.eq.18)
    460493     $          CALL locwrite(cl_writ(jf),rivrunoff, imjm,
    461494     $          file_unit_field(jf), ierror)
    462             IF (jf.eq.12)
     495            IF (jf.eq.19)
    463496     $          CALL locwrite(cl_writ(jf),calving, imjm,
    464497     $          file_unit_field(jf), ierror)
    … …  
    475508c$$$     $          CALL locwrite(cl_writ(jf),tauyu, imjm,
    476509c$$$     $          file_unit_field(jf), ierror)
    477             IF (jf.eq.13)
     510            IF (jf.eq.1)
    478511     $          CALL locwrite(cl_writ(jf),tauxx_u, imjm,
    479512     $          file_unit_field(jf),ierror)
    480             IF (jf.eq.14)
     513            IF (jf.eq.2)
    481514     $          CALL locwrite(cl_writ(jf),tauyy_u, imjm,
    482515     $          file_unit_field(jf),ierror)
    483             IF (jf.eq.15)
     516            IF (jf.eq.3)
    484517     $          CALL locwrite(cl_writ(jf),tauzz_u, imjm,
    485518     $          file_unit_field(jf),ierror)
    486             IF (jf.eq.16)
     519            IF (jf.eq.4)
    487520     $          CALL locwrite(cl_writ(jf),tauxx_v, imjm,
    488521     $          file_unit_field(jf),ierror)
    489             IF (jf.eq.17)
     522            IF (jf.eq.5)
    490523     $          CALL locwrite(cl_writ(jf),tauyy_v, imjm,
    491524     $          file_unit_field(jf),ierror)
    492             IF (jf.eq.18)
     525            IF (jf.eq.6)
    493526     $          CALL locwrite(cl_writ(jf),tauzz_v, imjm,
    494527     $          file_unit_field(jf),ierror)
     528c -- LOOP
     529            IF (jf.eq.7)
     530               CALL locwrite(cl_writ(jf),windsp, imjm,
     531     $         file_unit_field(jf),ierror)
     532c -- LOOP
     533
    495534          END DO
    496535C
    … …  
    522561          DO jn=1, jpflda2o1 + jpflda2o2
    523562C           
    524           IF (jn.eq.1) CALL CLIM_Export(cl_writ(jn), kt, fsolice, info)
    525           IF (jn.eq.2) CALL CLIM_Export(cl_writ(jn), kt, fsolwat, info)
    526           IF (jn.eq.3) CALL CLIM_Export(cl_writ(jn), kt, fnsolice, info)
    527           IF (jn.eq.4) CALL CLIM_Export(cl_writ(jn), kt, fnsolwat, info)
    528           IF (jn.eq.5) CALL CLIM_Export(cl_writ(jn), kt, fnsicedt, info)
     563          IF (jn.eq.8) CALL CLIM_Export(cl_writ(jn), kt, fsolice, info)
     564          IF (jn.eq.9) CALL CLIM_Export(cl_writ(jn), kt, fsolwat, info)
     565         IF (jn.eq.10) CALL CLIM_Export(cl_writ(jn), kt, fnsolice, info)
     566         IF (jn.eq.11) CALL CLIM_Export(cl_writ(jn), kt, fnsolwat, info)
     567         IF (jn.eq.12) CALL CLIM_Export(cl_writ(jn), kt, fnsicedt, info)
    529568c          IF (jn.eq.6) CALL CLIM_Export(cl_writ(jn), kt, ictemp, info)
    530           IF (jn.eq.6) CALL CLIM_Export(cl_writ(jn), kt, evice, info)
    531           IF (jn.eq.7) CALL CLIM_Export(cl_writ(jn), kt, evwat, info)
    532           IF (jn.eq.8) CALL CLIM_Export(cl_writ(jn), kt, lpre, info)
    533           IF (jn.eq.9) CALL CLIM_Export(cl_writ(jn), kt, spre, info)
    534           IF (jn.eq.10) CALL CLIM_Export(cl_writ(jn),kt,dirunoff, info)
    535           IF (jn.eq.11) CALL CLIM_Export(cl_writ(jn),kt,rivrunoff,info)
    536           IF (jn.eq.12) CALL CLIM_Export(cl_writ(jn),kt,calving,info)
     569          IF (jn.eq.13) CALL CLIM_Export(cl_writ(jn), kt, evice, info)
     570          IF (jn.eq.14) CALL CLIM_Export(cl_writ(jn), kt, evwat, info)
     571          IF (jn.eq.15) CALL CLIM_Export(cl_writ(jn), kt, lpre, info)
     572          IF (jn.eq.16) CALL CLIM_Export(cl_writ(jn), kt, spre, info)
     573          IF (jn.eq.17) CALL CLIM_Export(cl_writ(jn),kt,dirunoff, info)
     574          IF (jn.eq.18) CALL CLIM_Export(cl_writ(jn),kt,rivrunoff,info)
     575          IF (jn.eq.19) CALL CLIM_Export(cl_writ(jn),kt,calving,info)
    537576c$$$          IF (jn.eq.13) CALL CLIM_Export(cl_writ(jn), kt, tauxu, info)
    538577c$$$          IF (jn.eq.14) CALL CLIM_Export(cl_writ(jn), kt, tauxv, info)
    539578c$$$          IF (jn.eq.15) CALL CLIM_Export(cl_writ(jn), kt, tauyv, info)
    540579c$$$          IF (jn.eq.16) CALL CLIM_Export(cl_writ(jn), kt, tauyu, info)
    541           IF (jn.eq.13) CALL CLIM_Export(cl_writ(jn), kt, tauxx_u, info)
    542           IF (jn.eq.14) CALL CLIM_Export(cl_writ(jn), kt, tauyy_u, info)
    543           IF (jn.eq.15) CALL CLIM_Export(cl_writ(jn), kt, tauzz_u, info)
    544           IF (jn.eq.16) CALL CLIM_Export(cl_writ(jn), kt, tauxx_v, info)
    545           IF (jn.eq.17) CALL CLIM_Export(cl_writ(jn), kt, tauyy_v, info)
    546           IF (jn.eq.18) CALL CLIM_Export(cl_writ(jn), kt, tauzz_v, info)
    547          
     580          IF (jn.eq.1) CALL CLIM_Export(cl_writ(jn), kt, tauxx_u, info)
     581          IF (jn.eq.2) CALL CLIM_Export(cl_writ(jn), kt, tauyy_u, info)
     582          IF (jn.eq.3) CALL CLIM_Export(cl_writ(jn), kt, tauzz_u, info)
     583          IF (jn.eq.4) CALL CLIM_Export(cl_writ(jn), kt, tauxx_v, info)
     584          IF (jn.eq.5) CALL CLIM_Export(cl_writ(jn), kt, tauyy_v, info)
     585          IF (jn.eq.6) CALL CLIM_Export(cl_writ(jn), kt, tauzz_v, info)
     586c -- LOOP
     587          IF (jn.eq.7) CALL CLIM_Export(cl_writ(jn), kt, windsp, info)
     588c -- LOOP         
    548589            IF (info .NE. CLIM_Ok) THEN
    549590                WRITE (nuout,*) 'STEP : Pb giving ',cl_writ(jn), ':',jn

  • LMDZ4/trunk/libf/phylmd/oasis.true

    r524 r590  
    22! $Header$
    33!
    4 C $Id$
    54C****
    65C
    … …  
    2019c     INCLUDE "param.h"
    2120c
    22       INTEGER kastp, kexch, kstep,imjm
     21      INTEGER kastp, kexch, kstep,imjm,klon
    2322      INTEGER iparal(3)
    2423      INTEGER ifcpl, idt, info, imxtag, istep, jf
     24
     25c -- LOOP
     26c
     27#include "dimensions.h"
     28      INTEGER jjmp1
     29      PARAMETER (jjmp1=jjm+1-1/jjm)
     30#include "dimphy.h"
     31c      REAL zu10m(klon), zv10m(klon)
     32       REAL zwindsp(klon)
     33c
     34c -- LOOP
    2535c
    2636#include "param_cou.h"
    … …  
    4353#include "mpiclim.h"
    4454c
    45 #include "oasis.h"      ! contains the name of communication technique. Here
     55#include "oasis.h"     
     56                        ! contains the name of communication technique. Here
    4657                        ! cchan=CLIM only is possible.
    4758c                       ! ctype=MPI2
    … …  
    7182c         must be the same as (1) of the field  definition in namcouple:
    7283c
    73       cl_writ(1)='COSHFICE'
    74       cl_writ(2)='COSHFOCE'
    75       cl_writ(3)='CONSFICE'
    76       cl_writ(4)='CONSFOCE'
    77       cl_writ(5)='CODFLXDT'
     84      cl_writ(8)='COSHFICE'
     85      cl_writ(9)='COSHFOCE'
     86      cl_writ(10)='CONSFICE'
     87      cl_writ(11)='CONSFOCE'
     88      cl_writ(12)='CODFLXDT'
    7889c      cl_writ(6)='COICTEMP'
    79       cl_writ(6)='COTFSICE'
    80       cl_writ(7)='COTFSOCE'
    81       cl_writ(8)='COTOLPSU'
    82       cl_writ(9)='COTOSPSU'
    83       cl_writ(10)='CORUNCOA'
    84       cl_writ(11)='CORIVFLU'
    85       cl_writ(12)='COCALVIN'
     90      cl_writ(13)='COTFSICE'
     91      cl_writ(14)='COTFSOCE'
     92      cl_writ(15)='COTOLPSU'
     93      cl_writ(16)='COTOSPSU'
     94      cl_writ(17)='CORUNCOA'
     95      cl_writ(18)='CORIVFLU'
     96      cl_writ(19)='COCALVIN'
    8697c$$$      cl_writ(13)='COZOTAUX'
    8798c$$$      cl_writ(14)='COZOTAUV'
    8899c$$$      cl_writ(15)='COMETAUY'
    89100c$$$      cl_writ(16)='COMETAUU'
    90       cl_writ(13)='COTAUXXU'
    91       cl_writ(14)='COTAUYYU'
    92       cl_writ(15)='COTAUZZU'
    93       cl_writ(16)='COTAUXXV'
    94       cl_writ(17)='COTAUYYV'
    95       cl_writ(18)='COTAUZZV'
     101      cl_writ(1)='COTAUXXU'
     102      cl_writ(2)='COTAUYYU'
     103      cl_writ(3)='COTAUZZU'
     104      cl_writ(4)='COTAUXXV'
     105      cl_writ(5)='COTAUYYV'
     106      cl_writ(6)='COTAUZZV'
     107c -- LOOP
     108      cl_writ(7)='COWINDSP'
     109c -- LOOP
    96110c
    97111c     Define files name for fields exchanged from atmos to coupler,
    … …  
    116130      cl_f_writ(17)='flxatmos'
    117131      cl_f_writ(18)='flxatmos'
     132c -- LOOP
     133      cl_f_writ(19)='flxatmos'
     134c -- LOOP
    118135
    119136c
    … …  
    311328      SUBROUTINE intocpl(kt, imjm, fsolice, fsolwat, fnsolice, fnsolwat,
    312329     $    fnsicedt, evice, evwat, lpre, spre, dirunoff, rivrunoff,
    313      $    calving, tauxx_u, tauyy_u, tauzz_u, tauxx_v, tauyy_v, tauzz_v
    314      $    , last)
     330     $    calving, tauxx_u, tauyy_u, tauzz_u, tauxx_v, tauyy_v, tauzz_v,
     331c -- LOOP
     332     $    windsp, last)
     333c -- LOOP
    315334c ======================================================================
    316335c S. Valcke (02/99) adapted From L.Z.X Li: this subroutine provides the
    … …  
    320339c ======================================================================
    321340      IMPLICIT NONE
     341c -- LOOP
     342c
     343#include "dimensions.h"
     344      INTEGER jjmp1
     345      PARAMETER (jjmp1=jjm+1-1/jjm)
     346#include "dimphy.h"
     347c      REAL zu10m(klon), zv10m(klon)
     348       REAL zwindsp(klon)
     349c
     350c -- LOOP
     351c
     352
    322353      INTEGER kt, imjm
    323354c
    … …  
    344375      REAL tauzz_u(imjm)
    345376      REAL tauzz_v(imjm)
     377c -- LOOP
     378       REAL windsp(imjm)
     379c -- LOOP
    346380      LOGICAL last
    347381c
    … …  
    423457c         WRITE fields to files         
    424458          DO jf=1, jpflda2o1 + jpflda2o2
    425             IF (jf.eq.1)
     459            IF (jf.eq.8)
    426460     $          CALL locwrite(cl_writ(jf),fsolice, imjm,
    427461     $          file_unit_field(jf), ierror)
    428             IF (jf.eq.2)
     462            IF (jf.eq.9)
    429463     $          CALL locwrite(cl_writ(jf),fsolwat, imjm,
    430464     $          file_unit_field(jf), ierror)
    431             IF (jf.eq.3)
     465            IF (jf.eq.10)
    432466     $          CALL locwrite(cl_writ(jf),fnsolice, imjm,
    433467     $          file_unit_field(jf), ierror)
    434             IF (jf.eq.4)
     468            IF (jf.eq.11)
    435469     $          CALL locwrite(cl_writ(jf),fnsolwat, imjm,
    436470     $          file_unit_field(jf), ierror)
    437             IF (jf.eq.5)
     471            IF (jf.eq.12)
    438472     $          CALL locwrite(cl_writ(jf),fnsicedt, imjm,
    439473     $          file_unit_field(jf), ierror)
    440 c            IF (jf.eq.6)
     474c            IF (jf.eq.13)
    441475c     $          CALL locwrite(cl_writ(jf),ictemp, imjm,
    442476c     $          file_unit_field(jf), ierror)
    443             IF (jf.eq.6)
     477            IF (jf.eq.13)
    444478     $          CALL locwrite(cl_writ(jf),evice, imjm,
    445479     $          file_unit_field(jf), ierror)
    446             IF (jf.eq.7)
     480            IF (jf.eq.14)
    447481     $          CALL locwrite(cl_writ(jf),evwat, imjm,
    448482     $          file_unit_field(jf), ierror)
    449             IF (jf.eq.8)
     483            IF (jf.eq.15)
    450484     $          CALL locwrite(cl_writ(jf),lpre, imjm,
    451485     $          file_unit_field(jf), ierror)
    452             IF (jf.eq.9)
     486            IF (jf.eq.16)
    453487     $          CALL locwrite(cl_writ(jf),spre, imjm,
    454488     $          file_unit_field(jf), ierror)
    455             IF (jf.eq.10)
     489            IF (jf.eq.17)
    456490     $          CALL locwrite(cl_writ(jf),dirunoff, imjm,
    457491     $          file_unit_field(jf), ierror)
    458             IF (jf.eq.11)
     492            IF (jf.eq.18)
    459493     $          CALL locwrite(cl_writ(jf),rivrunoff, imjm,
    460494     $          file_unit_field(jf), ierror)
    461             IF (jf.eq.12)
     495            IF (jf.eq.19)
    462496     $          CALL locwrite(cl_writ(jf),calving, imjm,
    463497     $          file_unit_field(jf), ierror)
    … …  
    474508c$$$     $          CALL locwrite(cl_writ(jf),tauyu, imjm,
    475509c$$$     $          file_unit_field(jf), ierror)
    476             IF (jf.eq.13)
     510            IF (jf.eq.1)
    477511     $          CALL locwrite(cl_writ(jf),tauxx_u, imjm,
    478512     $          file_unit_field(jf),ierror)
    479             IF (jf.eq.14)
     513            IF (jf.eq.2)
    480514     $          CALL locwrite(cl_writ(jf),tauyy_u, imjm,
    481515     $          file_unit_field(jf),ierror)
    482             IF (jf.eq.15)
     516            IF (jf.eq.3)
    483517     $          CALL locwrite(cl_writ(jf),tauzz_u, imjm,
    484518     $          file_unit_field(jf),ierror)
    485             IF (jf.eq.16)
     519            IF (jf.eq.4)
    486520     $          CALL locwrite(cl_writ(jf),tauxx_v, imjm,
    487521     $          file_unit_field(jf),ierror)
    488             IF (jf.eq.17)
     522            IF (jf.eq.5)
    489523     $          CALL locwrite(cl_writ(jf),tauyy_v, imjm,
    490524     $          file_unit_field(jf),ierror)
    491             IF (jf.eq.18)
     525            IF (jf.eq.6)
    492526     $          CALL locwrite(cl_writ(jf),tauzz_v, imjm,
    493527     $          file_unit_field(jf),ierror)
     528c -- LOOP
     529            IF (jf.eq.7)
     530     $          CALL locwrite(cl_writ(jf),windsp, imjm,
     531     $         file_unit_field(jf),ierror)
     532c -- LOOP
     533
    494534          END DO
    495535C
    … …  
    521561          DO jn=1, jpflda2o1 + jpflda2o2
    522562C           
    523           IF (jn.eq.1) CALL CLIM_Export(cl_writ(jn), kt, fsolice, info)
    524           IF (jn.eq.2) CALL CLIM_Export(cl_writ(jn), kt, fsolwat, info)
    525           IF (jn.eq.3) CALL CLIM_Export(cl_writ(jn), kt, fnsolice, info)
    526           IF (jn.eq.4) CALL CLIM_Export(cl_writ(jn), kt, fnsolwat, info)
    527           IF (jn.eq.5) CALL CLIM_Export(cl_writ(jn), kt, fnsicedt, info)
     563          IF (jn.eq.8) CALL CLIM_Export(cl_writ(jn), kt, fsolice, info)
     564          IF (jn.eq.9) CALL CLIM_Export(cl_writ(jn), kt, fsolwat, info)
     565         IF (jn.eq.10) CALL CLIM_Export(cl_writ(jn), kt, fnsolice, info)
     566         IF (jn.eq.11) CALL CLIM_Export(cl_writ(jn), kt, fnsolwat, info)
     567         IF (jn.eq.12) CALL CLIM_Export(cl_writ(jn), kt, fnsicedt, info)
    528568c          IF (jn.eq.6) CALL CLIM_Export(cl_writ(jn), kt, ictemp, info)
    529           IF (jn.eq.6) CALL CLIM_Export(cl_writ(jn), kt, evice, info)
    530           IF (jn.eq.7) CALL CLIM_Export(cl_writ(jn), kt, evwat, info)
    531           IF (jn.eq.8) CALL CLIM_Export(cl_writ(jn), kt, lpre, info)
    532           IF (jn.eq.9) CALL CLIM_Export(cl_writ(jn), kt, spre, info)
    533           IF (jn.eq.10) CALL CLIM_Export(cl_writ(jn),kt,dirunoff, info)
    534           IF (jn.eq.11) CALL CLIM_Export(cl_writ(jn),kt,rivrunoff,info)
    535           IF (jn.eq.12) CALL CLIM_Export(cl_writ(jn),kt,calving,info)
     569          IF (jn.eq.13) CALL CLIM_Export(cl_writ(jn), kt, evice, info)
     570          IF (jn.eq.14) CALL CLIM_Export(cl_writ(jn), kt, evwat, info)
     571          IF (jn.eq.15) CALL CLIM_Export(cl_writ(jn), kt, lpre, info)
     572          IF (jn.eq.16) CALL CLIM_Export(cl_writ(jn), kt, spre, info)
     573          IF (jn.eq.17) CALL CLIM_Export(cl_writ(jn),kt,dirunoff, info)
     574          IF (jn.eq.18) CALL CLIM_Export(cl_writ(jn),kt,rivrunoff,info)
     575          IF (jn.eq.19) CALL CLIM_Export(cl_writ(jn),kt,calving,info)
    536576c$$$          IF (jn.eq.13) CALL CLIM_Export(cl_writ(jn), kt, tauxu, info)
    537577c$$$          IF (jn.eq.14) CALL CLIM_Export(cl_writ(jn), kt, tauxv, info)
    538578c$$$          IF (jn.eq.15) CALL CLIM_Export(cl_writ(jn), kt, tauyv, info)
    539579c$$$          IF (jn.eq.16) CALL CLIM_Export(cl_writ(jn), kt, tauyu, info)
    540           IF (jn.eq.13) CALL CLIM_Export(cl_writ(jn), kt, tauxx_u, info)
    541           IF (jn.eq.14) CALL CLIM_Export(cl_writ(jn), kt, tauyy_u, info)
    542           IF (jn.eq.15) CALL CLIM_Export(cl_writ(jn), kt, tauzz_u, info)
    543           IF (jn.eq.16) CALL CLIM_Export(cl_writ(jn), kt, tauxx_v, info)
    544           IF (jn.eq.17) CALL CLIM_Export(cl_writ(jn), kt, tauyy_v, info)
    545           IF (jn.eq.18) CALL CLIM_Export(cl_writ(jn), kt, tauzz_v, info)
    546          
     580          IF (jn.eq.1) CALL CLIM_Export(cl_writ(jn), kt, tauxx_u, info)
     581          IF (jn.eq.2) CALL CLIM_Export(cl_writ(jn), kt, tauyy_u, info)
     582          IF (jn.eq.3) CALL CLIM_Export(cl_writ(jn), kt, tauzz_u, info)
     583          IF (jn.eq.4) CALL CLIM_Export(cl_writ(jn), kt, tauxx_v, info)
     584          IF (jn.eq.5) CALL CLIM_Export(cl_writ(jn), kt, tauyy_v, info)
     585          IF (jn.eq.6) CALL CLIM_Export(cl_writ(jn), kt, tauzz_v, info)
     586c -- LOOP
     587          IF (jn.eq.7) CALL CLIM_Export(cl_writ(jn), kt, windsp, info)
     588c -- LOOP         
    547589            IF (info .NE. CLIM_Ok) THEN
    548590                WRITE (nuout,*) 'STEP : Pb giving ',cl_writ(jn), ':',jn
Note: See TracChangeset for help on using the changeset viewer.