Changeset 6507 for branches/UKMO/dev_r6501_GO6_package_trunk/NEMOGCM/NEMO/OPA_SRC/DOM/dommsk.F90

Timestamp:

2016-05-03T14:28:12+02:00 (8 years ago)

Author:

timgraham

Message:

First attempt at merging in science changes from GO6 package branch at v3.6 stable (Note-namelists not yet dealt with)

File:

• branches/UKMO/dev_r6501_GO6_package_trunk/NEMOGCM/NEMO/OPA_SRC/DOM/dommsk.F90 (modified) (2 diffs)

branches/UKMO/dev_r6501_GO6_package_trunk/NEMOGCM/NEMO/OPA_SRC/DOM/dommsk.F90

@@ -105 +105 @@
       INTEGER  ::   isrow                    ! index for ORCA1 starting row
       INTEGER , POINTER, DIMENSION(:,:) ::  imsk
+      REAL(wp) ::  zphi_drake_passage, zshlat_antarc
       REAL(wp), POINTER, DIMENSION(:,:) ::  zwf
       !!
@@ -345 +346 @@
       ENDIF
       !
+      IF( cp_cfg == "orca" .AND. jp_cfg == 025 .AND. rn_shlat == 0.0 ) THEN   
+         !                                              ! ORCA_R025 configuration
+         !                                              ! Increased lateral friction on parts of Antarctic coastline
+         !                                              ! for increased stability
+         !                                              ! NB. This only works to do this here if we have free slip 
+         !                                              ! generally, so fmask is zero at coast points.
+         IF(lwp) WRITE(numout,*)
+         IF(lwp) WRITE(numout,*) '   orca_r025: increase friction in following regions : '
+         IF(lwp) WRITE(numout,*) '      whole Antarctic coastline: partial slip shlat=1 '
+
+         zphi_drake_passage = -58.0_wp
+         zshlat_antarc = 1.0_wp
+         zwf(:,:) = fmask(:,:,1)         
+         DO jj = 2, jpjm1
+            DO ji = fs_2, fs_jpim1   ! vector opt.
+               IF( gphif(ji,jj) .lt. zphi_drake_passage .and. fmask(ji,jj,1) == 0._wp ) THEN
+                  fmask(ji,jj,:) = zshlat_antarc * MIN( 1._wp , MAX( zwf(ji+1,jj), zwf(ji,jj+1),   &
+                     &                                           zwf(ji-1,jj), zwf(ji,jj-1)  )  )
+               ENDIF
+            END DO
+         END DO
+      END IF
+      !
       CALL lbc_lnk( fmask, 'F', 1._wp )      ! Lateral boundary conditions on fmask
       !