Changeset 3804
- Timestamp:
- 2013-02-12T14:16:51+01:00 (11 years ago)
- Location:
- branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPA_SRC/ZDF
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfevd.F90
r3294 r3804 19 19 USE zdf_oce ! ocean vertical physics variables 20 20 USE zdfkpp ! KPP vertical mixing 21 USE zdfgls ! GLS vertical mixing 21 22 USE in_out_manager ! I/O manager 22 23 USE iom ! for iom_put … … 67 68 IF(lwp) WRITE(numout,*) '~~~~~~~ ' 68 69 IF(lwp) WRITE(numout,*) 70 ! 71 IF(lwp .AND. lk_zdfgls ) CALL ctl_warn(' No need zdf_evd with GLS closures ') 72 ! 69 73 ENDIF 70 74 -
branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfgls.F90
r3558 r3804 45 45 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: mxln !: now mixing length 46 46 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: zwall !: wall function 47 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: avt_k ! not enhanced Kz 48 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: avm_k ! not enhanced Kz 49 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: avmu_k ! not enhanced Kz 50 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: avmv_k ! not enhanced Kz 47 51 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: ustars2 !: Squared surface velocity scale at T-points 48 52 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: ustarb2 !: Squared bottom velocity scale at T-points … … 120 124 !!---------------------------------------------------------------------- 121 125 ALLOCATE( en(jpi,jpj,jpk), mxln(jpi,jpj,jpk), zwall(jpi,jpj,jpk) , & 126 & avt_k (jpi,jpj,jpk) , avm_k (jpi,jpj,jpk), & 127 & avmu_k(jpi,jpj,jpk) , avmv_k(jpi,jpj,jpk), & 122 128 & ustars2(jpi,jpj), ustarb2(jpi,jpj) , STAT= zdf_gls_alloc ) 123 129 ! … … 160 166 161 167 ustars2 = 0._wp ; ustarb2 = 0._wp ; psi = 0._wp ; zwall_psi = 0._wp 168 169 IF( kt /= nit000 ) THEN ! restore before value to compute tke 170 avt (:,:,:) = avt_k (:,:,:) 171 avm (:,:,:) = avm_k (:,:,:) 172 avmu(:,:,:) = avmu_k(:,:,:) 173 avmv(:,:,:) = avmv_k(:,:,:) 174 ENDIF 162 175 163 176 ! Compute surface and bottom friction at T-points … … 883 896 ENDIF 884 897 ! 898 avt_k (:,:,:) = avt (:,:,:) 899 avm_k (:,:,:) = avm (:,:,:) 900 avmu_k(:,:,:) = avmu(:,:,:) 901 avmv_k(:,:,:) = avmv(:,:,:) 902 ! 885 903 CALL wrk_dealloc( jpi,jpj, zdep, zflxs, zhsro ) 886 904 CALL wrk_dealloc( jpi,jpj,jpk, eb, mxlb, shear, eps, zwall_psi, z_elem_a, z_elem_b, z_elem_c, psi ) … … 1246 1264 IF(lwp) WRITE(numout,*) '---- gls-rst ----' 1247 1265 CALL iom_rstput( kt, nitrst, numrow, 'en' , en ) 1248 CALL iom_rstput( kt, nitrst, numrow, 'avt' , avt 1249 CALL iom_rstput( kt, nitrst, numrow, 'avm' , avm 1250 CALL iom_rstput( kt, nitrst, numrow, 'avmu' , avmu 1251 CALL iom_rstput( kt, nitrst, numrow, 'avmv' , avmv 1266 CALL iom_rstput( kt, nitrst, numrow, 'avt' , avt_k ) 1267 CALL iom_rstput( kt, nitrst, numrow, 'avm' , avm_k ) 1268 CALL iom_rstput( kt, nitrst, numrow, 'avmu' , avmu_k ) 1269 CALL iom_rstput( kt, nitrst, numrow, 'avmv' , avmv_k ) 1252 1270 CALL iom_rstput( kt, nitrst, numrow, 'mxln' , mxln ) 1253 1271 !
Note: See TracChangeset
for help on using the changeset viewer.