Changeset 7037 for branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90

Timestamp:

2016-10-18T15:32:04+02:00 (8 years ago)

Author:

mocavero

Message:

ORCA2_LIM_PISCES hybrid version update

File:

• branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90 (modified) (5 diffs)

branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90

@@ -131 +131 @@
          SELECT CASE( cp_ice_msh )
          CASE( 'I' )                  !== B-grid ice dynamics :   I-point (i.e. F-point with sea-ice indexation)
+!$OMP PARALLEL DO schedule(static) private(jj, ji)
             DO jj = 2, jpj
                DO ji = 2, jpi   ! NO vector opt. possible
@@ -143 +144 @@
             !
          CASE( 'C' )                  !== C-grid ice dynamics :   U & V-points (same as ocean)
+!$OMP PARALLEL WORKSHARE
             u_oce(:,:) = ssu_m(:,:) * umask(:,:,1)                     ! mean surface ocean current at ice velocity point
             v_oce(:,:) = ssv_m(:,:) * vmask(:,:,1)
+!$OMP END PARALLEL WORKSHARE
             !
          END SELECT
@@ -150 +153 @@
          ! ... masked sea surface freezing temperature [Kelvin] (set to rt0 over land)
          CALL eos_fzp( sss_m(:,:), tfu(:,:) )
+!$OMP PARALLEL WORKSHARE
          tfu(:,:) = tfu(:,:) + rt0
 
          zsist (:,:,1) = sist (:,:) + rt0 * ( 1. - tmask(:,:,1) )
+!$OMP END PARALLEL WORKSHARE
 
          ! Ice albedo
@@ -164 +169 @@
 
             ! albedo depends on cloud fraction because of non-linear spectral effects
+!$OMP PARALLEL WORKSHARE
             zalb_ice(:,:,:) = ( 1. - cldf_ice ) * zalb_cs(:,:,:) + cldf_ice * zalb_os(:,:,:)
+!$OMP END PARALLEL WORKSHARE
             ! In CLIO the cloud fraction is read in the climatology and the all-sky albedo 
             ! (zalb_ice) is computed within the bulk routine
@@ -203 +210 @@
          IF( ln_mixcpl) THEN
             CALL sbc_cpl_ice_tau( zutau_ice , zvtau_ice )
+!$OMP PARALLEL WORKSHARE
             utau_ice(:,:) = utau_ice(:,:) * xcplmask(:,:,0) + zutau_ice(:,:) * ( 1. - xcplmask(:,:,0) )
             vtau_ice(:,:) = vtau_ice(:,:) * xcplmask(:,:,0) + zvtau_ice(:,:) * ( 1. - xcplmask(:,:,0) )
+!$OMP END PARALLEL WORKSHARE
          ENDIF