Changeset 888 for trunk/NEMO/C1D_SRC
- Timestamp:
- 2008-04-11T19:05:03+02:00 (16 years ago)
- Location:
- trunk/NEMO/C1D_SRC
- Files:
-
- 3 edited
-
diawri1d.F90 (modified) (14 diffs)
-
icestp1d.F90 (modified) (7 diffs)
-
step1d.F90 (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/C1D_SRC/diawri1d.F90
r833 r888 13 13 USE dom_oce ! ocean space and time domain 14 14 USE zdf_oce ! ocean vertical physics 15 USE sbc_oce ! surface boundary condition: ocean 16 USE sbc_ice ! surface boundary condition: ice 15 17 USE zdftke ! TKE vertical mixing 16 18 USE zdfkpp ! KPP vertical mixing … … 19 21 USE phycst ! physical constants 20 22 USE ocfzpt ! ??? 21 USE ocesbc ! surface thermohaline fluxes22 USE taumod ! surface stress23 USE flxrnf ! ???24 23 USE zdfmxl ! mixed layer 25 24 USE daymod ! calendar … … 49 48 !!---------------------------------------------------------------------- 50 49 !! OPA 9.0 , LOCEAN-IPSL (2005) 51 !! $ Header$50 !! $Id$ 52 51 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 53 52 !!---------------------------------------------------------------------- … … 194 193 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 195 194 196 #if ! defined key_dynspg_rl && defined key_lim3 197 ! sowaflup = sowaflep + sorunoff + sowafldp + a term associated to198 ! internal damping to Levitus that can be diagnosed from others199 ! sowaflcd = sowaflep + sorunoff + sowafldp + iowaflup200 CALL histdef( nid_T, "iowaflup", "Ice=>ocean net freshwater" , "kg/m2/s", & ! fsalt201 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )202 CALL histdef( nid_T, "sowaflep", "atmos=>ocean net freshwater" , "kg/m2/s", & ! fmass203 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )204 #endif195 !!$#if ! defined key_dynspg_rl && ( defined key_lim2 || defined key_lim2 ) 196 !!$ ! sowaflup = sowaflep + sorunoff + sowafldp + a term associated to 197 !!$ ! internal damping to Levitus that can be diagnosed from others 198 !!$ ! sowaflcd = sowaflep + sorunoff + sowafldp + iowaflup 199 !!$ CALL histdef( nid_T, "iowaflup", "Ice=>ocean net freshwater" , "kg/m2/s", & ! fsalt 200 !!$ & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 201 !!$ CALL histdef( nid_T, "sowaflep", "atmos=>ocean net freshwater" , "kg/m2/s", & ! fmass 202 !!$ & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 203 !!$#endif 205 204 CALL histdef( nid_T, "sowaflup", "Net Upward Water Flux" , "Kg/m2/s", & ! emp 206 205 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) … … 211 210 CALL histdef( nid_T, "sosalflx", "Surface Salt Flux" , "Kg/m2/s", & ! emps * sn 212 211 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 213 CALL histdef( nid_T, "sohefldo", "Net Downward Heat Flux" , "W/m2" , & ! q t212 CALL histdef( nid_T, "sohefldo", "Net Downward Heat Flux" , "W/m2" , & ! qsr + qns 214 213 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 215 214 CALL histdef( nid_T, "soshfldo", "Shortwave Radiation" , "W/m2" , & ! qsr … … 238 237 #endif 239 238 240 #if ( defined key_coupled && ! defined key_lim3)239 #if ( defined key_coupled && ! ( defined key_lim3 || defined key_lim2 ) ) 241 240 CALL histdef( nid_T, "sohefldp", "Surface Heat Flux: Damping" , "W/m2" , & ! qrp 242 241 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) … … 260 259 #endif 261 260 262 #if defined key_lim3&& defined key_coupled261 #if ( defined key_lim3 || defined key_lim2 ) && defined key_coupled 263 262 CALL histdef( nid_T,"soicetem" , "Ice Surface Temperature" , "K" , & ! tn_ice 264 263 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) … … 275 274 #endif 276 275 ! !!! nid_U : 2D 277 CALL histdef( nid_T, "sozotaux", "Wind Stress along i-axis" , "N/m2" , & ! taux276 CALL histdef( nid_T, "sozotaux", "Wind Stress along i-axis" , "N/m2" , & ! utau 278 277 & jpi, jpj, nh_T, 1 , 1, 1 , - 99, 32, clop, zsto, zout ) 279 278 … … 286 285 #endif 287 286 ! !!! nid_V : 2D 288 CALL histdef( nid_T, "sometauy", "Wind Stress along j-axis" , "N/m2" , & ! tauy287 CALL histdef( nid_T, "sometauy", "Wind Stress along j-axis" , "N/m2" , & ! vtau 289 288 & jpi, jpj, nh_T, 1 , 1, 1 , - 99, 32, clop, zsto, zout ) 290 289 #if defined key_zdftke … … 365 364 CALL histwrite( nid_T, "sosstsst", it, tn(:,:,1) , ndim_hT, ndex_hT ) ! sea surface temperature 366 365 CALL histwrite( nid_T, "sosaline", it, sn(:,:,1) , ndim_hT, ndex_hT ) ! sea surface salinity 367 #if ! defined key_dynspg_rl && defined key_lim3366 #if ! defined key_dynspg_rl && ( defined key_lim3 || defined key_lim2 ) 368 367 CALL histwrite( nid_T, "iowaflup", it, fsalt(:,:) , ndim_hT, ndex_hT ) ! ice=>ocean water flux 369 368 CALL histwrite( nid_T, "sowaflep", it, fmass(:,:) , ndim_hT, ndex_hT ) ! atmos=>ocean water flux … … 374 373 zw2d(:,:) = emps(:,:) * sn(:,:,1) * tmask(:,:,1) 375 374 CALL histwrite( nid_T, "sosalflx", it, zw2d , ndim_hT, ndex_hT ) ! c/d salt flux 376 CALL histwrite( nid_T, "sohefldo", it, q t, ndim_hT, ndex_hT ) ! total heat flux375 CALL histwrite( nid_T, "sohefldo", it, qsr + qns , ndim_hT, ndex_hT ) ! total heat flux 377 376 CALL histwrite( nid_T, "soshfldo", it, qsr , ndim_hT, ndex_hT ) ! solar heat flux 378 377 CALL histwrite( nid_T, "somxl010", it, hmlp , ndim_hT, ndex_hT ) ! mixed layer depth … … 397 396 CALL histwrite( nid_T, "sosafldp", it, zw2d , ndim_hT, ndex_hT ) ! salt flux damping 398 397 #endif 399 #if ( defined key_coupled && ! defined key_lim3)398 #if ( defined key_coupled && ! ( defined key_lim3 || defined key_lim2 ) ) 400 399 CALL histwrite( nid_T, "sohefldp", it, qrp , ndim_hT, ndex_hT ) ! heat flux damping 401 400 CALL histwrite( nid_T, "sowafldp", it, erp , ndim_hT, ndex_hT ) ! freshwater flux damping … … 412 411 CALL histwrite( nid_T, "sohtc300", it, htc3 , ndim_hT, ndex_hT ) ! first 300m heaat content 413 412 #endif 414 #if defined key_lim3&& defined key_coupled413 #if ( defined key_lim3 || defined key_lim2 ) && defined key_coupled 415 414 CALL histwrite( nid_T, "soicetem", it, tn_ice , ndim_hT, ndex_hT ) ! surf. ice temperature 416 415 CALL histwrite( nid_T, "soicealb", it, alb_ice , ndim_hT, ndex_hT ) ! ice albedo … … 418 417 419 418 CALL histwrite( nid_T, "vozocrtx", it, un , ndim_T , ndex_T ) ! i-current 420 CALL histwrite( nid_T, "sozotaux", it, taux, ndim_hT, ndex_hT ) ! i-wind stress419 CALL histwrite( nid_T, "sozotaux", it, utau , ndim_hT, ndex_hT ) ! i-wind stress 421 420 CALL histwrite( nid_T, "vomecrty", it, vn , ndim_T , ndex_T ) ! j-current 422 CALL histwrite( nid_T, "sometauy", it, tauy, ndim_hT, ndex_hT ) ! j-wind stress421 CALL histwrite( nid_T, "sometauy", it, vtau , ndim_hT, ndex_hT ) ! j-wind stress 423 422 #if defined key_zdftke 424 423 CALL histwrite( nid_T, "votlsdis", it, e_dis , ndim_T , ndex_T ) ! Diss. Turb. lenght scale -
trunk/NEMO/C1D_SRC/icestp1d.F90
r833 r888 6 6 !! History : 9.0 ! 04-10 (C. Ethe) from icestp, 1D configuration 7 7 !!---------------------------------------------------------------------- 8 #if defined key_cfg_1d && defined key_lim38 #if defined key_cfg_1d && ( defined key_lim3 || defined key_lim2 ) 9 9 !!---------------------------------------------------------------------- 10 10 !! 'key_cfg_1d' .AND. 1D Configuration 11 !! 'key_lim 3' Limsea-ice model11 !! 'key_lim2' OR 'key_lim3' : LIM 2.0 or 3.0 sea-ice model 12 12 !!---------------------------------------------------------------------- 13 13 !!---------------------------------------------------------------------- … … 18 18 USE in_out_manager ! I/O manager 19 19 USE ice_oce ! ice variables 20 USE flx_oce ! forcings variables 21 USE dom_ice ! LIM sea-ice domain 22 USE cpl_oce ! coupled ocean-atmosphere variables 23 USE blk_oce ! bulk variables 20 USE dom_ice_2 ! LIM sea-ice domain 21 USE sbc_oce ! surface boundary condition: ocean 22 USE sbc_ice ! surface boundary condition: ice 24 23 USE daymod ! calendar 25 24 USE phycst ! Define parameters for the routines 26 USE taumod ! surface stress forcing 27 USE ice ! ice variables 25 USE ice_2 ! ice variables 28 26 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 29 USE limthd 30 USE limflx 31 USE limwri 32 USE limrst 33 34 USE ocesbc ! thermohaline fluxes 35 USE flxmod ! thermohaline forcing 36 USE flxrnf ! runoffs forcing 27 USE limthd_2 28 USE limwri_2 29 USE limrst_2 30 37 31 USE tradmp ! damping salinity trend 38 32 USE dtatem ! ocean temperature data … … 52 46 !!---------------------------------------------------------------------- 53 47 !! LIM 2.0 , UCL-LOCEAN-IPSL (2006) 54 !! $ Header$48 !! $Id$ 55 49 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 56 50 !!---------------------------------------------------------------------- … … 109 103 u_io (:,:) = u_io (:,:) / FLOAT( nfice ) 110 104 v_io (:,:) = v_io (:,:) / FLOAT( nfice ) 111 gtaux (:,:) = taux(:,:)112 gtauy (:,:) = tauy(:,:)105 gtaux (:,:) = utau (:,:) 106 gtauy (:,:) = vtau (:,:) 113 107 114 108 zsss_io (:,:) = SQRT( sss_io(:,:) ) … … 220 214 IF( kt == nit000 ) THEN 221 215 qsr (:,:) = 0.e0 222 q t(:,:) = 0.e0216 qns (:,:) = 0.e0 223 217 qrp (:,:) = 0.e0 224 218 emp (:,:) = 0.e0 … … 238 232 ! ----------------- 239 233 240 q t (:,:) = fnsolar(:,:) + fsolar(:,:) ! non solar heat flux + solarflux234 qns (:,:) = fnsolar(:,:) ! non solar heat flux 241 235 qsr (:,:) = fsolar(:,:) ! solar flux 242 236 … … 261 255 DO ji = 1, fs_jpim1 ! vertor opt. 262 256 ztxy = freezn(ji,jj) ! ice/ocean indicator at T-points 263 taux(ji,jj) = (1.-ztxy) * taux(ji,jj) + ztxy * ftaux(ji,jj) ! stress at the ocean surface264 tauy(ji,jj) = (1.-ztxy) * tauy(ji,jj) + ztxy * ftauy(ji,jj)265 END DO 266 END DO 267 268 ! boundary condition on the stress ( taux,tauy)269 CALL lbc_lnk( taux, 'U', -1. )270 CALL lbc_lnk( tauy, 'V', -1. )257 utau(ji,jj) = (1.-ztxy) * utau(ji,jj) + ztxy * ftaux(ji,jj) ! stress at the ocean surface 258 vtau(ji,jj) = (1.-ztxy) * vtau(ji,jj) + ztxy * ftauy(ji,jj) 259 END DO 260 END DO 261 262 ! boundary condition on the stress (utau,vtau) 263 CALL lbc_lnk( utau, 'U', -1. ) 264 CALL lbc_lnk( vtau, 'V', -1. ) 271 265 272 266 ! Re-initialization of fluxes -
trunk/NEMO/C1D_SRC/step1d.F90
r719 r888 15 15 USE dom_oce ! ocean space and time domain variables 16 16 USE zdf_oce ! ocean vertical physics variables 17 USE sbc_oce ! surface boundary condition: ocean 17 18 USE ldftra_oce 18 19 USE ldfdyn_oce … … 24 25 USE dtatem ! ocean temperature data (dta_tem routine) 25 26 USE dtasal ! ocean salinity data (dta_sal routine) 26 USE dtasst ! ocean sea surface temerature (dta_sst routine)27 USE taumod ! surface stress (tau routine)28 USE flxmod ! thermohaline fluxes (flx routine)29 USE ocesbc ! thermohaline fluxes (oce_sbc routine)30 USE flxrnf ! runoffs (flx_rnf routine)31 USE flxfwb ! freshwater budget correction (flx_fwb routine)32 27 USE ocfzpt ! surface ocean freezing point (oc_fz_pt routine) 33 28 … … 75 70 !!---------------------------------------------------------------------- 76 71 !! OPA 9.0 , LOCEAN-IPSL (2005) 77 !! $ Header$72 !! $Id$ 78 73 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 79 74 !!---------------------------------------------------------------------- … … 157 152 CALL prt_ctl(tab2d_1=emp , clinfo1=' emp - : ', mask1=tmask, ovlap=1) 158 153 CALL prt_ctl(tab2d_1=emps , clinfo1=' emps - : ', mask1=tmask, ovlap=1) 159 CALL prt_ctl(tab2d_1=q t , clinfo1=' qt- : ', mask1=tmask, ovlap=1)154 CALL prt_ctl(tab2d_1=qns , clinfo1=' qns - : ', mask1=tmask, ovlap=1) 160 155 CALL prt_ctl(tab2d_1=qsr , clinfo1=' qsr - : ', mask1=tmask, ovlap=1) 161 CALL prt_ctl(tab2d_1=runoff , clinfo1=' runoff : ', mask1=tmask, ovlap=1)162 156 CALL prt_ctl(tab3d_1=tmask , clinfo1=' tmask : ', mask1=tmask, ovlap=1, kdim=jpk) 163 157 CALL prt_ctl(tab3d_1=tn , clinfo1=' sst - : ', mask1=tmask, ovlap=1, kdim=1) 164 158 CALL prt_ctl(tab3d_1=sn , clinfo1=' sss - : ', mask1=tmask, ovlap=1, kdim=1) 165 CALL prt_ctl(tab2d_1= taux , clinfo1=' tau - x: ', mask1=umask, &166 & tab2d_2= tauy , clinfo2=' - y: ', mask2=vmask, ovlap=1)159 CALL prt_ctl(tab2d_1=utau , clinfo1=' tau - u : ', mask1=umask, & 160 & tab2d_2=vtau , clinfo2=' - v : ', mask2=vmask, ovlap=1) 167 161 ENDIF 168 162
Note: See TracChangeset
for help on using the changeset viewer.
