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 12744 for NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icedyn_rdgrft.F90 – NEMO

Ignore:
Timestamp:
2020-04-12T16:11:29+02:00 (4 years ago)
Author:
clem
Message:

make sure all pond lids are set to 0 when not using this option

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icedyn_rdgrft.F90

    r12720 r12744  
    569569                  aprdg2(ji) = a_ip_2d(ji,jl1) * afrdg * hi_hrdg(ji,jl1) 
    570570                  vprdg (ji) = v_ip_2d(ji,jl1) * afrdg 
    571                   vlrdg (ji) = v_il_2d(ji,jl1) * afrdg 
    572571                  aprft1     = a_ip_2d(ji,jl1) * afrft 
    573572                  aprft2(ji) = a_ip_2d(ji,jl1) * afrft * hi_hrft 
    574573                  vprft (ji) = v_ip_2d(ji,jl1) * afrft 
    575                   vlrft (ji) = v_il_2d(ji,jl1) * afrft 
     574                  IF ( ln_pnd_lids ) THEN 
     575                     vlrdg (ji) = v_il_2d(ji,jl1) * afrdg 
     576                     vlrft (ji) = v_il_2d(ji,jl1) * afrft 
     577                  ENDIF 
    576578               ENDIF 
    577579 
     
    603605                  a_ip_2d(ji,jl1) = a_ip_2d(ji,jl1) - aprdg1    - aprft1 
    604606                  v_ip_2d(ji,jl1) = v_ip_2d(ji,jl1) - vprdg(ji) - vprft(ji) 
    605                   v_il_2d(ji,jl1) = v_il_2d(ji,jl1) - vlrdg(ji) - vlrft(ji) 
     607                  IF ( ln_pnd_lids ) THEN 
     608                     v_il_2d(ji,jl1) = v_il_2d(ji,jl1) - vlrdg(ji) - vlrft(ji) 
     609                  ENDIF 
    606610               ENDIF 
    607611            ENDIF 
     
    700704                     a_ip_2d (ji,jl2) = a_ip_2d(ji,jl2) + (   aprdg2(ji) * rn_fpndrdg * farea         &  
    701705                        &                                   + aprft2(ji) * rn_fpndrft * zswitch(ji)   ) 
    702                      v_il_2d (ji,jl2) = v_il_2d(ji,jl2) + (   vlrdg (ji) * rn_fpndrdg * fvol   (ji)   & 
    703                         &                                   + vlrft (ji) * rn_fpndrft * zswitch(ji)   ) 
     706                     IF ( ln_pnd_lids ) THEN 
     707                        v_il_2d (ji,jl2) = v_il_2d(ji,jl2) + (   vlrdg (ji) * rn_fpndrdg * fvol   (ji)   & 
     708                           &                                   + vlrft (ji) * rn_fpndrft * zswitch(ji)   ) 
     709                     ENDIF 
    704710                  ENDIF 
    705711                   
Note: See TracChangeset for help on using the changeset viewer.