Changeset 710 for trunk/NEMO/OPA_SRC
- Timestamp:
- 2007-10-10T17:36:03+02:00 (17 years ago)
- Location:
- trunk/NEMO/OPA_SRC
- Files:
-
- 11 edited
-
DIA/diawri.F90 (modified) (5 diffs)
-
DIA/diawri_dimg.h90 (modified) (4 diffs)
-
DYN/dynzdf_exp.F90 (modified) (1 diff)
-
DYN/dynzdf_imp.F90 (modified) (2 diffs)
-
DYN/dynzdf_imp_jki.F90 (modified) (2 diffs)
-
SBC/sbcblk_clio.F90 (modified) (1 diff)
-
SBC/sbcblk_core.F90 (modified) (1 diff)
-
TRD/trdmod.F90 (modified) (2 diffs)
-
ZDF/zdfkpp.F90 (modified) (1 diff)
-
ZDF/zdftke.F90 (modified) (3 diffs)
-
ZDF/zdftke_jki.F90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DIA/diawri.F90
r708 r710 331 331 #endif 332 332 ! !!! nid_U : 2D 333 CALL histdef( nid_U, "sozotaux", "Wind Stress along i-axis" , "N/m2" , & ! taux333 CALL histdef( nid_U, "sozotaux", "Wind Stress along i-axis" , "N/m2" , & ! utau 334 334 & jpi, jpj, nh_U, 1 , 1, 1 , - 99, 32, clop, zsto, zout ) 335 335 #if defined key_dynspg_rl … … 348 348 #endif 349 349 ! !!! nid_V : 2D 350 CALL histdef( nid_V, "sometauy", "Wind Stress along j-axis" , "N/m2" , & ! tauy350 CALL histdef( nid_V, "sometauy", "Wind Stress along j-axis" , "N/m2" , & ! vtau 351 351 & jpi, jpj, nh_V, 1 , 1, 1 , - 99, 32, clop, zsto, zout ) 352 352 #if defined key_dynspg_rl … … 472 472 CALL histwrite( nid_U, "vozoeivu", it, u_eiv , ndim_U , ndex_U ) ! i-eiv current 473 473 #endif 474 CALL histwrite( nid_U, "sozotaux", it, taux, ndim_hU, ndex_hU ) ! i-wind stress474 CALL histwrite( nid_U, "sozotaux", it, utau , ndim_hU, ndex_hU ) ! i-wind stress 475 475 #if defined key_dynspg_rl 476 476 CALL lbc_lnk( spgu, 'U', -1. ) … … 483 483 CALL histwrite( nid_V, "vomeeivv", it, v_eiv , ndim_V , ndex_V ) ! j-eiv current 484 484 #endif 485 CALL histwrite( nid_V, "sometauy", it, tauy, ndim_hV, ndex_hV ) ! j-wind stress485 CALL histwrite( nid_V, "sometauy", it, vtau , ndim_hV, ndex_hV ) ! j-wind stress 486 486 #if defined key_dynspg_rl 487 487 CALL lbc_lnk( spgv, 'V', -1. ) … … 651 651 CALL histwrite( id_i, "soshfldo", 1, qsr , jpi*jpj , idex ) ! total heat flux 652 652 CALL histwrite( id_i, "soicecov", 1, freeze, jpi*jpj , idex ) ! ice cover 653 CALL histwrite( id_i, "sozotaux", 1, taux, jpi*jpj , idex ) ! i-wind stress654 CALL histwrite( id_i, "sometauy", 1, tauy, jpi*jpj , idex ) ! j-wind stress653 CALL histwrite( id_i, "sozotaux", 1, utau , jpi*jpj , idex ) ! i-wind stress 654 CALL histwrite( id_i, "sometauy", 1, vtau , jpi*jpj , idex ) ! j-wind stress 655 655 656 656 ! 3. Close the file -
trunk/NEMO/OPA_SRC/DIA/diawri_dimg.h90
r699 r710 41 41 !! To be tested with a lot of procs !!!! 42 42 !! 43 !! level 1: taux(:,:) * umask(:,:,1) zonal stress in N.m-244 !! level 2: tauy(:,:) * vmask(:,:,1) meridional stress in N. m-243 !! level 1: utau(:,:) * umask(:,:,1) zonal stress in N.m-2 44 !! level 2: vtau(:,:) * vmask(:,:,1) meridional stress in N. m-2 45 45 !! level 3: qt (:,:) total heat flux (W/m2) 46 46 !! level 4: emp (:,:) E-P flux (mm/day) … … 167 167 sm(:,:,:)=sm(:,:,:) + sn (:,:,:) 168 168 ! 169 fsel(:,:,1 ) = fsel(:,:,1 ) + taux(:,:) * umask(:,:,1)170 fsel(:,:,2 ) = fsel(:,:,2 ) + tauy(:,:) * vmask(:,:,1)169 fsel(:,:,1 ) = fsel(:,:,1 ) + utau(:,:) * umask(:,:,1) 170 fsel(:,:,2 ) = fsel(:,:,2 ) + vtau(:,:) * vmask(:,:,1) 171 171 fsel(:,:,3 ) = fsel(:,:,3 ) + qt (:,:) 172 172 fsel(:,:,4 ) = fsel(:,:,4 ) + emp (:,:) … … 231 231 fsel(:,:,20)= spgv(:,:) 232 232 #endif 233 ! mask mean field with tmask except taux tauy(1,2)233 ! mask mean field with tmask except utau vtau (1,2) 234 234 DO jk=3,inbsel 235 235 fsel(:,:,jk)=fsel(:,:,jk)*tmask(:,:,1) … … 256 256 fsel(:,:,:) = 0._wp 257 257 ! 258 fsel(:,:,1 ) = taux(:,:) * umask(:,:,1)259 fsel(:,:,2 ) = tauy(:,:) * vmask(:,:,1)258 fsel(:,:,1 ) = utau(:,:) * umask(:,:,1) 259 fsel(:,:,2 ) = vtau(:,:) * vmask(:,:,1) 260 260 fsel(:,:,3 ) = qt (:,:) * tmask(:,:,1) 261 261 fsel(:,:,4 ) = emp (:,:) * tmask(:,:,1) -
trunk/NEMO/OPA_SRC/DYN/dynzdf_exp.F90
r708 r710 81 81 ! Surface boundary condition 82 82 DO ji = 2, jpim1 83 zwy(ji,1) = taux(ji,jj) * zrau0r84 zww(ji,1) = tauy(ji,jj) * zrau0r83 zwy(ji,1) = utau(ji,jj) * zrau0r 84 zww(ji,1) = vtau(ji,jj) * zrau0r 85 85 END DO 86 86 -
trunk/NEMO/OPA_SRC/DYN/dynzdf_imp.F90
r708 r710 141 141 !!! change les resultats (derniers digit, pas significativement + rapide 1* de moins) 142 142 !!! ua(ji,jj,1) = ub(ji,jj,1) & 143 !!! + p2dt * ( ua(ji,jj,1) + taux(ji,jj) / ( fse3u(ji,jj,1)*rau0 ) )143 !!! + p2dt * ( ua(ji,jj,1) + utau(ji,jj) / ( fse3u(ji,jj,1)*rau0 ) ) 144 144 z2dtf = p2dt / ( fse3u(ji,jj,1)*rau0 ) 145 145 ua(ji,jj,1) = ub(ji,jj,1) & 146 + p2dt * ua(ji,jj,1) + z2dtf * taux(ji,jj)146 + p2dt * ua(ji,jj,1) + z2dtf * utau(ji,jj) 147 147 END DO 148 148 END DO … … 236 236 !!! change les resultats (derniers digit, pas significativement + rapide 1* de moins) 237 237 !!! va(ji,jj,1) = vb(ji,jj,1) & 238 !!! + p2dt * ( va(ji,jj,1) + tauy(ji,jj) / ( fse3v(ji,jj,1)*rau0 ) )238 !!! + p2dt * ( va(ji,jj,1) + vtau(ji,jj) / ( fse3v(ji,jj,1)*rau0 ) ) 239 239 z2dtf = p2dt / ( fse3v(ji,jj,1)*rau0 ) 240 240 va(ji,jj,1) = vb(ji,jj,1) & 241 + p2dt * va(ji,jj,1) + z2dtf * tauy(ji,jj)241 + p2dt * va(ji,jj,1) + z2dtf * vtau(ji,jj) 242 242 END DO 243 243 END DO -
trunk/NEMO/OPA_SRC/DYN/dynzdf_imp_jki.F90
r708 r710 105 105 zwi(ji,1) = 0. 106 106 zwd(ji,1) = 1. - zws(ji,1) 107 zwy(ji,1) = zwy(ji,1) + z2dtf * taux(ji,jj)107 zwy(ji,1) = zwy(ji,1) + z2dtf * utau(ji,jj) 108 108 END DO 109 109 … … 192 192 zwi(ji,1) = 0.e0 193 193 zwd(ji,1) = 1. - zws(ji,1) 194 zwy(ji,1) = zwy(ji,1) + z2dtf * tauy(ji,jj)194 zwy(ji,1) = zwy(ji,1) + z2dtf * vtau(ji,jj) 195 195 END DO 196 196 -
trunk/NEMO/OPA_SRC/SBC/sbcblk_clio.F90
r702 r710 94 94 !! 95 95 !! ** Method : READ each fluxes in NetCDF files 96 !! The i-component of the stress taux(N/m2)97 !! The j-component of the stress tauy(N/m2)96 !! The i-component of the stress utau (N/m2) 97 !! The j-component of the stress vtau (N/m2) 98 98 !! the net downward heat flux qtot (watt/m2) 99 99 !! the net downward radiative flux qsr (watt/m2) -
trunk/NEMO/OPA_SRC/SBC/sbcblk_core.F90
r702 r710 76 76 !! 77 77 !! ** Method : READ each fluxes in NetCDF files 78 !! The i-component of the stress taux(N/m2)79 !! The j-component of the stress tauy(N/m2)78 !! The i-component of the stress utau (N/m2) 79 !! The j-component of the stress vtau (N/m2) 80 80 !! the net downward heat flux qtot (watt/m2) 81 81 !! the net downward radiative flux qsr (watt/m2) -
trunk/NEMO/OPA_SRC/TRD/trdmod.F90
r708 r710 124 124 DO ji = fs_2, fs_jpim1 ! vector opt. 125 125 ! save the surface forcing momentum fluxes 126 ztswu(ji,jj) = taux(ji,jj) / ( fse3u(ji,jj,1)*rau0 )127 ztswv(ji,jj) = tauy(ji,jj) / ( fse3v(ji,jj,1)*rau0 )126 ztswu(ji,jj) = utau(ji,jj) / ( fse3u(ji,jj,1)*rau0 ) 127 ztswv(ji,jj) = vtau(ji,jj) / ( fse3v(ji,jj,1)*rau0 ) 128 128 ! save bottom friction momentum fluxes 129 129 ikbu = MIN( mbathy(ji+1,jj ), mbathy(ji,jj) ) … … 175 175 DO ji = fs_2, fs_jpim1 ! vector opt. 176 176 ! save the surface forcing momentum fluxes 177 ztswu(ji,jj) = taux(ji,jj) / ( fse3u(ji,jj,1)*rau0 )178 ztswv(ji,jj) = tauy(ji,jj) / ( fse3v(ji,jj,1)*rau0 )177 ztswu(ji,jj) = utau(ji,jj) / ( fse3u(ji,jj,1)*rau0 ) 178 ztswv(ji,jj) = vtau(ji,jj) / ( fse3v(ji,jj,1)*rau0 ) 179 179 ! save bottom friction momentum fluxes 180 180 ikbu = MIN( mbathy(ji+1,jj ), mbathy(ji,jj) ) -
trunk/NEMO/OPA_SRC/ZDF/zdfkpp.F90
r708 r710 475 475 zrhos = rhop(ji,jj,1) + zflageos * rau0 * ( 1. - tmask(ji,jj,1) ) 476 476 ! Friction velocity (zustar), at T-point : LMD94 eq. 2 477 ztx = 0.5 * ( taux(ji,jj) + taux(ji - 1, jj ) )478 zty = 0.5 * ( tauy(ji,jj) + tauy(ji , jj - 1) )477 ztx = 0.5 * ( utau(ji,jj) + utau(ji - 1, jj ) ) 478 zty = 0.5 * ( vtau(ji,jj) + vtau(ji , jj - 1) ) 479 479 ztau = SQRT( ztx * ztx + zty * zty ) 480 480 zustar(ji,jj) = SQRT( ztau / ( zrhos + epsln ) ) -
trunk/NEMO/OPA_SRC/ZDF/zdftke.F90
r708 r710 101 101 !! - ediss / emxl en**(2/3) ! dissipation 102 102 !! with the boundary conditions: 103 !! surface: en = max( emin0,ebb sqrt( taux^2 + tauy^2) )103 !! surface: en = max( emin0,ebb sqrt(utau^2 + vtau^2) ) 104 104 !! bottom : en = emin 105 105 !! -1- The dissipation and mixing turbulent lengh scales are computed … … 301 301 ! 2. Surface boundary condition on tke and its eddy viscosity (zmxlm) 302 302 ! ------------------------------------------------- 303 ! en(1) = ebb sqrt( taux^2+tauy^2) / rau0 (min value emin0)303 ! en(1) = ebb sqrt(utau^2+vtau^2) / rau0 (min value emin0) 304 304 ! zmxlm(1) = avmb(1) and zmxlm(jpk) = 0. 305 305 !CDIR NOVERRCHK … … 307 307 !CDIR NOVERRCHK 308 308 DO ji = fs_2, fs_jpim1 ! vector opt. 309 ztx2 = taux(ji-1,jj ) + taux(ji,jj)310 zty2 = tauy(ji ,jj-1) + tauy(ji,jj)309 ztx2 = utau(ji-1,jj ) + utau(ji,jj) 310 zty2 = vtau(ji ,jj-1) + vtau(ji,jj) 311 311 zesurf = zbbrau * SQRT( ztx2 * ztx2 + zty2 * zty2 ) 312 312 en (ji,jj,1) = MAX( zesurf, emin0 ) * tmask(ji,jj,1) -
trunk/NEMO/OPA_SRC/ZDF/zdftke_jki.F90
r708 r710 57 57 !! - ediss / emxl en**(2/3) ! dissipation 58 58 !! with the boundary conditions: 59 !! surface: en = max( emin0,ebb sqrt( taux^2 + tauy^2) )59 !! surface: en = max( emin0,ebb sqrt(utau^2 + vtau^2) ) 60 60 !! bottom : en = emin 61 61 !! -1- The dissipation and mixing turbulent lengh scales are computed … … 253 253 ! 2. Surface boundary condition on tke and its eddy viscosity (zmxlm) 254 254 ! ------------------------------------------------- 255 ! en(1) = ebb sqrt( taux^2+tauy^2) / rau0 (min value emin0)255 ! en(1) = ebb sqrt(utau^2+vtau^2) / rau0 (min value emin0) 256 256 ! zmxlm(1) = avmb(1) and zmxlm(jpk) = 0. 257 257 !CDIR NOVERRCHK 258 258 DO ji = 2, jpim1 259 ztx2 = taux(ji-1,jj ) + taux(ji,jj)260 zty2 = tauy(ji ,jj-1) + tauy(ji,jj)259 ztx2 = utau(ji-1,jj ) + utau(ji,jj) 260 zty2 = vtau(ji ,jj-1) + vtau(ji,jj) 261 261 zesurf = zbbrau * SQRT( ztx2 * ztx2 + zty2 * zty2 ) 262 262 en (ji,jj,1) = MAX( zesurf, emin0 ) * tmask(ji,jj,1)
Note: See TracChangeset
for help on using the changeset viewer.
