- Timestamp:
- 2014-01-16T15:27:29+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_ECMWF_waves/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_oce.F90
r4306 r4351 53 53 LOGICAL , PUBLIC :: ln_sdw !: true if 3d stokes drift from wave model 54 54 ! 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 55 60 LOGICAL , PUBLIC :: ln_icebergs !: Icebergs 56 61 ! … … 69 74 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: vtau , vtau_b !: sea surface j-stress (ocean referential) [N/m2] 70 75 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 71 79 !! wndm is used onmpute surface gases exchanges in ice-free ocean or leads 72 80 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: wndm !: wind speed module at T-point (=|U10m-Uoce|) [m/s] … … 131 139 ALLOCATE( rnf (jpi,jpj) , sbc_tsc (jpi,jpj,jpts) , qsr_hc (jpi,jpj,jpk) , & 132 140 & 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 133 149 ! 134 150 ALLOCATE( tprecip(jpi,jpj) , sprecip(jpi,jpj) , fr_i(jpi,jpj) , & … … 138 154 & ssu_m (jpi,jpj) , sst_m(jpi,jpj) , & 139 155 & 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) ) 140 158 ! 141 159 #if defined key_vvl
Note: See TracChangeset
for help on using the changeset viewer.