Changeset 11466


Ignore:
Timestamp:
2019-08-21T09:18:21+02:00 (13 months ago)
Author:
davestorkey
Message:

dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap: Bug fixes for ssh_atf routine in sshwzv.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/OCE/DYN/sshwzv.F90

    r11426 r11466  
    226226      !! ** Purpose :   achieve the sea surface  height time stepping by  
    227227      !!              applying Asselin time filter and swapping the arrays 
    228       !!              ssh(:,:,Kaa)  already computed in ssh_nxt   
     228      !!              pssh(:,:,Kaa)  already computed in ssh_nxt   
    229229      !! 
    230230      !! ** Method  : - apply Asselin time fiter to now ssh (excluding the forcing 
    231231      !!              from the filter, see Leclair and Madec 2010) and swap : 
    232       !!                ssh(:,:,Kmm) = ssh(:,:,Kaa) + atfp * ( ssh(:,:,Kbb) -2 ssh(:,:,Kmm) + ssh(:,:,Kaa) ) 
     232      !!                pssh(:,:,Kmm) = pssh(:,:,Kaa) + atfp * ( pssh(:,:,Kbb) -2 pssh(:,:,Kmm) + pssh(:,:,Kaa) ) 
    233233      !!                            - atfp * rdt * ( emp_b - emp ) / rau0 
    234       !!                ssh(:,:,Kmm) = ssh(:,:,Kaa) 
    235       !! 
    236       !! ** action  : - ssh(:,:,Kbb), ssh(:,:,Kmm)   : before & now sea surface height 
     234      !!                pssh(:,:,Kmm) = pssh(:,:,Kaa) 
     235      !! 
     236      !! ** action  : - pssh(:,:,Kbb), pssh(:,:,Kmm)   : before & now sea surface height 
    237237      !!                               ready for the next time step 
    238238      !! 
     
    256256      IF ( .NOT.( neuler == 0 .AND. kt == nit000 ) ) THEN   ! Only do time filtering for leapfrog timesteps 
    257257         !                                                  ! filtered "now" field 
    258          ssh(:,:,Kmm) = ssh(:,:,Kmm) + atfp * ( ssh(:,:,Kbb) - 2 * ssh(:,:,Kmm) + ssh(:,:,Kaa) ) 
     258         pssh(:,:,Kmm) = pssh(:,:,Kmm) + atfp * ( pssh(:,:,Kbb) - 2 * pssh(:,:,Kmm) + pssh(:,:,Kaa) ) 
    259259         IF( .NOT.ln_linssh ) THEN                          ! "now" <-- with forcing removed 
    260260            zcoef = atfp * rdt * r1_rau0 
    261             ssh(:,:,Kmm) = ssh(:,:,Kmm) - zcoef * (     emp_b(:,:) - emp   (:,:)   & 
     261            pssh(:,:,Kmm) = pssh(:,:,Kmm) - zcoef * (     emp_b(:,:) - emp   (:,:)   & 
    262262               &                             -    rnf_b(:,:) + rnf   (:,:)   & 
    263263               &                             + fwfisf_b(:,:) - fwfisf(:,:)   ) * ssmask(:,:) 
     
    265265      ENDIF 
    266266      ! 
    267       IF(ln_ctl)   CALL prt_ctl( tab2d_1=ssh(:,:,Kmm), clinfo1=' ssh(:,:,Kmm)  - : ', mask1=tmask ) 
     267      IF(ln_ctl)   CALL prt_ctl( tab2d_1=pssh(:,:,Kmm), clinfo1=' pssh(:,:,Kmm)  - : ', mask1=tmask ) 
    268268      ! 
    269269      IF( ln_timing )   CALL timing_stop('ssh_atf') 
Note: See TracChangeset for help on using the changeset viewer.