Changeset 5961

Timestamp:

2015-12-01T12:37:10+01:00 (9 years ago)

Author:

davestorkey

Message:

UKMO/antarctic_partial_slip branch: commit code change.

File:

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

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

@@ -136 +136 @@
       INTEGER  ::   isrow                    ! index for ORCA1 starting row
       INTEGER , POINTER, DIMENSION(:,:) ::  imsk
+      REAL(wp) ::  zphi_drake_passage, zshlat_antarc
       REAL(wp), POINTER, DIMENSION(:,:) ::  zwf
       !!
@@ -445 +446 @@
       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