Changeset 10803
- Timestamp:
- 2019-03-26T11:46:36+01:00 (5 years ago)
- Location:
- NEMO/branches/UKMO/r8395_cpl-pressure/NEMOGCM/NEMO/OPA_SRC
- Files:
-
- 5 edited
-
SBC/sbcapr.F90 (modified) (4 diffs)
-
SBC/sbccpl.F90 (modified) (4 diffs)
-
SBC/sbcmod.F90 (modified) (3 diffs)
-
SBC/sbcssm.F90 (modified) (8 diffs)
-
step.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/r8395_cpl-pressure/NEMOGCM/NEMO/OPA_SRC/SBC/sbcapr.F90
r10797 r10803 26 26 27 27 ! !!* namsbc_apr namelist (Atmospheric PRessure) * 28 LOGICAL, PUBLIC :: cpl_mslp = .FALSE. ! Presure is passed via coupling 28 29 LOGICAL, PUBLIC :: ln_apr_obc !: inverse barometer added to OBC ssh data 29 30 LOGICAL, PUBLIC :: ln_ref_apr !: ref. pressure: global mean Patm (F) or a constant (F) 30 REAL(wp) 31 REAL(wp),PUBLIC :: rn_pref ! reference atmospheric pressure [N/m2] 31 32 32 33 REAL(wp), ALLOCATABLE, SAVE, PUBLIC, DIMENSION(:,:) :: ssh_ib ! Inverse barometer now sea surface height [m] … … 34 35 REAL(wp), ALLOCATABLE, SAVE, PUBLIC, DIMENSION(:,:) :: apr ! atmospheric pressure at kt [N/m2] 35 36 36 REAL(wp) :: tarea! whole domain mean masked ocean surface37 REAL(wp) :: r1_grau! = 1.e0 / (grav * rau0)37 REAL(wp), PUBLIC :: tarea ! whole domain mean masked ocean surface 38 REAL(wp), PUBLIC :: r1_grau ! = 1.e0 / (grav * rau0) 38 39 39 40 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_apr ! structure of input fields (file informations, fields read) … … 82 83 IF(lwm) WRITE ( numond, namsbc_apr ) 83 84 ! 84 ALLOCATE( sf_apr(1), STAT=ierror ) !* allocate and fill sf_sst (forcing structure) with sn_sst 85 IF( ierror > 0 ) CALL ctl_stop( 'STOP', 'sbc_apr: unable to allocate sf_apr structure' ) 86 ! 87 CALL fld_fill( sf_apr, (/ sn_apr /), cn_dir, 'sbc_apr', 'Atmospheric pressure ', 'namsbc_apr' ) 88 ALLOCATE( sf_apr(1)%fnow(jpi,jpj,1) ) 89 IF( sn_apr%ln_tint ) ALLOCATE( sf_apr(1)%fdta(jpi,jpj,1,2) ) 90 ALLOCATE( ssh_ib(jpi,jpj) , ssh_ibb(jpi,jpj) ) 91 ALLOCATE( apr (jpi,jpj) ) 85 IF( .NOT. cpl_mslp ) THEN 86 ALLOCATE( sf_apr(1), STAT=ierror ) !* allocate and fill sf_sst (forcing structure) with sn_sst 87 IF( ierror > 0 ) CALL ctl_stop( 'STOP', 'sbc_apr: unable to allocate sf_apr structure' ) 88 ! 89 CALL fld_fill( sf_apr, (/ sn_apr /), cn_dir, 'sbc_apr', 'Atmospheric pressure ', 'namsbc_apr' ) 90 ALLOCATE( sf_apr(1)%fnow(jpi,jpj,1) ) 91 IF( sn_apr%ln_tint ) ALLOCATE( sf_apr(1)%fdta(jpi,jpj,1,2) ) 92 ENDIF 93 ALLOCATE( ssh_ib(jpi,jpj) , ssh_ibb(jpi,jpj) ) 94 ALLOCATE( apr (jpi,jpj) ) 92 95 ! 93 96 IF(lwp) THEN !* control print 94 97 WRITE(numout,*) 95 WRITE(numout,*) ' Namelist namsbc_apr : Atmospheric PRessure as extrenal forcing' 98 IF( cpl_mslp ) THEN 99 WRITE(numout,*) ' Namelist namsbc_apr : Atmospheric Pressure as extrenal coupling' 100 ELSE 101 WRITE(numout,*) ' Namelist namsbc_apr : Atmospheric Pressure as extrenal forcing' 102 ENDIF 96 103 WRITE(numout,*) ' ref. pressure: global mean Patm (T) or a constant (F) ln_ref_apr = ', ln_ref_apr 97 104 ENDIF … … 115 122 ENDIF 116 123 117 ! ! ========================== ! 118 IF( MOD( kt-1, nn_fsbc ) == 0 ) THEN ! At each sbc time-step ! 119 ! ! ===========+++============ ! 120 ! 121 IF( kt /= nit000 ) ssh_ibb(:,:) = ssh_ib(:,:) !* Swap of ssh_ib fields 122 ! 123 CALL fld_read( kt, nn_fsbc, sf_apr ) !* input Patm provided at kt + nn_fsbc/2 124 ! 125 ! !* update the reference atmospheric pressure (if necessary) 126 IF( ln_ref_apr ) rn_pref = glob_sum( sf_apr(1)%fnow(:,:,1) * e1e2t(:,:) ) / tarea 127 ! 128 ! !* Patm related forcing at kt 129 ssh_ib(:,:) = - ( sf_apr(1)%fnow(:,:,1) - rn_pref ) * r1_grau ! equivalent ssh (inverse barometer) 130 apr (:,:) = sf_apr(1)%fnow(:,:,1) ! atmospheric pressure 131 ! 132 CALL iom_put( "ssh_ib", ssh_ib ) !* output the inverse barometer ssh 133 ENDIF 134 135 ! ! ---------------------------------------- ! 136 IF( kt == nit000 ) THEN ! set the forcing field at nit000 - 1 ! 137 ! ! ---------------------------------------- ! 138 ! !* Restart: read in restart file 139 IF( ln_rstart .AND. iom_varid( numror, 'ssh_ibb', ldstop = .FALSE. ) > 0 ) THEN 140 IF(lwp) WRITE(numout,*) 'sbc_apr: ssh_ibb read in the restart file' 141 CALL iom_get( numror, jpdom_autoglo, 'ssh_ibb', ssh_ibb ) ! before inv. barometer ssh 124 IF( .NOT. cpl_mslp ) THEN 125 ! ! ========================== ! 126 IF( MOD( kt-1, nn_fsbc ) == 0 ) THEN ! At each sbc time-step ! 127 ! ! ===========+++============ ! 142 128 ! 143 ELSE !* no restart: set from nit000 values 144 IF(lwp) WRITE(numout,*) 'sbc_apr: ssh_ibb set to nit000 values' 145 ssh_ibb(:,:) = ssh_ib(:,:) 129 IF( kt /= nit000 ) ssh_ibb(:,:) = ssh_ib(:,:) !* Swap of ssh_ib fields 130 ! 131 CALL fld_read( kt, nn_fsbc, sf_apr ) !* input Patm provided at kt + nn_fsbc/2 132 ! 133 ! !* update the reference atmospheric pressure (if necessary) 134 IF( ln_ref_apr ) rn_pref = glob_sum( sf_apr(1)%fnow(:,:,1) * e1e2t(:,:) ) / tarea 135 ! 136 ! !* Patm related forcing at kt 137 ssh_ib(:,:) = - ( sf_apr(1)%fnow(:,:,1) - rn_pref ) * r1_grau ! equivalent ssh (inverse barometer) 138 apr (:,:) = sf_apr(1)%fnow(:,:,1) ! atmospheric pressure 139 ! 140 CALL iom_put( "ssh_ib", ssh_ib ) !* output the inverse barometer ssh 146 141 ENDIF 147 ENDIF 148 ! ! ---------------------------------------- ! 149 IF( lrst_oce ) THEN ! Write in the ocean restart file ! 150 ! ! ---------------------------------------- ! 151 IF(lwp) WRITE(numout,*) 152 IF(lwp) WRITE(numout,*) 'sbc_apr : ssh_ib written in ocean restart file at it= ', kt,' date= ', ndastp 153 IF(lwp) WRITE(numout,*) '~~~~' 154 CALL iom_rstput( kt, nitrst, numrow, 'ssh_ibb' , ssh_ib ) 142 143 ! ! ---------------------------------------- ! 144 IF( kt == nit000 ) THEN ! set the forcing field at nit000 - 1 ! 145 ! ! ---------------------------------------- ! 146 ! !* Restart: read in restart file 147 IF( ln_rstart .AND. iom_varid( numror, 'ssh_ibb', ldstop = .FALSE. ) > 0 ) THEN 148 IF(lwp) WRITE(numout,*) 'sbc_apr: ssh_ibb read in the restart file' 149 CALL iom_get( numror, jpdom_autoglo, 'ssh_ibb', ssh_ibb ) ! before inv. barometer ssh 150 ! 151 ELSE !* no restart: set from nit000 values 152 IF(lwp) WRITE(numout,*) 'sbc_apr: ssh_ibb set to nit000 values' 153 ssh_ibb(:,:) = ssh_ib(:,:) 154 ENDIF 155 ENDIF 156 ! ! ---------------------------------------- ! 157 IF( lrst_oce ) THEN ! Write in the ocean restart file ! 158 ! ! ---------------------------------------- ! 159 IF(lwp) WRITE(numout,*) 160 IF(lwp) WRITE(numout,*) 'sbc_apr : ssh_ib written in ocean restart file at it= ', kt,' date= ', ndastp 161 IF(lwp) WRITE(numout,*) '~~~~' 162 CALL iom_rstput( kt, nitrst, numrow, 'ssh_ibb' , ssh_ib ) 163 ENDIF 155 164 ENDIF 156 165 ! -
NEMO/branches/UKMO/r8395_cpl-pressure/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90
r10797 r10803 182 182 183 183 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: albedo_oce_mix ! ocean albedo sent to atmosphere (mix clear/overcast sky) 184 185 REAL(wp) :: rpref = 101000._wp ! reference atmospheric pressure[N/m2]186 REAL(wp) :: r1_grau ! = 1.e0 / (grav * rau0)187 184 188 185 INTEGER , ALLOCATABLE, SAVE, DIMENSION( :) :: nrcvinfo ! OASIS info argument … … 535 532 ! ! Mean Sea Level Pressure ! 536 533 ! ! ------------------------- ! 537 srcv(jpr_mslp)%clname = 'O_MSLP' ; IF( TRIM(sn_rcv_mslp%cldes ) == 'coupled' ) srcv(jpr_mslp)%laction = .TRUE. 538 534 srcv(jpr_mslp)%clname = 'O_MSLP' 535 IF( TRIM(sn_rcv_mslp%cldes ) == 'coupled' ) THEN 536 srcv(jpr_mslp)%laction = .TRUE. 537 cpl_mslp = .TRUE. 538 ENDIF 539 539 ! ! ------------------------- ! 540 540 ! ! topmelt and botmelt ! … … 975 975 !!---------------------------------------------------------------------- 976 976 USE zdf_oce, ONLY : ln_zdfqiao 977 USE sbcssm , ONLY : sbc_ssm_cpl 978 USE lib_fortran ! distributed memory computing library 977 979 978 980 IMPLICIT NONE … … 1138 1140 IF( kt /= nit000 ) ssh_ibb(:,:) = ssh_ib(:,:) !* Swap of ssh_ib fields 1139 1141 1140 r1_grau = 1.e0 / (grav * rau0) !* constant for optimization 1141 ssh_ib(:,:) = - ( frcv(jpr_mslp)%z3(:,:,1) - rpref ) * r1_grau ! equivalent ssh (inverse barometer) 1142 apr (:,:) = frcv(jpr_mslp)%z3(:,:,1) !atmospheric pressure 1143 1144 IF( kt == nit000 ) ssh_ibb(:,:) = ssh_ib(:,:) ! correct this later (read from restart if possible) 1142 ! !* update the reference atmospheric pressure (if necessary) 1143 IF( ln_ref_apr ) rn_pref = glob_sum( frcv(jpr_mslp)%z3(:,:,1) * e1e2t(:,:) ) / tarea 1144 1145 ssh_ib(:,:) = - ( frcv(jpr_mslp)%z3(:,:,1) - rn_pref ) * r1_grau ! equivalent ssh (inverse barometer) 1146 apr (:,:) = frcv(jpr_mslp)%z3(:,:,1) !atmospheric pressure 1147 ! 1148 CALL iom_put( "ssh_ib", ssh_ib ) !* output the inverse barometer ssh 1149 ! ! ---------------------------------------- ! 1150 IF( kt == nit000 ) THEN ! set the forcing field at nit000 - 1 ! 1151 ! ! ---------------------------------------- ! 1152 !* Restart: read in restart file 1153 IF( ln_rstart .AND. iom_varid( numror, 'ssh_ibb', ldstop = .FALSE. ) > 0 ) THEN 1154 IF(lwp) WRITE(numout,*) 'sbc_cpl: ssh_ibb read in the restart file' 1155 CALL iom_get( numror, jpdom_autoglo, 'ssh_ibb', ssh_ibb ) ! before inv. barometer ssh 1156 ELSE !* no restart: set from nit000 values 1157 IF(lwp) WRITE(numout,*) 'sbc_cpl: ssh_ibb set to nit000 values' 1158 ssh_ibb(:,:) = ssh_ib(:,:) 1159 ENDIF 1160 ENDIF 1161 ! ! ---------------------------------------- ! 1162 IF( lrst_oce ) THEN ! Write in the ocean restart file ! 1163 ! ! ---------------------------------------- ! 1164 IF(lwp) WRITE(numout,*) 1165 IF(lwp) WRITE(numout,*) 'sbc_cpl : ssh_ib written in ocean restart file at it= ', kt,' date= ', ndastp 1166 IF(lwp) WRITE(numout,*) '~~~~' 1167 CALL iom_rstput( kt, nitrst, numrow, 'ssh_ibb' , ssh_ib ) 1168 ENDIF 1169 1170 ! Update mean ssh 1171 CALL sbc_ssm_cpl( kt ) 1145 1172 END IF 1146 1173 ! -
NEMO/branches/UKMO/r8395_cpl-pressure/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90
r10797 r10803 357 357 !! - updte the ice fraction : fr_i 358 358 !!---------------------------------------------------------------------- 359 USE sbcapr, ONLY: sbc_apr 360 USE bdydta, ONLY: bdy_dta 361 ! 359 362 INTEGER, INTENT(in) :: kt ! ocean time step 360 363 ! … … 385 388 ! ! forcing field computation ! 386 389 ! ! ---------------------------------------- ! 390 ! 391 IF( ln_apr_dyn ) CALL sbc_apr( kt ) ! atmospheric pressure provided at kt+0.5*nn_fsbc 392 ! (caution called before sbc_ssm) 387 393 ! 388 394 ll_sas = nn_components == jp_iam_sas ! component flags … … 420 426 END IF 421 427 IF( ln_mixcpl ) CALL sbc_cpl_rcv ( kt, nn_fsbc, nn_ice ) ! forced-coupled mixed formulation after forcing 422 428 IF( ln_bdy ) CALL bdy_dta ( kt, time_offset=+1 ) ! update dynamic & tracer data at open boundaries 429 ! (caution called after sbc_ssm[_cpl] and before ice) 423 430 ! 424 431 ! !== Misc. Options ==! -
NEMO/branches/UKMO/r8395_cpl-pressure/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssm.F90
r10797 r10803 27 27 28 28 PUBLIC sbc_ssm ! routine called by step.F90 29 PUBLIC sbc_ssm_cpl ! routine called by sbccpl.F90 29 30 PUBLIC sbc_ssm_init ! routine called by sbcmod.F90 30 31 … … 75 76 sss_m(:,:) = zts(:,:,jp_sal) 76 77 ! ! removed inverse barometer ssh when Patm forcing is used (for sea-ice dynamics) 77 IF( ln_apr_dyn ) THEN ; ssh_m(:,:) = sshn(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) ) 78 ELSE ; ssh_m(:,:) = sshn(:,:) 78 IF( .NOT. cpl_mslp ) THEN 79 IF( ln_apr_dyn ) THEN ; ssh_m(:,:) = sshn(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) ) 80 ELSE ; ssh_m(:,:) = sshn(:,:) 81 ENDIF 79 82 ENDIF 80 83 ! … … 98 101 sss_m(:,:) = zcoef * zts(:,:,jp_sal) 99 102 ! ! removed inverse barometer ssh when Patm forcing is used (for sea-ice dynamics) 100 IF( ln_apr_dyn ) THEN ; ssh_m(:,:) = zcoef * ( sshn(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) ) ) 101 ELSE ; ssh_m(:,:) = zcoef * sshn(:,:) 103 IF( .NOT. cpl_mslp ) THEN 104 IF( ln_apr_dyn ) THEN ; ssh_m(:,:) = zcoef * ( sshn(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) ) ) 105 ELSE ; ssh_m(:,:) = zcoef * sshn(:,:) 106 ENDIF 102 107 ENDIF 103 108 ! … … 112 117 sst_m(:,:) = 0._wp 113 118 sss_m(:,:) = 0._wp 114 ssh_m(:,:) = 0._wp119 IF( .NOT. cpl_mslp ) ssh_m(:,:) = 0._wp 115 120 e3t_m(:,:) = 0._wp 116 121 frq_m(:,:) = 0._wp … … 126 131 sss_m(:,:) = sss_m(:,:) + zts(:,:,jp_sal) 127 132 ! ! removed inverse barometer ssh when Patm forcing is used (for sea-ice dynamics) 128 IF( ln_apr_dyn ) THEN ; ssh_m(:,:) = ssh_m(:,:) + sshn(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) ) 129 ELSE ; ssh_m(:,:) = ssh_m(:,:) + sshn(:,:) 133 IF( .NOT. cpl_mslp ) THEN 134 IF( ln_apr_dyn ) THEN ; ssh_m(:,:) = ssh_m(:,:) + sshn(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) ) 135 ELSE ; ssh_m(:,:) = ssh_m(:,:) + sshn(:,:) 136 ENDIF 130 137 ENDIF 131 138 ! … … 142 149 ssu_m(:,:) = ssu_m(:,:) * zcoef ! mean suface current [m/s] 143 150 ssv_m(:,:) = ssv_m(:,:) * zcoef ! 144 ssh_m(:,:) = ssh_m(:,:) * zcoef ! mean SSH [m]151 IF( .NOT. cpl_mslp ) ssh_m(:,:) = ssh_m(:,:) * zcoef ! mean SSH [m] 145 152 e3t_m(:,:) = e3t_m(:,:) * zcoef ! mean vertical scale factor [m] 146 153 frq_m(:,:) = frq_m(:,:) * zcoef ! mean fraction of solar net radiation absorbed in the 1st T level [-] … … 160 167 CALL iom_rstput( kt, nitrst, numrow, 'sst_m' , sst_m ) 161 168 CALL iom_rstput( kt, nitrst, numrow, 'sss_m' , sss_m ) 162 CALL iom_rstput( kt, nitrst, numrow, 'ssh_m' , ssh_m )169 IF( .NOT. cpl_mslp ) CALL iom_rstput( kt, nitrst, numrow, 'ssh_m' , ssh_m ) 163 170 CALL iom_rstput( kt, nitrst, numrow, 'e3t_m' , e3t_m ) 164 171 CALL iom_rstput( kt, nitrst, numrow, 'frq_m' , frq_m ) … … 180 187 END SUBROUTINE sbc_ssm 181 188 189 SUBROUTINE sbc_ssm_cpl( kt ) 190 !!--------------------------------------------------------------------- 191 !! *** ROUTINE sbc_ssm_cpl *** 192 !! 193 !! ** Purpose : provide ocean surface variable to sea-surface boundary 194 !! condition computation when pressure is read from coupling 195 !! 196 !! ** Method : The inverse barometer ssh (i.e. ssh associated with Patm) 197 !! is added to ssh_m when ln_apr_dyn = T. Required for sea-ice dynamics. 198 !!--------------------------------------------------------------------- 199 INTEGER, INTENT(in) :: kt ! ocean time step 200 ! 201 REAL(wp) :: zcoef ! local scalar 202 !!--------------------------------------------------------------------- 203 ! 204 IF( nn_fsbc == 1 ) THEN ! Instantaneous surface fields ! 205 ! ! ---------------------------------------- ! 206 IF( ln_apr_dyn ) THEN ; ssh_m(:,:) = sshn(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) ) 207 ELSE ; ssh_m(:,:) = sshn(:,:) 208 ENDIF 209 ELSE 210 ! ! ----------------------------------------------- ! 211 IF( kt == nit000 .AND. .NOT. l_ssm_mean ) THEN ! Initialisation: 1st time-step, no input means ! 212 ! ! ----------------------------------------------- ! 213 IF(lwp) WRITE(numout,*) 214 IF(lwp) WRITE(numout,*) '~~~~~~~ mean ssh field initialised to instantaneous values' 215 zcoef = REAL( nn_fsbc - 1, wp ) 216 zcoef = REAL( nn_fsbc - 1, wp ) 217 IF( ln_apr_dyn ) THEN ; ssh_m(:,:) = zcoef * ( sshn(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) ) ) 218 ELSE ; ssh_m(:,:) = zcoef * sshn(:,:) 219 ENDIF 220 ! ! ---------------------------------------- ! 221 ELSEIF( MOD( kt - 2 , nn_fsbc ) == 0 ) THEN ! Initialisation: New mean computation ! 222 ! ! ---------------------------------------- ! 223 ssh_m(:,:) = 0.e0 224 ENDIF 225 226 IF( ln_apr_dyn ) THEN ; ssh_m(:,:) = ssh_m(:,:) + sshn(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) ) 227 ELSE ; ssh_m(:,:) = ssh_m(:,:) + sshn(:,:) 228 ENDIF 229 ! ! ---------------------------------------- ! 230 IF( MOD( kt - 1 , nn_fsbc ) == 0 ) THEN ! Mean value at each nn_fsbc time-step ! 231 ! ! ---------------------------------------- ! 232 zcoef = 1. / REAL( nn_fsbc, wp ) 233 ssh_m(:,:) = ssh_m(:,:) * zcoef ! mean SSH [m] 234 ENDIF 235 ! ! ---------------------------------------- ! 236 IF( lrst_oce ) THEN ! Write in the ocean restart file ! 237 ! ! ---------------------------------------- ! 238 IF(lwp) WRITE(numout,*) 239 IF(lwp) WRITE(numout,*) 'sbc_ssm_cpl : ssh mean field written in ocean restart file ', & 240 & 'at it= ', kt,' date= ', ndastp 241 IF(lwp) WRITE(numout,*) '~~~~~~~' 242 CALL iom_rstput( kt, nitrst, numrow, 'ssh_m' , ssh_m ) 243 ENDIF 244 ENDIF 245 ! 246 IF( MOD( kt - 1 , nn_fsbc ) == 0 ) THEN ! Mean value at each nn_fsbc time-step ! 247 CALL iom_put( 'ssh_m', ssh_m ) 248 ENDIF 249 ! 250 END SUBROUTINE sbc_ssm_cpl 182 251 183 252 SUBROUTINE sbc_ssm_init -
NEMO/branches/UKMO/r8395_cpl-pressure/NEMOGCM/NEMO/OPA_SRC/step.F90
r10797 r10803 106 106 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 107 107 IF( ln_tide ) CALL sbc_tide( kstp ) ! update tide potential 108 IF( ln_apr_dyn ) CALL sbc_apr ( kstp ) ! atmospheric pressure (NB: call before bdy_dta which needs ssh_ib)109 IF( ln_bdy ) CALL bdy_dta ( kstp, time_offset=+1 ) ! update dynamic & tracer data at open boundaries110 108 CALL sbc ( kstp ) ! Sea Boundary Condition (including sea-ice) 111 109
Note: See TracChangeset
for help on using the changeset viewer.
