Changeset 9288 for branches/UKMO/dev_r5518_GO6_package_FOAMv14
- Timestamp:
- 2018-01-30T14:59:05+01:00 (7 years ago)
- Location:
- branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/NEMO/OPA_SRC
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90
r6823 r9288 91 91 REAL(wp) :: rn_zqt ! z(q,t) : height of humidity and temperature measurements 92 92 REAL(wp) :: rn_zu ! z(u) : height of wind measurements 93 REAL(wp), PUBLIC :: rn_sfac ! multiplication factor for snow precipitation over sea-ice 93 94 94 95 !! * Substitutions … … 151 152 & sn_wndi, sn_wndj, sn_humi , sn_qsr , & 152 153 & sn_qlw , sn_tair, sn_prec , sn_snow, & 153 & sn_tdif, rn_zqt, rn_zu 154 & sn_tdif, rn_zqt, rn_zu, rn_sfac 154 155 !!--------------------------------------------------------------------- 155 156 ! … … 158 159 ! ! ====================== ! 159 160 ! 161 rn_sfac = 1._wp ! Default to one if missing from namelist 160 162 REWIND( numnam_ref ) ! Namelist namsbc_core in reference namelist : CORE bulk parameters 161 163 READ ( numnam_ref, namsbc_core, IOSTAT = ios, ERR = 901) -
branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90
r8400 r9288 474 474 ! Ensure fsnow is positive (as in CICE routine prepare_forcing) 475 475 IF( iom_use('snowpre') ) CALL iom_put('snowpre',MAX( (1.0-fr_i(:,:))*sprecip(:,:) ,0.0)) !!Joakim edit 476 ztmp(:,:)=MAX(fr_i(:,:)*sprecip(:,:),0.0) 476 IF( kt == nit000 .AND. lwp ) THEN 477 WRITE(numout,*) 'sprecip weight, rn_sfac=', rn_sfac 478 ENDIF 479 ztmp(:,:)=MAX(fr_i(:,:)*rn_sfac*sprecip(:,:),0.0) 477 480 CALL nemo2cice(ztmp,fsnow,'T', 1. ) 478 481 -
branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/NEMO/OPA_SRC/step.F90
r8400 r9288 167 167 ! 168 168 IF( lk_ldfslp ) THEN ! slope of lateral mixing 169 169 CALL eos( tsb, rhd, gdept_0(:,:,:) ) ! before in situ density 170 170 IF( ln_zps .AND. .NOT. ln_isfcav) & 171 171 & CALL zps_hde ( kstp, jpts, tsb, gtsu, gtsv, & ! Partial steps: before horizontal gradient … … 202 202 ! Note that the computation of vertical velocity above, hence "after" sea level 203 203 ! is necessary to compute momentum advection for the rhs of barotropic loop: 204 204 CALL eos ( tsn, rhd, rhop, fsdept_n(:,:,:) ) ! now in situ density for hpg computation 205 205 IF( ln_zps .AND. .NOT. ln_isfcav) & 206 206 & CALL zps_hde ( kstp, jpts, tsn, gtsu, gtsv, & ! Partial steps: before horizontal gradient
Note: See TracChangeset
for help on using the changeset viewer.