Changeset 9739 for NEMO/trunk/cfgs/SPITZ12
- Timestamp:
- 2018-06-05T16:25:10+02:00 (6 years ago)
- Location:
- NEMO/trunk/cfgs/SPITZ12
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/cfgs/SPITZ12/EXPREF/namelist_cfg
r9657 r9739 1 1 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 !! NEMO SPITZ namelist : overwrite some defaultsvalues defined in SHARED/namelist_ref2 !! NEMO/OPA Configuration namelist : overwrite default values defined in SHARED/namelist_ref 3 3 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 4 4 !! SPITZ 1/12 configuration !! … … 23 23 nn_it000 = 1 ! first time step 24 24 nn_itend = 100 ! last time step (std 5475) 25 nn_date0 = 20020101 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)26 25 / 27 26 !----------------------------------------------------------------------- … … 80 79 ln_blk = .true. ! Bulk formulation (T => fill namsbc_blk ) 81 80 ! Sea-ice : 82 nn_ice = 2 ! LIM381 nn_ice = 2 ! SI3 83 82 ! Misc. options of sbc : 84 83 ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr ) … … 190 189 nn_ice_dta = 1 ! = 0, bdy data are equal to the initial state 191 190 ! ! = 1, bdy data are read in 'bdydata .nc' files 192 rn_ice_tem = 267. ! lim3 only: arbitrary temperature of incoming sea ice193 rn_ice_sal = 6. ! lim3 only: -- salinity --194 rn_ice_age = 365. ! lim3 only: -- age --191 rn_ice_tem = 267. ! si3 only: arbitrary temperature of incoming sea ice 192 rn_ice_sal = 6. ! si3 only: -- salinity -- 193 rn_ice_age = 365. ! si3 only: -- age -- 195 194 ! 196 195 nn_rimwidth = 1 ! width of the relaxation zone … … 332 331 !----------------------------------------------------------------------- 333 332 ln_dynspg_ts = .true. ! split-explicit free surface 334 ln_bt_fw = .false. ! Forward integration of barotropic Eqs.333 !!clem ln_bt_fw = .false. ! Forward integration of barotropic Eqs. 335 334 / 336 335 !----------------------------------------------------------------------- -
NEMO/trunk/cfgs/SPITZ12/MY_SRC/sbcblk.F90
r9665 r9739 135 135 # include "vectopt_loop_substitute.h90" 136 136 !!---------------------------------------------------------------------- 137 !! NEMO/O PA 3.7 , NEMO-consortium (2014)137 !! NEMO/OCE 4.0 , NEMO Consortium (2018) 138 138 !! $Id: sbcblk.F90 6416 2016-04-01 12:22:17Z clem $ 139 139 !! Software governed by the CeCILL licence (./LICENSE) … … 496 496 497 497 ! zqla used as temporary array, for rho*U (common term of bulk formulae): 498 zqla(:,:) = zrhoa(:,:) * zU_zu(:,:) 498 zqla(:,:) = zrhoa(:,:) * zU_zu(:,:) * tmask(:,:,1) 499 499 500 500 IF( ABS( rn_zu - rn_zqt) < 0.01_wp ) THEN … … 536 536 & * ( sf(jp_tair)%fnow(:,:,1) - rt0 ) * rcp & 537 537 & + sf(jp_snow)%fnow(:,:,1) * rn_pfac & ! add solid precip heat content at min(Tair,Tsnow) 538 & * ( MIN( sf(jp_tair)%fnow(:,:,1), rt0_snow ) - rt0 ) * cpic * tmask(:,:,1) 538 & * ( MIN( sf(jp_tair)%fnow(:,:,1), rt0_snow ) - rt0 ) * cpic 539 qns(:,:) = qns(:,:) * tmask(:,:,1) 539 540 ! 540 541 #if defined key_si3 541 qns_oce(:,:) = zqlw(:,:) - zqsb(:,:) - zqla(:,:) ! non solar without emp (only needed by LIM3)542 qns_oce(:,:) = zqlw(:,:) - zqsb(:,:) - zqla(:,:) ! non solar without emp (only needed by SI3) 542 543 qsr_oce(:,:) = qsr(:,:) 543 544 #endif … … 551 552 CALL iom_put( "qsr_oce" , qsr ) ! output downward solar heat over the ocean 552 553 CALL iom_put( "qt_oce" , qns+qsr ) ! output total downward heat over the ocean 553 tprecip(:,:) = sf(jp_prec)%fnow(:,:,1) * rn_pfac 554 sprecip(:,:) = sf(jp_snow)%fnow(:,:,1) * rn_pfac 554 tprecip(:,:) = sf(jp_prec)%fnow(:,:,1) * rn_pfac * tmask(:,:,1) ! output total precipitation [kg/m2/s] 555 sprecip(:,:) = sf(jp_snow)%fnow(:,:,1) * rn_pfac * tmask(:,:,1) ! output solid precipitation [kg/m2/s] 555 556 CALL iom_put( 'snowpre', sprecip ) ! Snow 556 557 CALL iom_put( 'precip' , tprecip ) ! Total precipitation … … 842 843 REAL(wp), DIMENSION(jpi,jpj,jpl) :: z_dqlw ! long wave heat sensitivity over ice 843 844 REAL(wp), DIMENSION(jpi,jpj,jpl) :: z_dqsb ! sensible heat sensitivity over ice 844 REAL(wp), DIMENSION(jpi,jpj) :: zevap, zsnw ! evaporation and snw distribution after wind blowing ( LIM3)845 REAL(wp), DIMENSION(jpi,jpj) :: zevap, zsnw ! evaporation and snw distribution after wind blowing (SI3) 845 846 REAL(wp), DIMENSION(jpi,jpj) :: zrhoa 846 847 !!--------------------------------------------------------------------- … … 906 907 END DO 907 908 ! 908 tprecip(:,:) = sf(jp_prec)%fnow(:,:,1) * rn_pfac 909 sprecip(:,:) = sf(jp_snow)%fnow(:,:,1) * rn_pfac 909 tprecip(:,:) = sf(jp_prec)%fnow(:,:,1) * rn_pfac * tmask(:,:,1) ! total precipitation [kg/m2/s] 910 sprecip(:,:) = sf(jp_snow)%fnow(:,:,1) * rn_pfac * tmask(:,:,1) ! solid precipitation [kg/m2/s] 910 911 CALL iom_put( 'snowpre', sprecip ) ! Snow precipitation 911 912 CALL iom_put( 'precip' , tprecip ) ! Total precipitation
Note: See TracChangeset
for help on using the changeset viewer.