Changeset 3489
- Timestamp:
- 2012-10-05T17:10:13+02:00 (11 years ago)
- Location:
- branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/iodef.xml
r3294 r3489 124 124 125 125 <field id="empmr" description="Net Upward Water Flux" unit="kg/m2/s" /> 126 <field id=" empsmr" description="concentration/dilution water flux" unit="kg/m2/s"/>126 <field id="saltflx" description="Downward Salt Flux" unit="PSU/m2/s" /> 127 127 <field id="snowpre" description="Snow precipitation" unit="kg/m2/s" /> 128 128 <field id="runoffs" description="River Runoffs" unit="Kg/m2/s" /> … … 145 145 <field id="qsb_oce" description="Sensible Downward Heat Flux over open ocean" unit="W/m2" /> 146 146 <field id="qla_oce" description="Latent Downward Heat Flux over open ocean" unit="W/m2" /> 147 <field id="qhc_oce" description="Downward Heat Content of E-P over open ocean" unit="W/m2" /> 147 148 <field id="taum_oce" description="wind stress module over open ocean" unit="N/m2" /> 148 149 … … 173 174 <field id="v_imasstr" description="Sea-ice mass transport along j-axis" unit="kg/s" /> 174 175 176 <!-- available if not defined key_vvl --> 177 <field id="emp_x_sst" description="Concentration/Dilution term on SST" unit="kgC/m2/s" /> 178 <field id="emp_x_sss" description="Concentration/Dilution term on SSS" unit="kgPSU/m2/s" /> 175 179 <!-- available key_coupled --> 176 180 <field id="snow_ao_cea" description="Snow over ice-free ocean (cell average)" unit="kg/m2/s" /> … … 1016 1020 <field ref="empmr" name="sowaflup" /> 1017 1021 <field ref="qsr" name="soshfldo" /> 1018 <field ref=" empsmr" name="sowaflcd" />1022 <field ref="saltflx" name="sosfldow" /> 1019 1023 <field ref="qt" name="sohefldo" /> 1020 1024 <field ref="mldr10_1" name="somxl010" /> -
branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/LIM_SRC_2/limsbc_2.F90
r3488 r3489 192 192 193 193 fsbbq(ji,jj) = ( 1.0 - ( ifvt + iflt ) ) * fscmbq(ji,jj) ! store residual heat flux (to put into the ocean at the next time-step) 194 zqhc = ( rdq_snw(ji,jj) + rdq_ice(ji,jj) ) * r1_rdtice ! heat flux due to snow & ice heat content exchange 195 ! 194 zqhc = ( rdq_snw(ji,jj) & 195 & + rdq_ice(ji,jj) * ( 1.- zswitch) ) * r1_rdtice ! heat flux due to snow ( & ice heat content, 196 ! ! if ice/ocean mass exchange active) 196 197 qsr (ji,jj) = zqsr ! solar heat flux 197 qns (ji,jj) = zqns - fdtcn(ji,jj) + zqhc ! non solar heat flux 198 qns (ji,jj) = zqns - fdtcn(ji,jj) + zqhc ! non solar heat flux 199 ! 198 200 ! !------------------------------------------! 199 201 ! ! mass flux at the ocean surface ! … … 223 225 ! salt and mass fluxes : non concentration dilution term in salt flux and add F-M term in volume flux 224 226 sfx (ji,jj) = zfsalt + zswitch * zcd ! salt flux (+ C/D if no ice/ocean mass exchange) 225 emp (ji,jj) = zemp + zemp_snw + ( 1.- zswitch) * zfmm ! mass flux ( - F/M mass flux if noice/ocean mass exchange)227 emp (ji,jj) = zemp + zemp_snw + ( 1.- zswitch) * zfmm ! mass flux (+ F/M mass flux if ice/ocean mass exchange) 226 228 ! 227 229 END DO -
branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90
r3488 r3489 372 372 CALL histdef( nid_T, "sossheig", "Sea Surface Height" , "m" , & ! ssh 373 373 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 374 >>>>>>>>>>> 375 CALL histwrite( nid_T, "sowaflup", it, ( emp-rnf ) , ndim_hT, ndex_hT ) ! upward water flux 376 CALL histwrite( nid_T, "sosfldow", it, sfx , ndim_hT, ndex_hT ) ! downward salt flux 377 ! (includes virtual salt flux beneath ice 378 ! in linear free surface case) 374 CALL histdef( nid_T, "sowaflup", "Net Upward Water Flux" , "Kg/m2/s", & ! (emp-rnf) 375 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 376 CALL histdef( nid_T, "sosfldow", "downward salt flux" , "PSU/m2/s", & ! sfx 377 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 379 378 #if ! defined key_vvl 380 zw2d(:,:) = emp (:,:) * tsn(:,:,1,jp_tem) * tmask(:,:,1) 381 CALL histwrite( nid_T, "sosst_cd", it, zw2d, ndim_hT, ndex_hT ) ! c/d term on sst 382 zw2d(:,:) = emp (:,:) * tsn(:,:,1,jp_sal) * tmask(:,:,1) 383 CALL histwrite( nid_T, "sosss_cd", it, zw2d, ndim_hT, ndex_hT ) ! c/d term on sss 384 #endif 385 <<<<<<<<<< 386 CALL histdef( nid_T, "sowaflup", "Net Upward Water Flux" , "Kg/m2/s", & ! (emp-rnf) 387 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 388 CALL histdef( nid_T, "sofldow", "downward salt flux" , "PSU/m2/s", & ! sfx 389 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 390 CALL histdef( nid_T, "sosalflx", "Surface Salt Flux" , "Kg/m2/s", & ! (emps-rnf) * sn 391 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 379 CALL histdef( nid_T, "sosst_cd", "Concentration/Dilution term on temperature" & ! emp * tsn(:,:,1,jp_tem) 380 & , "KgC/m2/s", & ! sosst_cd 381 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 382 CALL histdef( nid_T, "sosss_cd", "Concentration/Dilution term on salinity" & ! emp * tsn(:,:,1,jp_sal) 383 & , "KgPSU/m2/s",& ! sosss_cd 384 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) 385 #endif 392 386 CALL histdef( nid_T, "sohefldo", "Net Downward Heat Flux" , "W/m2" , & ! qns + qsr 393 387 & jpi, jpj, nh_T, 1 , 1, 1 , -99 , 32, clop, zsto, zout ) -
branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/OPA_SRC/DIA/diawri_dimg.h90
r3402 r3489 54 54 !! level 14: qct(:,:) equivalent flux due to treshold SST 55 55 !! level 15: fbt(:,:) feedback term . 56 !! level 16: ( sfx (:,:) - rnf(:,:) ) concentration/dilution water flux 56 !! level 16: ( emp * sss ) concentration/dilution term on salinity 57 !! level 17: ( emp * sst ) concentration/dilution term on temperature 57 58 !! level 17: fsalt(:,:) Ice=>ocean net freshwater 58 59 !! level 18: gps(:,:) the surface pressure (m). … … 62 63 !! History: OPA ! 1997-02 ( Clipper Group ) dimg files 63 64 !! - ! 2003-12 ( J.M. Molines) f90, mpp output for OPA9.0 64 !! NEMO 1.0 ! 2005-05 (S. Theetten) add sfxfsalt move gps spgu spgv 2 lines below65 !! NEMO 1.0 ! 2005-05 (S. Theetten) add emps fsalt move gps spgu spgv 2 lines below 65 66 !! - ! 2005-11 (V. Garnier) Surface pressure gradient organization 66 67 !!---------------------------------------------------------------------- … … 107 108 108 109 109 inbsel = 1 7110 inbsel = 18 110 111 111 112 IF( inbsel > jpk ) THEN … … 172 173 ! fsel(:,:,14) = fsel(:,:,14) + qct(:,:) 173 174 ! fsel(:,:,15) = fsel(:,:,15) + fbt(:,:) 174 fsel(:,:,16) = fsel(:,:,16) + ( sfx (:,:)-rnf(:,:) ) 175 fsel(:,:,16) = fsel(:,:,16) + ( emp(:,:)*tsn(:,:,1,jp_sal) ) 176 fsel(:,:,17) = fsel(:,:,17) + ( emp(:,:)*tsn(:,:,1,jp_tem) ) 175 177 ! 176 178 ! Output of dynamics and tracer fields and selected fields … … 240 242 ! fsel(:,:,14) = qct(:,:) 241 243 ! fsel(:,:,15) = fbt(:,:) 242 fsel(:,:,16) = ( sfx (:,:)-rnf(:,:) ) * tmask(:,:,1) 244 fsel(:,:,16) = ( emp(:,:)-tsn(:,:,1,jp_sal) ) * tmask(:,:,1) 245 fsel(:,:,17) = ( emp(:,:)-tsn(:,:,1,jp_tem) ) * tmask(:,:,1) 243 246 ! 244 247 ! qct(:,:) = 0._wp -
branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90
r3402 r3489 392 392 CALL iom_put( "qsb_oce", - zqsb ) ! output downward sensible heat over the ocean 393 393 CALL iom_put( "qla_oce", - zqla ) ! output downward latent heat over the ocean 394 CALL iom_put( "qhc_oce", qns-zqlw+zqsb+zqla ) ! output downward heat content of E-P over the ocean 394 395 CALL iom_put( "qns_oce", qns ) ! output downward non solar heat over the ocean 395 396 !
Note: See TracChangeset
for help on using the changeset viewer.