Changeset 5381
- Timestamp:
- 2015-06-09T00:48:54+02:00 (9 years ago)
- Location:
- branches/2015/dev_r5218_CNRS17_coupling/NEMOGCM
- Files:
-
- 9 edited
-
CONFIG/SHARED/field_def.xml (modified) (4 diffs)
-
NEMO/LIM_SRC_3/limsbc.F90 (modified) (1 diff)
-
NEMO/OPA_SRC/IOM/restart.F90 (modified) (4 diffs)
-
NEMO/OPA_SRC/SBC/sbcblk_clio.F90 (modified) (1 diff)
-
NEMO/OPA_SRC/SBC/sbcblk_core.F90 (modified) (1 diff)
-
NEMO/OPA_SRC/SBC/sbccpl.F90 (modified) (10 diffs)
-
NEMO/OPA_SRC/SBC/sbcmod.F90 (modified) (1 diff)
-
NEMO/OPA_SRC/SBC/sbcrnf.F90 (modified) (7 diffs)
-
NEMO/SAS_SRC/sbcssm.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5218_CNRS17_coupling/NEMOGCM/CONFIG/SHARED/field_def.xml
r5352 r5381 192 192 <field id="qsb_oce" long_name="Sensible Downward Heat Flux over open ocean" unit="W/m2" /> 193 193 <field id="qla_oce" long_name="Latent Downward Heat Flux over open ocean" unit="W/m2" /> 194 <field id="q hc_oce"long_name="Downward Heat Content of E-P over open ocean" unit="W/m2" />194 <field id="qemp_oce" long_name="Downward Heat Content of E-P over open ocean" unit="W/m2" /> 195 195 <field id="taum_oce" long_name="wind stress module over open ocean" unit="N/m2" /> 196 196 … … 241 241 <field id="uice_ipa" long_name="Ice velocity along i-axis at I-point (ice presence average)" unit="m/s" /> 242 242 <field id="vice_ipa" long_name="Ice velocity along j-axis at I-point (ice presence average)" unit="m/s" /> 243 243 244 244 <field id="utau_ice" long_name="Wind stress along i-axis over the ice at i-point" unit="N/m2" /> 245 245 <field id="vtau_ice" long_name="Wind stress along j-axis over the ice at i-point" unit="N/m2" /> … … 251 251 252 252 <field id="iceconc" long_name="ice concentration" unit="%" /> 253 <field id="uice_ipa" long_name="Ice velocity along i-axis at I-point (ice presence average)" unit="m/s" />254 <field id="vice_ipa" long_name="Ice velocity along j-axis at I-point (ice presence average)" unit="m/s" />255 253 <field id="isst" long_name="sea surface temperature" unit="degC" /> 256 254 <field id="isss" long_name="sea surface salinity" unit="psu" /> … … 262 260 <field id="qns_ice" long_name="non-solar heat flux at ice surface" unit="W/m2" /> 263 261 <field id="qtr_ice" long_name="solar heat flux transmitted thru the ice" unit="W/m2" /> 262 <field id="qemp_ice" long_name="Downward Heat Content of E-P over ice" unit="W/m2" /> 264 263 <field id="utau_ice" long_name="Wind stress along i-axis over the ice at i-point" unit="N/m2" /> 265 264 <field id="vtau_ice" long_name="Wind stress along j-axis over the ice at i-point" unit="N/m2" /> -
branches/2015/dev_r5218_CNRS17_coupling/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90
r5376 r5381 110 110 111 111 ! make calls for heat fluxes before it is modified 112 !!$IF( iom_use('qsr_oce') ) CALL iom_put( "qsr_oce" , qsr_oce(:,:) * pfrld(:,:) ) ! solar flux at ocean surface113 !!$IF( iom_use('qns_oce') ) CALL iom_put( "qns_oce" , qns_oce(:,:) * pfrld(:,:) + qemp_oce(:,:) ) ! non-solar flux at ocean surface112 IF( iom_use('qsr_oce') ) CALL iom_put( "qsr_oce" , qsr_oce(:,:) * pfrld(:,:) ) ! solar flux at ocean surface 113 IF( iom_use('qns_oce') ) CALL iom_put( "qns_oce" , qns_oce(:,:) * pfrld(:,:) + qemp_oce(:,:) ) ! non-solar flux at ocean surface 114 114 IF( iom_use('qsr_ice') ) CALL iom_put( "qsr_ice" , SUM( qsr_ice(:,:,:) * a_i_b(:,:,:), dim=3 ) ) ! solar flux at ice surface 115 115 IF( iom_use('qns_ice') ) CALL iom_put( "qns_ice" , SUM( qns_ice(:,:,:) * a_i_b(:,:,:), dim=3 ) + qemp_ice(:,:) ) ! non-solar flux at ice surface 116 116 IF( iom_use('qtr_ice') ) CALL iom_put( "qtr_ice" , SUM( ftr_ice(:,:,:) * a_i_b(:,:,:), dim=3 ) ) ! solar flux transmitted thru ice 117 !!$IF( iom_use('qt_oce' ) ) CALL iom_put( "qt_oce" , ( qsr_oce(:,:) + qns_oce(:,:) ) * pfrld(:,:) + qemp_oce(:,:) )117 IF( iom_use('qt_oce' ) ) CALL iom_put( "qt_oce" , ( qsr_oce(:,:) + qns_oce(:,:) ) * pfrld(:,:) + qemp_oce(:,:) ) 118 118 IF( iom_use('qt_ice' ) ) CALL iom_put( "qt_ice" , SUM( ( qns_ice(:,:,:) + qsr_ice(:,:,:) ) & 119 119 & * a_i_b(:,:,:), dim=3 ) + qemp_ice(:,:) ) 120 IF( iom_use('qemp_oce' ) ) CALL iom_put( "qemp_oce" , qemp_oce(:,:) ) 121 IF( iom_use('qemp_ice' ) ) CALL iom_put( "qemp_ice" , qemp_ice(:,:) ) 120 122 121 123 ! pfrld is the lead fraction at the previous time step (actually between TRP and THD) -
branches/2015/dev_r5218_CNRS17_coupling/NEMOGCM/NEMO/OPA_SRC/IOM/restart.F90
r5357 r5381 24 24 USE trdmxl_oce ! ocean active mixed layer tracers trends variables 25 25 USE divcur ! hor. divergence and curl (div & cur routines) 26 USE sbc_ice, ONLY : lk_lim327 USE sbc_oce, ONLY : nn_components, jp_iam_opa28 26 29 27 IMPLICIT NONE … … 121 119 CALL iom_rstput( kt, nitrst, numrow, 'hdivb' , hdivb ) 122 120 CALL iom_rstput( kt, nitrst, numrow, 'sshb' , sshb ) 123 !124 IF( lk_lim3 .OR. ( nn_components == jp_iam_opa ) ) CALL iom_rstput( kt, nitrst, numrow, 'fse3t_b', fse3t_b(:,:,:) )125 121 ! 126 122 CALL iom_rstput( kt, nitrst, numrow, 'un' , un ) ! now fields … … 212 208 CALL iom_get( numror, jpdom_autoglo, 'hdivb' , hdivb ) 213 209 CALL iom_get( numror, jpdom_autoglo, 'sshb' , sshb ) 214 ! EM Attention Ceci doit etre reimplemente correctement215 IF( lk_lim3 .OR. ( nn_components == jp_iam_opa ) ) CALL iom_get( numror, jpdom_autoglo, 'fse3t_b', fse3t_b(:,:,:) )216 !clem CALL iom_get( numror, jpdom_autoglo, 'fse3t_b', fse3t_b(:,:,:) )217 210 ELSE 218 211 neuler = 0 … … 257 250 ENDIF 258 251 259 IF( ( lk_lim3 .OR. ( nn_components == jp_iam_opa ) ) .AND. .NOT. lk_vvl ) THEN260 DO jk = 1, jpk261 fse3t_b(:,:,jk) = fse3t_n(:,:,jk)262 END DO263 ENDIF264 265 252 ENDIF 266 253 ! -
branches/2015/dev_r5218_CNRS17_coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_clio.F90
r5364 r5381 390 390 ! NB: if sea-ice model, the snow precip are computed and the associated heat is added to qns (see blk_ice_clio) 391 391 392 CALL iom_put( "qlw_oce", zqlw ) ! output downward longwave heat over the ocean 393 CALL iom_put( "qsb_oce", - zqsb ) ! output downward sensible heat over the ocean 394 CALL iom_put( "qla_oce", - zqla ) ! output downward latent heat over the ocean 395 CALL iom_put( "qns_oce", qns ) ! output downward non solar heat over the ocean 392 IF ( nn_ice == 0 ) THEN 393 CALL iom_put( "qlw_oce" , zqlw ) ! output downward longwave heat over the ocean 394 CALL iom_put( "qsb_oce" , - zqsb ) ! output downward sensible heat over the ocean 395 CALL iom_put( "qla_oce" , - zqla ) ! output downward latent heat over the ocean 396 CALL iom_put( "qemp_oce", qns-zqlw+zqsb+zqla ) ! output downward heat content of E-P over the ocean 397 CALL iom_put( "qns_oce" , qns ) ! output downward non solar heat over the ocean 398 CALL iom_put( "qsr_oce" , qsr ) ! output downward solar heat over the ocean 399 CALL iom_put( "qt_oce" , qns+qsr ) ! output total downward heat over the ocean 400 ENDIF 396 401 397 402 IF(ln_ctl) THEN -
branches/2015/dev_r5218_CNRS17_coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90
r5376 r5381 400 400 #endif 401 401 ! 402 CALL iom_put( "qlw_oce", zqlw ) ! output downward longwave heat over the ocean 403 CALL iom_put( "qsb_oce", - zqsb ) ! output downward sensible heat over the ocean 404 CALL iom_put( "qla_oce", - zqla ) ! output downward latent heat over the ocean 405 CALL iom_put( "qhc_oce", qns-zqlw+zqsb+zqla ) ! output downward heat content of E-P over the ocean 406 CALL iom_put( "qns_oce", qns ) ! output downward non solar heat over the ocean 407 CALL iom_put( "qsr_oce", qsr ) ! output downward solar heat over the ocean 408 CALL iom_put( "qt_oce" , qns+qsr ) ! output total downward heat over the ocean 402 IF ( nn_ice == 0 ) THEN 403 CALL iom_put( "qlw_oce" , zqlw ) ! output downward longwave heat over the ocean 404 CALL iom_put( "qsb_oce" , - zqsb ) ! output downward sensible heat over the ocean 405 CALL iom_put( "qla_oce" , - zqla ) ! output downward latent heat over the ocean 406 CALL iom_put( "qemp_oce", qns-zqlw+zqsb+zqla ) ! output downward heat content of E-P over the ocean 407 CALL iom_put( "qns_oce" , qns ) ! output downward non solar heat over the ocean 408 CALL iom_put( "qsr_oce" , qsr ) ! output downward solar heat over the ocean 409 CALL iom_put( "qt_oce" , qns+qsr ) ! output total downward heat over the ocean 410 ENDIF 409 411 ! 410 412 IF(ln_ctl) THEN -
branches/2015/dev_r5218_CNRS17_coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90
r5376 r5381 43 43 USE eosbn2 44 44 USE traqsr , ONLY : fraqsr_1lev 45 USE sbcrnf , ONLY : l_rnfcpl, ln_rnf_emp 45 46 #if defined key_cpl_carbon_cycle 46 47 USE p4zflx, ONLY : oce_co2 … … 401 402 ! ! ------------------------- ! 402 403 srcv(jpr_rnf )%clname = 'O_Runoff' ; IF( TRIM( sn_rcv_rnf%cldes ) == 'coupled' ) srcv(jpr_rnf)%laction = .TRUE. 403 ! This isn't right - really just want ln_rnf_emp changed 404 ! IF( TRIM( sn_rcv_rnf%cldes ) == 'climato' ) THEN ; ln_rnf = .TRUE. 405 ! ELSE ; ln_rnf = .FALSE. 406 ! ENDIF 404 l_rnfcpl = srcv(jpr_rnf)%laction 405 ! 407 406 srcv(jpr_cal )%clname = 'OCalving' ; IF( TRIM( sn_rcv_cal%cldes ) == 'coupled' ) srcv(jpr_cal)%laction = .TRUE. 408 407 … … 536 535 ! Vectors: change of sign at north fold ONLY if on the local grid 537 536 srcv(jpr_ocx1:jpr_ocy1)%nsgn = -1. 538 ! Change first letter to couple with atmosphere if already coupled with sea_ice 537 ! Change first letter to couple with atmosphere if already coupled OPA 538 ! this is nedeed as each variable name used in the namcouple must be unique: 539 ! for example O_Runoff received by OPA from SAS and therefore O_Runoff received by SAS from the Atmosphere 539 540 DO jn = 1, jprcv 540 541 IF ( srcv(jn)%clname(1:1) == "O" ) srcv(jn)%clname = "S"//srcv(jn)%clname(2:LEN(srcv(jn)%clname)) … … 777 778 xcplmask(:,:,0) = 1. - SUM( xcplmask(:,:,1:nn_cplmodel), dim = 3 ) 778 779 ! 779 IF( ln_dm2dc .AND. &780 IF( ln_dm2dc .AND. ln_cpl .AND. & 780 781 & ( cpl_freq( 'O_QsrOce' ) + cpl_freq( 'O_QsrMix' ) + cpl_freq( 'S_QsrOce' ) + cpl_freq( 'S_QsrMix' ) ) /= 86400 ) & 781 782 & CALL ctl_stop( 'sbc_cpl_init: diurnal cycle reconstruction (ln_dm2dc) needs daily couping for solar radiation' ) … … 1066 1067 ! 1067 1068 ! ! runoffs and calving (added in emp) 1068 IF( srcv(jpr_rnf)%laction ) zemp(:,:) = zemp(:,:) - frcv(jpr_rnf)%z3(:,:,1) 1069 IF( srcv(jpr_cal)%laction ) zemp(:,:) = zemp(:,:) - frcv(jpr_cal)%z3(:,:,1) 1069 IF( srcv(jpr_rnf)%laction ) THEN 1070 IF( ln_rnf_emp ) THEN 1071 zemp(:,:) = zemp(:,:) - frcv(jpr_rnf)%z3(:,:,1) 1072 CALL iom_put( 'runoffs', zemp(:,:) ) ! rivers 1073 ELSE 1074 rnf(:,:) = frcv(jpr_rnf)%z3(:,:,1) 1075 ENDIF 1076 ENDIF 1077 IF( srcv(jpr_cal)%laction ) zemp(:,:) = zemp(:,:) - frcv(jpr_cal)%z3(:,:,1) 1070 1078 1071 1079 IF( ln_mixcpl ) THEN ; emp(:,:) = emp(:,:) * xcplmask(:,:,0) + zemp(:,:) * zmsk(:,:) … … 1094 1102 ELSE ; zqsr(:,:) = 0._wp 1095 1103 ENDIF 1096 IF( ln_dm2dc .AND. nn_components /= jp_iam_opa) zqsr(:,:) = sbc_dcy( zqsr ) ! modify qsr to include the diurnal cycle1104 IF( ln_dm2dc .AND. ln_cpl ) zqsr(:,:) = sbc_dcy( zqsr ) ! modify qsr to include the diurnal cycle 1097 1105 IF( ln_mixcpl ) THEN ; qsr(:,:) = qsr(:,:) * xcplmask(:,:,0) + zqsr(:,:) * zmsk(:,:) 1098 1106 ELSE ; qsr(:,:) = zqsr(:,:) … … 1422 1430 ! ! runoffs and calving (put in emp_tot) 1423 1431 IF( srcv(jpr_rnf)%laction ) THEN 1424 zemp_tot(:,:) = zemp_tot(:,:) - frcv(jpr_rnf)%z3(:,:,1) 1425 CALL iom_put( 'runoffs' , frcv(jpr_rnf)%z3(:,:,1) ) ! rivers 1426 IF( iom_use('hflx_rnf_cea') ) & 1427 CALL iom_put( 'hflx_rnf_cea' , frcv(jpr_rnf)%z3(:,:,1) * zcptn(:,:) ) ! heat flux from rivers 1432 IF( ln_rnf_emp ) THEN 1433 zemp_tot(:,:) = zemp_tot(:,:) - frcv(jpr_rnf)%z3(:,:,1) 1434 CALL iom_put( 'runoffs' , frcv(jpr_rnf)%z3(:,:,1) ) ! rivers 1435 IF( iom_use('hflx_rnf_cea') ) CALL iom_put( 'hflx_rnf_cea' , frcv(jpr_rnf)%z3(:,:,1) * zcptn(:,:) ) ! river heat flx 1436 ELSE 1437 rnf(:,:) = frcv(jpr_rnf)%z3(:,:,1) 1438 ENDIF 1428 1439 ENDIF 1429 1440 IF( srcv(jpr_cal)%laction ) THEN … … 1546 1557 ENDDO 1547 1558 qprec_ice(:,:) = qprec_ice(:,:) * xcplmask(:,:,0) + zqprec_ice(:,:)* zmsk(:,:) 1548 qemp_oce (:,:) = qemp_oce(:,:) * xcplmask(:,:,0) +zqemp_oce(:,:)* zmsk(:,:)1559 qemp_oce (:,:) = qemp_oce(:,:) * xcplmask(:,:,0) + zqemp_oce(:,:)* zmsk(:,:) 1549 1560 !!clem evap_ice(:,:) = evap_ice(:,:) * xcplmask(:,:,0) 1550 1561 ELSE … … 1560 1571 #else 1561 1572 1562 ! clem: this formulation is certainly wrong. 1573 ! clem: this formulation is certainly wrong... but better than it was... 1563 1574 zqns_tot(:,:) = zqns_tot(:,:) & ! zqns_tot update over free ocean with: 1564 1575 & - ztmp(:,:) & ! remove the latent heat flux of solid precip. melting … … 1619 1630 & + palbi (:,:,1) * zicefr(:,:) ) ) 1620 1631 END SELECT 1621 IF( ln_dm2dc ) THEN ! modify qsr to include the diurnal cycle1632 IF( ln_dm2dc .AND. ln_cpl ) THEN ! modify qsr to include the diurnal cycle 1622 1633 zqsr_tot(:,: ) = sbc_dcy( zqsr_tot(:,: ) ) 1623 1634 DO jl=1,jpl -
branches/2015/dev_r5218_CNRS17_coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90
r5376 r5381 225 225 IF( ln_dm2dc ) nday_qsr = -1 ! initialisation flag 226 226 227 IF( ln_dm2dc .AND. .NOT.( ln_flx .OR. ln_blk_core ) ) &227 IF( ln_dm2dc .AND. .NOT.( ln_flx .OR. ln_blk_core ) .AND. nn_components /= jp_iam_opa ) & 228 228 & CALL ctl_stop( 'diurnal cycle into qsr field from daily values requires a flux or core-bulk formulation' ) 229 229 -
branches/2015/dev_r5218_CNRS17_coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcrnf.F90
r4990 r5381 51 51 REAL(wp) , PUBLIC :: rn_rfact !: multiplicative factor for runoff 52 52 53 LOGICAL , PUBLIC :: l_rnfcpl = .false. ! runoffs recieved from oasis 54 53 55 INTEGER , PUBLIC :: nkrnf = 0 !: nb of levels over which Kz is increased at river mouths 54 56 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: rnfmsk !: river mouth mask (hori.) … … 120 122 ! !-------------------! 121 123 ! 122 123 IF( ln_rnf_tem) CALL fld_read ( kt, nn_fsbc, sf_t_rnf ) ! idem for runoffs temperature if required124 IF( ln_rnf_sal) CALL fld_read ( kt, nn_fsbc, sf_s_rnf ) ! idem for runoffs salinity if required124 IF( .NOT. l_rnfcpl ) CALL fld_read ( kt, nn_fsbc, sf_rnf ) ! Read Runoffs data and provide it at kt 125 IF( ln_rnf_tem ) CALL fld_read ( kt, nn_fsbc, sf_t_rnf ) ! idem for runoffs temperature if required 126 IF( ln_rnf_sal ) CALL fld_read ( kt, nn_fsbc, sf_s_rnf ) ! idem for runoffs salinity if required 125 127 ! 126 128 ! Runoff reduction only associated to the ORCA2_LIM configuration 127 129 ! when reading the NetCDF file runoff_1m_nomask.nc 128 IF( cp_cfg == 'orca' .AND. jp_cfg == 2 ) THEN130 IF( cp_cfg == 'orca' .AND. jp_cfg == 2 .AND. .NOT. l_rnfcpl ) THEN 129 131 WHERE( 40._wp < gphit(:,:) .AND. gphit(:,:) < 65._wp ) 130 132 sf_rnf(1)%fnow(:,:,1) = 0.85 * sf_rnf(1)%fnow(:,:,1) … … 134 136 IF( MOD( kt - 1, nn_fsbc ) == 0 ) THEN 135 137 ! 136 rnf(:,:) = rn_rfact * ( sf_rnf(1)%fnow(:,:,1) ) ! updated runoff value at time step kt138 IF( .NOT. l_rnfcpl ) rnf(:,:) = rn_rfact * ( sf_rnf(1)%fnow(:,:,1) ) ! updated runoff value at time step kt 137 139 ! 138 140 ! ! set temperature & salinity content of runoffs … … 152 154 IF( ln_rnf_sal ) rnf_tsc(:,:,jp_sal) = ( sf_s_rnf(1)%fnow(:,:,1) ) * rnf(:,:) * r1_rau0 153 155 ! ! else use S=0 for runoffs (done one for all in the init) 154 IF ( ANY( rnf(:,:) < 0._wp ) ) z_err=1155 IF(lk_mpp) CALL mpp_sum(z_err)156 IF( z_err > 0 ) CALL ctl_stop( 'sbc_rnf : negative runnoff values exist' )157 !158 156 CALL iom_put( "runoffs", rnf ) ! output runoffs arrays 159 157 ENDIF … … 161 159 ENDIF 162 160 ! 161 ! ! ---------------------------------------- ! 163 162 IF( kt == nit000 ) THEN ! set the forcing field at nit000 - 1 ! 164 163 ! ! ---------------------------------------- ! … … 290 289 ENDIF 291 290 ! 291 IF( ln_rnf_emp .AND. nn_components == jp_iam_opa ) THEN 292 CALL ctl_stop( 'sbc_rnf_init: ln_rnf_emp must be false in case of SAS-OPA coupling' ) ; RETURN 293 ENDIF 292 294 ! ! ================== 293 295 ! ! Type of runoff … … 306 308 ELSE !== runoffs read in a file : set sf_rnf structure ==! 307 309 ! 308 ALLOCATE( sf_rnf(1), STAT=ierror ) ! Create sf_rnf structure (runoff inflow) 309 IF(lwp) WRITE(numout,*) 310 IF(lwp) WRITE(numout,*) ' runoffs inflow read in a file' 311 IF( ierror > 0 ) THEN 312 CALL ctl_stop( 'sbc_rnf: unable to allocate sf_rnf structure' ) ; RETURN 313 ENDIF 314 ALLOCATE( sf_rnf(1)%fnow(jpi,jpj,1) ) 315 IF( sn_rnf%ln_tint ) ALLOCATE( sf_rnf(1)%fdta(jpi,jpj,1,2) ) 316 ! ! fill sf_rnf with the namelist (sn_rnf) and control print 317 CALL fld_fill( sf_rnf, (/ sn_rnf /), cn_dir, 'sbc_rnf_init', 'read runoffs data', 'namsbc_rnf' ) 310 IF( .NOT. l_rnfcpl ) THEN 311 ALLOCATE( sf_rnf(1), STAT=ierror ) ! Create (if required) sf_rnf structure (runoff inflow) 312 IF(lwp) WRITE(numout,*) 313 IF(lwp) WRITE(numout,*) ' runoffs inflow read in a file' 314 IF( ierror > 0 ) THEN 315 CALL ctl_stop( 'sbc_rnf: unable to allocate sf_rnf structure' ) ; RETURN 316 ENDIF 317 ALLOCATE( sf_rnf(1)%fnow(jpi,jpj,1) ) 318 IF( sn_rnf%ln_tint ) ALLOCATE( sf_rnf(1)%fdta(jpi,jpj,1,2) ) 319 CALL fld_fill( sf_rnf, (/ sn_rnf /), cn_dir, 'sbc_rnf_init', 'read runoffs data', 'namsbc_rnf' ) 320 ENDIF 318 321 ! 319 322 IF( ln_rnf_tem ) THEN ! Create (if required) sf_t_rnf structure -
branches/2015/dev_r5218_CNRS17_coupling/NEMOGCM/NEMO/SAS_SRC/sbcssm.F90
r5376 r5381 180 180 WRITE(numout,*) 181 181 ENDIF 182 183 IF( nn_components /= jp_iam_sas ) THEN 184 ! 185 !! switch off stuff that isn't sensible with a standalone module 186 !! note that we need sbc_ssm called first in sbc 187 ! 188 IF( ln_apr_dyn ) THEN 189 IF( lwp ) WRITE(numout,*) 'No atmospheric gradient needed with StandAlone Surface scheme' 190 ln_apr_dyn = .FALSE. 191 ENDIF 192 IF( ln_dm2dc ) THEN 193 IF( lwp ) WRITE(numout,*) 'No diurnal cycle needed with StandAlone Surface scheme' 194 ln_dm2dc = .FALSE. 195 ENDIF 196 IF( ln_rnf ) THEN 197 IF( lwp ) WRITE(numout,*) 'No runoff needed with StandAlone Surface scheme' 198 ln_rnf = .FALSE. 199 ENDIF 200 IF( ln_ssr ) THEN 201 IF( lwp ) WRITE(numout,*) 'No surface relaxation needed with StandAlone Surface scheme' 202 ln_ssr = .FALSE. 203 ENDIF 204 IF( nn_fwb > 0 ) THEN 205 IF( lwp ) WRITE(numout,*) 'No freshwater budget adjustment needed with StandAlone Surface scheme' 206 nn_fwb = 0 207 ENDIF 208 IF( nn_closea > 0 ) THEN 209 IF( lwp ) WRITE(numout,*) 'No closed seas adjustment needed with StandAlone Surface scheme' 210 nn_closea = 0 211 ENDIF 212 182 ! 183 !! switch off stuff that isn't sensible with a standalone module 184 !! note that we need sbc_ssm called first in sbc 185 ! 186 IF( ln_apr_dyn ) THEN 187 IF( lwp ) WRITE(numout,*) 'No atmospheric gradient needed with StandAlone Surface scheme' 188 ln_apr_dyn = .FALSE. 189 ENDIF 190 IF( ln_rnf ) THEN 191 IF( lwp ) WRITE(numout,*) 'No runoff needed with StandAlone Surface scheme' 192 ln_rnf = .FALSE. 193 ENDIF 194 IF( ln_ssr ) THEN 195 IF( lwp ) WRITE(numout,*) 'No surface relaxation needed with StandAlone Surface scheme' 196 ln_ssr = .FALSE. 197 ENDIF 198 IF( nn_fwb > 0 ) THEN 199 IF( lwp ) WRITE(numout,*) 'No freshwater budget adjustment needed with StandAlone Surface scheme' 200 nn_fwb = 0 201 ENDIF 202 IF( nn_closea > 0 ) THEN 203 IF( lwp ) WRITE(numout,*) 'No closed seas adjustment needed with StandAlone Surface scheme' 204 nn_closea = 0 213 205 ENDIF 214 206 !
Note: See TracChangeset
for help on using the changeset viewer.
