Changeset 3798
- Timestamp:
- 2013-02-11T15:16:58+01:00 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_MERGE_2012/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfgls.F90
r3680 r3798 44 44 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: mxln !: now mixing length 45 45 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: zwall !: wall function 46 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: avt_k ! not enhanced Kz 47 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: avm_k ! not enhanced Kz 48 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: avmu_k ! not enhanced Kz 49 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: avmv_k ! not enhanced Kz 46 50 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: ustars2 !: Squared surface velocity scale at T-points 47 51 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: ustarb2 !: Squared bottom velocity scale at T-points … … 119 123 !!---------------------------------------------------------------------- 120 124 ALLOCATE( en(jpi,jpj,jpk), mxln(jpi,jpj,jpk), zwall(jpi,jpj,jpk) , & 125 & avt_k (jpi,jpj,jpk) , avm_k (jpi,jpj,jpk), & 126 & avmu_k(jpi,jpj,jpk) , avmv_k(jpi,jpj,jpk), & 121 127 & ustars2(jpi,jpj), ustarb2(jpi,jpj) , STAT= zdf_gls_alloc ) 122 128 ! … … 159 165 160 166 ustars2 = 0._wp ; ustarb2 = 0._wp ; psi = 0._wp ; zwall_psi = 0._wp 167 168 IF( kt /= nit000 ) THEN ! restore before value to compute tke 169 avt (:,:,:) = avt_k (:,:,:) 170 avm (:,:,:) = avm_k (:,:,:) 171 avmu(:,:,:) = avmu_k(:,:,:) 172 avmv(:,:,:) = avmv_k(:,:,:) 173 ENDIF 161 174 162 175 ! Compute surface and bottom friction at T-points … … 882 895 ENDIF 883 896 ! 897 avt_k (:,:,:) = avt (:,:,:) 898 avm_k (:,:,:) = avm (:,:,:) 899 avmu_k(:,:,:) = avmu(:,:,:) 900 avmv_k(:,:,:) = avmv(:,:,:) 901 ! 884 902 CALL wrk_dealloc( jpi,jpj, zdep, zflxs, zhsro ) 885 903 CALL wrk_dealloc( jpi,jpj,jpk, eb, mxlb, shear, eps, zwall_psi, z_elem_a, z_elem_b, z_elem_c, psi ) … … 1244 1262 ! ! ------------------- 1245 1263 IF(lwp) WRITE(numout,*) '---- gls-rst ----' 1246 CALL iom_rstput( kt, nitrst, numrow, 'en' , en )1247 CALL iom_rstput( kt, nitrst, numrow, 'avt' , avt 1248 CALL iom_rstput( kt, nitrst, numrow, 'avm' , avm 1249 CALL iom_rstput( kt, nitrst, numrow, 'avmu' , avmu 1250 CALL iom_rstput( kt, nitrst, numrow, 'avmv' , avmv 1251 CALL iom_rstput( kt, nitrst, numrow, 'mxln' , mxln )1264 CALL iom_rstput( kt, nitrst, numrow, 'en' , en ) 1265 CALL iom_rstput( kt, nitrst, numrow, 'avt' , avt_k ) 1266 CALL iom_rstput( kt, nitrst, numrow, 'avm' , avm_k ) 1267 CALL iom_rstput( kt, nitrst, numrow, 'avmu' , avmu_k ) 1268 CALL iom_rstput( kt, nitrst, numrow, 'avmv' , avmv_k ) 1269 CALL iom_rstput( kt, nitrst, numrow, 'mxln' , mxln ) 1252 1270 ! 1253 1271 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.