New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 10965 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/DIA/diaharm.F90 – NEMO

Ignore:
Timestamp:
2019-05-10T18:02:51+02:00 (5 years ago)
Author:
davestorkey
Message:

2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps : DIA and stpctl.F90. Just testing in ORCA1 so far.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/DIA/diaharm.F90

    r10425 r10965  
    162162 
    163163 
    164    SUBROUTINE dia_harm ( kt ) 
     164   SUBROUTINE dia_harm ( kt, Kmm ) 
    165165      !!---------------------------------------------------------------------- 
    166166      !!                 ***  ROUTINE dia_harm  *** 
     
    172172      !!-------------------------------------------------------------------- 
    173173      INTEGER, INTENT( IN ) :: kt 
     174      INTEGER, INTENT( IN ) :: Kmm     ! time level index 
    174175      ! 
    175176      INTEGER  :: ji, jj, jh, jc, nhc 
     
    194195                  DO ji = 1,jpi 
    195196                     ! Elevation 
    196                      ana_temp(ji,jj,nhc,1) = ana_temp(ji,jj,nhc,1) + ztemp*sshn(ji,jj)*ssmask (ji,jj)         
    197                      ana_temp(ji,jj,nhc,2) = ana_temp(ji,jj,nhc,2) + ztemp*un_b(ji,jj)*ssumask(ji,jj) 
    198                      ana_temp(ji,jj,nhc,3) = ana_temp(ji,jj,nhc,3) + ztemp*vn_b(ji,jj)*ssvmask(ji,jj) 
     197                     ana_temp(ji,jj,nhc,1) = ana_temp(ji,jj,nhc,1) + ztemp*ssh(ji,jj,Kmm)*ssmask (ji,jj)         
     198                     ana_temp(ji,jj,nhc,2) = ana_temp(ji,jj,nhc,2) + ztemp*uu_b(ji,jj,Kmm)*ssumask(ji,jj) 
     199                     ana_temp(ji,jj,nhc,3) = ana_temp(ji,jj,nhc,3) + ztemp*vv_b(ji,jj,Kmm)*ssvmask(ji,jj) 
    199200                  END DO 
    200201               END DO 
     
    521522   LOGICAL, PUBLIC, PARAMETER ::   lk_diaharm = .FALSE. 
    522523CONTAINS 
    523    SUBROUTINE dia_harm ( kt )     ! Empty routine 
     524   SUBROUTINE dia_harm ( kt, Kmm )     ! Empty routine 
    524525      INTEGER, INTENT( IN ) :: kt   
     526      INTEGER, INTENT( IN ) :: Kmm   
    525527      WRITE(*,*) 'dia_harm: you should not have seen this print' 
    526528   END SUBROUTINE dia_harm 
Note: See TracChangeset for help on using the changeset viewer.