Changeset 6748 for branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/SBC
- Timestamp:
- 2016-06-28T11:53:56+02:00 (8 years ago)
- Location:
- branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/SBC
- Files:
-
- 3 edited
-
sbcana.F90 (modified) (4 diffs)
-
sbcmod.F90 (modified) (6 diffs)
-
sbcssm.F90 (modified) (9 diffs)
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/SBC/sbcana.F90
r6140 r6748 203 203 ztrp= - 40.e0 ! retroaction term on heat fluxes (W/m2/K) 204 204 zconv = 3.16e-5 ! convertion factor: 1 m/yr => 3.16e-5 mm/s 205 !$OMP PARALLEL DO schedule(static) private(jj, ji, t_star) 205 206 DO jj = 1, jpj 206 207 DO ji = 1, jpi … … 238 239 239 240 ! freshwater (mass flux) and update of qns with heat content of emp 241 !$OMP PARALLEL WORKSHARE 240 242 emp (:,:) = emp(:,:) - zsumemp * tmask(:,:,1) ! freshwater flux (=0 in domain average) 241 243 sfx (:,:) = 0.0_wp ! no salt flux 242 244 qns (:,:) = qns(:,:) - emp(:,:) * sst_m(:,:) * rcp ! evap and precip are at SST 243 245 !$OMP END PARALLEL WORKSHARE 244 246 245 247 ! ---------------------------- ! … … 267 269 ztau_sais = 0.015 268 270 ztaun = ztau - ztau_sais * COS( (ztime - ztimemax) / (ztimemin - ztimemax) * rpi ) 271 !$OMP PARALLEL DO schedule(static) private(jj, ji) 269 272 DO jj = 1, jpj 270 273 DO ji = 1, jpi … … 278 281 ! module of wind stress and wind speed at T-point 279 282 zcoef = 1. / ( zrhoa * zcdrag ) 283 !$OMP PARALLEL DO schedule(static) private(jj, ji, ztx, zty, zmod) 280 284 DO jj = 2, jpjm1 281 285 DO ji = fs_2, fs_jpim1 ! vect. opt. -
branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90
r6460 r6748 187 187 IF( .NOT. ln_isf ) THEN ! variable initialisation if no ice shelf 188 188 IF( sbc_isf_alloc() /= 0 ) CALL ctl_stop( 'STOP', 'sbc_init : unable to allocate sbc_isf arrays' ) 189 !$OMP PARALLEL WORKSHARE 189 190 fwfisf (:,:) = 0.0_wp ; fwfisf_b (:,:) = 0.0_wp 190 191 risf_tsc(:,:,:) = 0.0_wp ; risf_tsc_b(:,:,:) = 0.0_wp 192 !$OMP END PARALLEL WORKSHARE 191 193 END IF 192 194 IF( nn_ice == 0 .AND. nn_components /= jp_iam_opa ) fr_i(:,:) = 0._wp ! no ice in the domain, ice fraction is always zero 193 195 !$OMP PARALLEL WORKSHARE 194 196 sfx(:,:) = 0._wp ! the salt flux due to freezing/melting will be computed (i.e. will be non-zero) 195 197 ! only if sea-ice is present … … 198 200 199 201 taum(:,:) = 0._wp ! Initialise taum for use in gls in case of reduced restart 200 202 !$OMP END PARALLEL WORKSHARE 201 203 ! ! restartability 202 204 IF( ( nn_ice == 2 .OR. nn_ice ==3 ) .AND. .NOT.( ln_blk_clio .OR. ln_blk_core .OR. ln_cpl ) ) & … … 318 320 IF( kt /= nit000 ) THEN ! Swap of forcing fields ! 319 321 ! ! ---------------------------------------- ! 322 !$OMP PARALLEL WORKSHARE 320 323 utau_b(:,:) = utau(:,:) ! Swap the ocean forcing fields 321 324 vtau_b(:,:) = vtau(:,:) ! (except at nit000 where before fields … … 323 326 emp_b (:,:) = emp (:,:) 324 327 sfx_b (:,:) = sfx (:,:) 328 !$OMP END PARALLEL WORKSHARE 325 329 IF ( ln_rnf ) THEN 330 !$OMP PARALLEL WORKSHARE 326 331 rnf_b (:,: ) = rnf (:,: ) 327 332 rnf_tsc_b(:,:,:) = rnf_tsc(:,:,:) 333 !$OMP END PARALLEL WORKSHARE 328 334 ENDIF 329 335 ENDIF … … 404 410 CALL iom_get( numror, jpdom_autoglo, 'sfx_b', sfx_b ) ! before salt flux (T-point) 405 411 ELSE 412 !$OMP PARALLEL WORKSHARE 406 413 sfx_b (:,:) = sfx(:,:) 414 !$OMP END PARALLEL WORKSHARE 407 415 ENDIF 408 416 ELSE !* no restart: set from nit000 values 409 417 IF(lwp) WRITE(numout,*) ' nit000-1 surface forcing fields set to nit000' 418 !$OMP PARALLEL WORKSHARE 410 419 utau_b(:,:) = utau(:,:) 411 420 vtau_b(:,:) = vtau(:,:) … … 413 422 emp_b (:,:) = emp(:,:) 414 423 sfx_b (:,:) = sfx(:,:) 424 !$OMP END PARALLEL WORKSHARE 415 425 ENDIF 416 426 ENDIF -
branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssm.F90
r6489 r6748 59 59 ! 60 60 ! !* surface T-, U-, V- ocean level variables (T, S, depth, velocity) 61 !$OMP PARALLEL DO schedule(static) private(jj, ji) 61 62 DO jj = 1, jpj 62 63 DO ji = 1, jpi … … 68 69 IF( nn_fsbc == 1 ) THEN ! Instantaneous surface fields ! 69 70 ! ! ---------------------------------------- ! 71 !$OMP PARALLEL WORKSHARE 70 72 ssu_m(:,:) = ub(:,:,1) 71 73 ssv_m(:,:) = vb(:,:,1) 74 !$OMP END PARALLEL WORKSHARE 72 75 IF( l_useCT ) THEN ; sst_m(:,:) = eos_pt_from_ct( zts(:,:,jp_tem), zts(:,:,jp_sal) ) 73 76 ELSE ; sst_m(:,:) = zts(:,:,jp_tem) 74 77 ENDIF 78 !$OMP PARALLEL WORKSHARE 75 79 sss_m(:,:) = zts(:,:,jp_sal) 80 !$OMP END PARALLEL WORKSHARE 76 81 ! ! 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(:,:) 79 ENDIF 80 ! 82 IF( ln_apr_dyn ) THEN 83 !$OMP PARALLEL WORKSHARE 84 ssh_m(:,:) = sshn(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) ) 85 !$OMP END PARALLEL WORKSHARE 86 ELSE 87 !$OMP PARALLEL WORKSHARE 88 ssh_m(:,:) = sshn(:,:) 89 !$OMP END PARALLEL WORKSHARE 90 ENDIF 91 ! 92 !$OMP PARALLEL WORKSHARE 81 93 e3t_m(:,:) = e3t_n(:,:,1) 82 94 ! 83 95 frq_m(:,:) = fraqsr_1lev(:,:) 96 !$OMP END PARALLEL WORKSHARE 84 97 ! 85 98 ELSE … … 90 103 IF(lwp) WRITE(numout,*) '~~~~~~~ mean fields initialised to instantaneous values' 91 104 zcoef = REAL( nn_fsbc - 1, wp ) 105 !$OMP PARALLEL WORKSHARE 92 106 ssu_m(:,:) = zcoef * ub(:,:,1) 93 107 ssv_m(:,:) = zcoef * vb(:,:,1) 108 !$OMP END PARALLEL WORKSHARE 94 109 IF( l_useCT ) THEN ; sst_m(:,:) = zcoef * eos_pt_from_ct( zts(:,:,jp_tem), zts(:,:,jp_sal) ) 95 110 ELSE ; sst_m(:,:) = zcoef * zts(:,:,jp_tem) 96 111 ENDIF 112 !$OMP PARALLEL WORKSHARE 97 113 sss_m(:,:) = zcoef * zts(:,:,jp_sal) 114 !$OMP END PARALLEL WORKSHARE 98 115 ! ! removed inverse barometer ssh when Patm forcing is used (for sea-ice dynamics) 99 IF( ln_apr_dyn ) THEN ; ssh_m(:,:) = zcoef * ( sshn(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) ) ) 100 ELSE ; ssh_m(:,:) = zcoef * sshn(:,:) 116 IF( ln_apr_dyn ) THEN 117 !$OMP PARALLEL WORKSHARE 118 ssh_m(:,:) = zcoef * ( sshn(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) ) ) 119 !$OMP END PARALLEL WORKSHARE 120 ELSE 121 !$OMP PARALLEL WORKSHARE 122 ssh_m(:,:) = zcoef * sshn(:,:) 123 !$OMP END PARALLEL WORKSHARE 101 124 ENDIF 102 125 ! 126 !$OMP PARALLEL WORKSHARE 103 127 e3t_m(:,:) = zcoef * e3t_n(:,:,1) 104 128 ! 105 129 frq_m(:,:) = zcoef * fraqsr_1lev(:,:) 130 !$OMP END PARALLEL WORKSHARE 106 131 ! ! ---------------------------------------- ! 107 132 ELSEIF( MOD( kt - 2 , nn_fsbc ) == 0 ) THEN ! Initialisation: New mean computation ! 108 133 ! ! ---------------------------------------- ! 134 !$OMP PARALLEL WORKSHARE 109 135 ssu_m(:,:) = 0._wp ! reset to zero ocean mean sbc fields 110 136 ssv_m(:,:) = 0._wp … … 114 140 e3t_m(:,:) = 0._wp 115 141 frq_m(:,:) = 0._wp 142 !$OMP END PARALLEL WORKSHARE 116 143 ENDIF 117 144 ! ! ---------------------------------------- ! 118 145 ! ! Cumulate at each time step ! 119 146 ! ! ---------------------------------------- ! 147 !$OMP PARALLEL WORKSHARE 120 148 ssu_m(:,:) = ssu_m(:,:) + ub(:,:,1) 121 149 ssv_m(:,:) = ssv_m(:,:) + vb(:,:,1) 150 !$OMP END PARALLEL WORKSHARE 122 151 IF( l_useCT ) THEN ; sst_m(:,:) = sst_m(:,:) + eos_pt_from_ct( zts(:,:,jp_tem), zts(:,:,jp_sal) ) 123 152 ELSE ; sst_m(:,:) = sst_m(:,:) + zts(:,:,jp_tem) 124 153 ENDIF 154 !$OMP PARALLEL WORKSHARE 125 155 sss_m(:,:) = sss_m(:,:) + zts(:,:,jp_sal) 156 !$OMP END PARALLEL WORKSHARE 126 157 ! ! removed inverse barometer ssh when Patm forcing is used (for sea-ice dynamics) 127 IF( ln_apr_dyn ) THEN ; ssh_m(:,:) = ssh_m(:,:) + sshn(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) ) 128 ELSE ; ssh_m(:,:) = ssh_m(:,:) + sshn(:,:) 129 ENDIF 130 ! 158 IF( ln_apr_dyn ) THEN 159 !$OMP PARALLEL WORKSHARE 160 ssh_m(:,:) = ssh_m(:,:) + sshn(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) ) 161 !$OMP END PARALLEL WORKSHARE 162 ELSE 163 !$OMP PARALLEL WORKSHARE 164 ssh_m(:,:) = ssh_m(:,:) + sshn(:,:) 165 !$OMP END PARALLEL WORKSHARE 166 ENDIF 167 ! 168 !$OMP PARALLEL WORKSHARE 131 169 e3t_m(:,:) = e3t_m(:,:) + e3t_n(:,:,1) 132 170 ! 133 171 frq_m(:,:) = frq_m(:,:) + fraqsr_1lev(:,:) 172 !$OMP END PARALLEL WORKSHARE 134 173 135 174 ! ! ---------------------------------------- ! … … 137 176 ! ! ---------------------------------------- ! 138 177 zcoef = 1. / REAL( nn_fsbc, wp ) 178 !$OMP PARALLEL WORKSHARE 139 179 sst_m(:,:) = sst_m(:,:) * zcoef ! mean SST [Celcius] 140 180 sss_m(:,:) = sss_m(:,:) * zcoef ! mean SSS [psu] … … 144 184 e3t_m(:,:) = e3t_m(:,:) * zcoef ! mean vertical scale factor [m] 145 185 frq_m(:,:) = frq_m(:,:) * zcoef ! mean fraction of solar net radiation absorbed in the 1st T level [-] 186 !$OMP END PARALLEL WORKSHARE 146 187 ! 147 188 ENDIF … … 223 264 & 'from ', zf_sbc, ' to ', nn_fsbc 224 265 zcoef = REAL( nn_fsbc - 1, wp ) / zf_sbc 266 !$OMP PARALLEL WORKSHARE 225 267 ssu_m(:,:) = zcoef * ssu_m(:,:) 226 268 ssv_m(:,:) = zcoef * ssv_m(:,:) … … 230 272 e3t_m(:,:) = zcoef * e3t_m(:,:) 231 273 frq_m(:,:) = zcoef * frq_m(:,:) 274 !$OMP END PARALLEL WORKSHARE 232 275 ELSE 233 276 IF(lwp) WRITE(numout,*) '~~~~~~~ mean fields read in the ocean restart file' … … 239 282 ! 240 283 IF(lwp) WRITE(numout,*) ' default initialisation of ss?_m arrays' 284 !$OMP PARALLEL WORKSHARE 241 285 ssu_m(:,:) = ub(:,:,1) 242 286 ssv_m(:,:) = vb(:,:,1) 287 !$OMP END PARALLEL WORKSHARE 243 288 IF( l_useCT ) THEN ; sst_m(:,:) = eos_pt_from_ct( tsn(:,:,1,jp_tem), tsn(:,:,1,jp_sal) ) 244 289 ELSE ; sst_m(:,:) = tsn(:,:,1,jp_tem) 245 290 ENDIF 291 !$OMP PARALLEL WORKSHARE 246 292 sss_m(:,:) = tsn (:,:,1,jp_sal) 247 293 ssh_m(:,:) = sshn (:,:) 248 294 e3t_m(:,:) = e3t_n(:,:,1) 249 295 frq_m(:,:) = 1._wp 296 !$OMP END PARALLEL WORKSHARE 250 297 ! 251 298 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.
