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 5866 for branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssm.F90 – NEMO

Ignore:
Timestamp:
2015-11-06T11:42:35+01:00 (8 years ago)
Author:
gm
Message:

#1613: vvl by default: add ln_linssh and remove key_vvl

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssm.F90

    r5845 r5866  
    7878         ENDIF 
    7979         ! 
    80          IF( lk_vvl )   e3t_m(:,:) = e3t_n(:,:,1) 
     80         IF( .NOT.ln_linssh )   e3t_m(:,:) = e3t_n(:,:,1) 
    8181         ! 
    8282         frq_m(:,:) = fraqsr_1lev(:,:) 
     
    100100            ENDIF 
    101101            ! 
    102             IF( lk_vvl )   e3t_m(:,:) = zcoef * e3t_n(:,:,1) 
     102            IF( .NOT.ln_linssh )   e3t_m(:,:) = zcoef * e3t_n(:,:,1) 
    103103            ! 
    104104            frq_m(:,:) = zcoef * fraqsr_1lev(:,:) 
     
    106106         ELSEIF( MOD( kt - 2 , nn_fsbc ) == 0 ) THEN      !   Initialisation: New mean computation   ! 
    107107            !                                             ! ---------------------------------------- ! 
    108             ssu_m(:,:) = 0.e0      ! reset to zero ocean mean sbc fields 
    109             ssv_m(:,:) = 0.e0 
    110             sst_m(:,:) = 0.e0 
    111             sss_m(:,:) = 0.e0 
    112             ssh_m(:,:) = 0.e0 
    113             IF( lk_vvl )   e3t_m(:,:) = 0.e0 
    114             frq_m(:,:) = 0.e0 
     108            ssu_m(:,:) = 0._wp     ! reset to zero ocean mean sbc fields 
     109            ssv_m(:,:) = 0._wp 
     110            sst_m(:,:) = 0._wp 
     111            sss_m(:,:) = 0._wp 
     112            ssh_m(:,:) = 0._wp 
     113            IF( .NOT.ln_linssh )   e3t_m(:,:) = 0._wp 
     114            frq_m(:,:) = 0._wp 
    115115         ENDIF 
    116116         !                                                ! ---------------------------------------- ! 
     
    128128         ENDIF 
    129129         ! 
    130          IF( lk_vvl )   e3t_m(:,:) = e3t_m(:,:) + e3t_n(:,:,1) 
     130         IF( .NOT.ln_linssh )   e3t_m(:,:) = e3t_m(:,:) + e3t_n(:,:,1) 
    131131         ! 
    132132         frq_m(:,:) =   frq_m(:,:) + fraqsr_1lev(:,:) 
     
    141141            ssv_m(:,:) = ssv_m(:,:) * zcoef           ! 
    142142            ssh_m(:,:) = ssh_m(:,:) * zcoef           ! mean SSH             [m] 
    143             IF( lk_vvl )   e3t_m(:,:) = e3t_m(:,:) * zcoef   ! mean vertical scale factor [m] 
     143            IF( .NOT.ln_linssh )   e3t_m(:,:) = e3t_m(:,:) * zcoef   ! mean vertical scale factor [m] 
    144144            frq_m(:,:) = frq_m(:,:) * zcoef   ! mean fraction of solar net radiation absorbed in the 1st T level [-] 
    145145            ! 
     
    159159            CALL iom_rstput( kt, nitrst, numrow, 'sss_m'  , sss_m  ) 
    160160            CALL iom_rstput( kt, nitrst, numrow, 'ssh_m'  , ssh_m  ) 
    161             IF( lk_vvl )   CALL iom_rstput( kt, nitrst, numrow, 'e3t_m'  , e3t_m  ) 
     161            IF( .NOT.ln_linssh )   CALL iom_rstput( kt, nitrst, numrow, 'e3t_m'  , e3t_m  ) 
    162162            CALL iom_rstput( kt, nitrst, numrow, 'frq_m'  , frq_m  ) 
    163163            ! 
     
    172172         CALL iom_put( 'sss_m', sss_m ) 
    173173         CALL iom_put( 'ssh_m', ssh_m ) 
    174          IF( lk_vvl )   CALL iom_put( 'e3t_m', e3t_m ) 
     174         IF( .NOT.ln_linssh )   CALL iom_put( 'e3t_m', e3t_m ) 
    175175         CALL iom_put( 'frq_m', frq_m ) 
    176176      ENDIF 
     
    209209            CALL iom_get( numror, jpdom_autoglo, 'sss_m'  , sss_m  )   !   "         "    salinity    (T-point) 
    210210            CALL iom_get( numror, jpdom_autoglo, 'ssh_m'  , ssh_m  )   !   "         "    height      (T-point) 
    211             IF( lk_vvl ) CALL iom_get( numror, jpdom_autoglo, 'e3t_m', e3t_m ) 
     211            IF( .NOT.ln_linssh )  CALL iom_get( numror, jpdom_autoglo, 'e3t_m', e3t_m ) 
    212212            ! fraction of solar net radiation absorbed in 1st T level 
    213213            IF( iom_varid( numror, 'frq_m', ldstop = .FALSE. ) > 0 ) THEN 
     
    226226               sss_m(:,:) = zcoef * sss_m(:,:) 
    227227               ssh_m(:,:) = zcoef * ssh_m(:,:) 
    228                IF( lk_vvl )   e3t_m(:,:) = zcoef * e3t_m(:,:) 
     228               IF( .NOT.ln_linssh )   e3t_m(:,:) = zcoef * e3t_m(:,:) 
    229229               frq_m(:,:) = zcoef * frq_m(:,:) 
    230230            ELSE 
     
    244244         sss_m(:,:) = tsn(:,:,1,jp_sal) 
    245245         ssh_m(:,:) = sshn(:,:) 
    246          IF( lk_vvl )   e3t_m(:,:) = e3t_n(:,:,1) 
     246         IF( .NOT.ln_linssh )   e3t_m(:,:) = e3t_n(:,:,1) 
    247247         frq_m(:,:) = 1._wp 
    248248         ! 
Note: See TracChangeset for help on using the changeset viewer.