Changeset 5961


Ignore:
Timestamp:
2015-12-01T12:37:10+01:00 (5 years ago)
Author:
davestorkey
Message:

UKMO/antarctic_partial_slip branch: commit code change.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/antarctic_partial_slip/NEMOGCM/NEMO/OPA_SRC/DOM/dommsk.F90

    r5958 r5961  
    136136      INTEGER  ::   isrow                    ! index for ORCA1 starting row 
    137137      INTEGER , POINTER, DIMENSION(:,:) ::  imsk 
     138      REAL(wp) ::  zphi_drake_passage, zshlat_antarc 
    138139      REAL(wp), POINTER, DIMENSION(:,:) ::  zwf 
    139140      !! 
     
    445446      ENDIF 
    446447      ! 
     448      IF( cp_cfg == "orca" .AND. jp_cfg == 025 .AND. rn_shlat == 0.0 ) THEN    
     449         !                                              ! ORCA_R025 configuration 
     450         !                                              ! Increased lateral friction on parts of Antarctic coastline 
     451         !                                              ! for increased stability 
     452         !                                              ! NB. This only works to do this here if we have free slip  
     453         !                                              ! generally, so fmask is zero at coast points. 
     454         IF(lwp) WRITE(numout,*) 
     455         IF(lwp) WRITE(numout,*) '   orca_r025: increase friction in following regions : ' 
     456         IF(lwp) WRITE(numout,*) '      whole Antarctic coastline: partial slip shlat=1 ' 
     457 
     458         zphi_drake_passage = -58.0_wp 
     459         zshlat_antarc = 1.0_wp 
     460         zwf(:,:) = fmask(:,:,1)          
     461         DO jj = 2, jpjm1 
     462            DO ji = fs_2, fs_jpim1   ! vector opt. 
     463               IF( gphif(ji,jj) .lt. zphi_drake_passage .and. fmask(ji,jj,1) == 0._wp ) THEN 
     464                  fmask(ji,jj,:) = zshlat_antarc * MIN( 1._wp , MAX( zwf(ji+1,jj), zwf(ji,jj+1),   & 
     465                     &                                           zwf(ji-1,jj), zwf(ji,jj-1)  )  ) 
     466               ENDIF 
     467            END DO 
     468         END DO 
     469      END IF 
     470      ! 
    447471      CALL lbc_lnk( fmask, 'F', 1._wp )      ! Lateral boundary conditions on fmask 
    448472 
Note: See TracChangeset for help on using the changeset viewer.