Changeset 1013


Ignore:
Timestamp:
2008-05-30T12:15:06+02:00 (13 years ago)
Author:
cetlod
Message:

adding wind speed module variable, see ticket 172

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OFF_SRC/dtadyn.F90

    r977 r1013  
    3838      jptaux = 1 , & ! indice in flux for taux 
    3939      jptauy = 2 , & ! indice in flux for tauy 
    40       jpemp = 3  , & ! indice in flux for E-P 
    41       jpice = 4  , & ! indice in flux for ice concentration 
    42       jpqsr = 5      ! indice in flux for shortwave heat flux  
     40      jpwind = 3 , & ! indice in flux for wind speed 
     41      jpemp = 4  , & ! indice in flux for E-P 
     42      jpice = 5  , & ! indice in flux for ice concentration 
     43      jpqsr = 6      ! indice in flux for shortwave heat flux  
    4344 
    4445   LOGICAL , PUBLIC :: & 
     
    561562       ! keep needed fluxes 
    562563       ! 
     564                    wndm(:,:) = flx(:,:,jpwind) 
    563565                    freeze(:,:) = flx(:,:,jpice) 
    564566                    emp(:,:) = flx(:,:,jpemp) 
     
    625627       ! keep needed fluxes 
    626628       ! 
     629                  wndm(:,:) = flx(:,:,jpwind) 
    627630                  freeze(:,:) = flx(:,:,jpice) 
    628631                  emp(:,:)    = flx(:,:,jpemp) 
     
    680683 
    681684      REAL(wp), DIMENSION(jpi,jpj) :: & 
    682          zemp, zqsr, zmld, zice  
     685         zemp, zqsr, zmld, zice, zwspd 
    683686#if defined key_trcbbl_dif   ||   defined key_trcbbl_adv 
    684687      REAL(wp), DIMENSION(jpi,jpj) :: & 
     
    756759      CALL iom_get ( numfl_t, jpdom_data, 'soshfldo', zqsr (:,:  ), jkenr ) 
    757760      CALL iom_get ( numfl_t, jpdom_data, 'soicecov', zice (:,:  ), jkenr ) 
     761      CALL iom_get ( numfl_t, jpdom_data, 'sowindsp', zwspd(:,:  ), jkenr ) 
    758762 
    759763      ! file grid-U 
     
    887891      ! flux : 
    888892      ! 
     893      flxdta(:,:,jpwind,2) = zwspd(:,:) * tmask(:,:,1) 
    889894      flxdta(:,:,jpice,2)  = MIN( 1., zice(:,:) ) * tmask(:,:,1) 
    890895      flxdta(:,:,jpemp,2)  = zemp(:,:) * tmask(:,:,1) 
Note: See TracChangeset for help on using the changeset viewer.