Changeset 888 for trunk/NEMO/C1D_SRC/diawri1d.F90
- Timestamp:
- 2008-04-11T19:05:03+02:00 (16 years ago)
- File:
-
- 1 edited
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
Note: See TracChangeset
for help on using the changeset viewer.