Changeset 719 for trunk/NEMO/OPA_SRC/DIA/diawri.F90
- Timestamp:
- 2007-10-16T16:59:56+02:00 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DIA/diawri.F90
- Property svn:keywords changed from Id to Author Date Id Revision
r717 r719 14 14 USE sol_oce ! solver variables 15 15 USE ice_oce ! ice variables 16 USE sbc_oce ! surface boundary condition: ocean17 USE sbc_ice ! surface boundary condition: ice18 USE sbcssr ! restoring term toward SST/SSS climatology19 16 USE phycst ! physical constants 20 17 USE ocfzpt ! ocean freezing point 18 USE ocesbc ! surface thermohaline fluxes 19 USE taumod ! surface stress 20 USE flxrnf ! ocean runoffs 21 21 USE zdfmxl ! mixed layer 22 22 USE daymod ! calendar … … 27 27 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 28 28 USE in_out_manager ! I/O manager 29 USE flx_oce ! sea-ice/ocean forcings variables 29 30 USE diadimg ! dimg direct access file format output 30 31 USE ioipsl … … 53 54 !!---------------------------------------------------------------------- 54 55 !! OPA 9.0 , LOCEAN-IPSL (2005) 55 !! $ Id$56 !! $Header$ 56 57 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 57 58 !!---------------------------------------------------------------------- … … 244 245 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 245 246 #endif 246 !!$#if ! defined key_dynspg_rl && defined key_ice_lim247 !!$! sowaflup = sowaflep + sorunoff + sowafldp + a term associated to248 !!$! internal damping to Levitus that can be diagnosed from others249 !!$! sowaflcd = sowaflep + sorunoff + sowafldp + iowaflup250 !!$CALL histdef( nid_T, "iowaflup", "Ice=>ocean net freshwater" , "kg/m2/s", & ! fsalt251 !!$& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )252 !!$CALL histdef( nid_T, "sowaflep", "atmos=>ocean net freshwater" , "kg/m2/s", & ! fmass253 !!$& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )254 !!$#endif247 #if ! defined key_dynspg_rl && defined key_ice_lim 248 ! sowaflup = sowaflep + sorunoff + sowafldp + a term associated to 249 ! internal damping to Levitus that can be diagnosed from others 250 ! sowaflcd = sowaflep + sorunoff + sowafldp + iowaflup 251 CALL histdef( nid_T, "iowaflup", "Ice=>ocean net freshwater" , "kg/m2/s", & ! fsalt 252 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 253 CALL histdef( nid_T, "sowaflep", "atmos=>ocean net freshwater" , "kg/m2/s", & ! fmass 254 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 255 #endif 255 256 CALL histdef( nid_T, "sowaflup", "Net Upward Water Flux" , "Kg/m2/s", & ! emp 256 257 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 257 !!$CALL histdef( nid_T, "sorunoff", "Runoffs" , "Kg/m2/s", & ! runoffs258 !!$& jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout )258 CALL histdef( nid_T, "sorunoff", "Runoffs" , "Kg/m2/s", & ! runoffs 259 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 259 260 CALL histdef( nid_T, "sowaflcd", "concentration/dilution water flux" , "kg/m2/s", & ! emps 260 261 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 261 262 CALL histdef( nid_T, "sosalflx", "Surface Salt Flux" , "Kg/m2/s", & ! emps * sn 262 263 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 263 CALL histdef( nid_T, "sohefldo", "Net Downward Heat Flux" , "W/m2" , & ! q ns + qsr264 CALL histdef( nid_T, "sohefldo", "Net Downward Heat Flux" , "W/m2" , & ! qt 264 265 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 265 266 CALL histdef( nid_T, "soshfldo", "Shortwave Radiation" , "W/m2" , & ! qsr … … 333 334 #endif 334 335 ! !!! nid_U : 2D 335 CALL histdef( nid_U, "sozotaux", "Wind Stress along i-axis" , "N/m2" , & ! utau336 CALL histdef( nid_U, "sozotaux", "Wind Stress along i-axis" , "N/m2" , & ! taux 336 337 & jpi, jpj, nh_U, 1 , 1, 1 , - 99, 32, clop, zsto, zout ) 337 338 #if defined key_dynspg_rl … … 350 351 #endif 351 352 ! !!! nid_V : 2D 352 CALL histdef( nid_V, "sometauy", "Wind Stress along j-axis" , "N/m2" , & ! vtau353 CALL histdef( nid_V, "sometauy", "Wind Stress along j-axis" , "N/m2" , & ! tauy 353 354 & jpi, jpj, nh_V, 1 , 1, 1 , - 99, 32, clop, zsto, zout ) 354 355 #if defined key_dynspg_rl … … 422 423 CALL histwrite( nid_T, "sossheig", it, sshn , ndim_hT, ndex_hT ) ! sea surface height 423 424 #endif 424 !!$#if ! defined key_dynspg_rl && defined key_ice_lim425 !!$CALL histwrite( nid_T, "iowaflup", it, fsalt(:,:) , ndim_hT, ndex_hT ) ! ice=>ocean water flux426 !!$CALL histwrite( nid_T, "sowaflep", it, fmass(:,:) , ndim_hT, ndex_hT ) ! atmos=>ocean water flux427 !!$#endif425 #if ! defined key_dynspg_rl && defined key_ice_lim 426 CALL histwrite( nid_T, "iowaflup", it, fsalt(:,:) , ndim_hT, ndex_hT ) ! ice=>ocean water flux 427 CALL histwrite( nid_T, "sowaflep", it, fmass(:,:) , ndim_hT, ndex_hT ) ! atmos=>ocean water flux 428 #endif 428 429 CALL histwrite( nid_T, "sowaflup", it, emp , ndim_hT, ndex_hT ) ! upward water flux 429 !!$CALL histwrite( nid_T, "sorunoff", it, runoff , ndim_hT, ndex_hT ) ! runoff430 CALL histwrite( nid_T, "sorunoff", it, runoff , ndim_hT, ndex_hT ) ! runoff 430 431 CALL histwrite( nid_T, "sowaflcd", it, emps , ndim_hT, ndex_hT ) ! c/d water flux 431 432 zw2d(:,:) = emps(:,:) * sn(:,:,1) * tmask(:,:,1) 432 433 CALL histwrite( nid_T, "sosalflx", it, zw2d , ndim_hT, ndex_hT ) ! c/d salt flux 433 CALL histwrite( nid_T, "sohefldo", it, q ns + qsr, ndim_hT, ndex_hT ) ! total heat flux434 CALL histwrite( nid_T, "sohefldo", it, qt , ndim_hT, ndex_hT ) ! total heat flux 434 435 CALL histwrite( nid_T, "soshfldo", it, qsr , ndim_hT, ndex_hT ) ! solar heat flux 435 436 CALL histwrite( nid_T, "somxl010", it, hmlp , ndim_hT, ndex_hT ) ! mixed layer depth … … 474 475 CALL histwrite( nid_U, "vozoeivu", it, u_eiv , ndim_U , ndex_U ) ! i-eiv current 475 476 #endif 476 CALL histwrite( nid_U, "sozotaux", it, utau, ndim_hU, ndex_hU ) ! i-wind stress477 CALL histwrite( nid_U, "sozotaux", it, taux , ndim_hU, ndex_hU ) ! i-wind stress 477 478 #if defined key_dynspg_rl 478 479 CALL lbc_lnk( spgu, 'U', -1. ) … … 485 486 CALL histwrite( nid_V, "vomeeivv", it, v_eiv , ndim_V , ndex_V ) ! j-eiv current 486 487 #endif 487 CALL histwrite( nid_V, "sometauy", it, vtau, ndim_hV, ndex_hV ) ! j-wind stress488 CALL histwrite( nid_V, "sometauy", it, tauy , ndim_hV, ndex_hV ) ! j-wind stress 488 489 #if defined key_dynspg_rl 489 490 CALL lbc_lnk( spgv, 'V', -1. ) … … 639 640 640 641 ! Write all fields on T grid 641 CALL histwrite( id_i, "votemper", 1, tn 642 CALL histwrite( id_i, "vosaline", 1, sn 643 #if defined key_dynspg_rl 644 CALL histwrite( id_i, "sobarstf", 1, bsfn 642 CALL histwrite( id_i, "votemper", 1, tn , jpi*jpj*jpk, idex ) ! now temperature 643 CALL histwrite( id_i, "vosaline", 1, sn , jpi*jpj*jpk, idex ) ! now salinity 644 #if defined key_dynspg_rl 645 CALL histwrite( id_i, "sobarstf", 1, bsfn , jpi*jpj , idex ) ! barotropic streamfunction 645 646 #else 646 CALL histwrite( id_i, "sossheig", 1, sshn 647 #endif 648 CALL histwrite( id_i, "vozocrtx", 1, un 649 CALL histwrite( id_i, "vomecrty", 1, vn 650 CALL histwrite( id_i, "vovecrtz", 1, wn 651 CALL histwrite( id_i, "sowaflup", 1, emp 652 CALL histwrite( id_i, "sohefldo", 1, q sr + qns, jpi*jpj , idex ) ! total heat flux653 CALL histwrite( id_i, "soshfldo", 1, qsr , jpi*jpj , idex ) ! solarheat flux654 CALL histwrite( id_i, "soicecov", 1, freeze 655 CALL histwrite( id_i, "sozotaux", 1, utau, jpi*jpj , idex ) ! i-wind stress656 CALL histwrite( id_i, "sometauy", 1, vtau, jpi*jpj , idex ) ! j-wind stress647 CALL histwrite( id_i, "sossheig", 1, sshn , jpi*jpj , idex ) ! sea surface height 648 #endif 649 CALL histwrite( id_i, "vozocrtx", 1, un , jpi*jpj*jpk, idex ) ! now i-velocity 650 CALL histwrite( id_i, "vomecrty", 1, vn , jpi*jpj*jpk, idex ) ! now j-velocity 651 CALL histwrite( id_i, "vovecrtz", 1, wn , jpi*jpj*jpk, idex ) ! now k-velocity 652 CALL histwrite( id_i, "sowaflup", 1, emp , jpi*jpj , idex ) ! freshwater budget 653 CALL histwrite( id_i, "sohefldo", 1, qt , jpi*jpj , idex ) ! total heat flux 654 CALL histwrite( id_i, "soshfldo", 1, qsr , jpi*jpj , idex ) ! total heat flux 655 CALL histwrite( id_i, "soicecov", 1, freeze, jpi*jpj , idex ) ! ice cover 656 CALL histwrite( id_i, "sozotaux", 1, taux , jpi*jpj , idex ) ! i-wind stress 657 CALL histwrite( id_i, "sometauy", 1, tauy , jpi*jpj , idex ) ! j-wind stress 657 658 658 659 ! 3. Close the file
Note: See TracChangeset
for help on using the changeset viewer.