New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 4351 for branches/2013/dev_ECMWF_waves/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_oce.F90 – NEMO

Ignore:
Timestamp:
2014-01-16T15:27:29+01:00 (10 years ago)
Author:
acc
Message:

Branch 2013/dev_ECMWF_waves ticket #1216. First set of code imports from ECMWF. Untested

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_ECMWF_waves/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_oce.F90

    r4306 r4351  
    5353   LOGICAL , PUBLIC ::   ln_sdw         !: true if 3d stokes drift from wave model 
    5454   ! 
     55   REAL(wp), PUBLIC ::   rn_rrelwind = 1.0       !: relative wind ratio (1.0 means surface currents are subtracted, 0.0 means no current) 
     56   LOGICAL , PUBLIC ::   ln_wavetke  = .FALSE.   !: true if wave parameters are read 
     57   LOGICAL , PUBLIC ::   ln_stcor    = .FALSE.   !: true if Stokes-Coriolis forcing is included  
     58   LOGICAL , PUBLIC ::   ln_tauoc    = .FALSE.   !: true if wave-modified water-side stress is used 
     59 
    5560   LOGICAL , PUBLIC ::   ln_icebergs    !: Icebergs 
    5661   ! 
     
    6974   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   vtau   , vtau_b   !: sea surface j-stress (ocean referential)     [N/m2] 
    7075   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   taum              !: module of sea surface stress (at T-point)    [N/m2]  
     76   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   cdn_wave          !: wave model drag coefficient [N/m2]  
     77   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   wspd_wavepar      !: wave model 10-m neutral wind [m/s] CCC 
     78   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   tauoc_wavepar     !: normalized stress into the ocean from wave model  CCC 
    7179   !! wndm is used onmpute surface gases exchanges in ice-free ocean or leads 
    7280   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   wndm              !: wind speed module at T-point (=|U10m-Uoce|)  [m/s] 
     
    131139      ALLOCATE( rnf  (jpi,jpj) , sbc_tsc  (jpi,jpj,jpts) , qsr_hc  (jpi,jpj,jpk) ,     & 
    132140         &      rnf_b(jpi,jpj) , sbc_tsc_b(jpi,jpj,jpts) , qsr_hc_b(jpi,jpj,jpk) , STAT=ierr(3) ) 
     141         ! Initialize these since that is not done in the code. 
     142      rnf        (:,:) = 0.0_wp 
     143      sbc_tsc  (:,:,:) = 0.0_wp 
     144      qsr_hc   (:,:,:) = 0.0_wp 
     145      rnf_b      (:,:) = 0.0_wp 
     146      sbc_tsc_b(:,:,:) = 0.0_wp 
     147      qsr_hc_b (:,:,:) = 0.0_wp 
     148 
    133149         ! 
    134150      ALLOCATE( tprecip(jpi,jpj) , sprecip(jpi,jpj) , fr_i(jpi,jpj) ,     & 
     
    138154         &      ssu_m  (jpi,jpj) , sst_m(jpi,jpj) ,                       & 
    139155         &      ssv_m  (jpi,jpj) , sss_m  (jpi,jpj), ssh_m(jpi,jpj) , STAT=ierr(4) ) 
     156      ALLOCATE( cdn_wave(jpi,jpj) ) 
     157      ALLOCATE( wspd_wavepar(jpi,jpj) ) 
    140158         ! 
    141159#if defined key_vvl 
Note: See TracChangeset for help on using the changeset viewer.