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 8316 for branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/OPA_SRC/TRA – NEMO

Ignore:
Timestamp:
2017-07-11T14:05:05+02:00 (7 years ago)
Author:
clem
Message:

STEP2 (3): remove obsolete features

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/OPA_SRC/TRA/traqsr.F90

    r8313 r8316  
    4848   LOGICAL , PUBLIC ::   ln_qsr_2bd   !: 2 band         light absorption flag 
    4949   LOGICAL , PUBLIC ::   ln_qsr_bio   !: bio-model      light absorption flag 
    50    LOGICAL , PUBLIC ::   ln_qsr_ice   !: light penetration for ice-model 
    5150   INTEGER , PUBLIC ::   nn_chldta    !: use Chlorophyll data (=1) or not (=0) 
    5251   REAL(wp), PUBLIC ::   rn_abs       !: fraction absorbed in the very near surface (RGB & 2 bands) 
     
    269268      END DO 
    270269      ! 
    271       IF( ln_qsr_ice ) THEN      ! sea-ice: store the 1st ocean level attenuation coefficient 
    272          DO jj = 2, jpjm1  
    273             DO ji = fs_2, fs_jpim1   ! vector opt. 
    274                IF( qsr(ji,jj) /= 0._wp ) THEN   ;   fraqsr_1lev(ji,jj) = qsr_hc(ji,jj,1) / ( r1_rau0_rcp * qsr(ji,jj) ) 
    275                ELSE                             ;   fraqsr_1lev(ji,jj) = 1._wp 
    276                ENDIF 
    277             END DO 
    278          END DO 
    279          ! Update haloes since lim_thd needs fraqsr_1lev to be defined everywhere 
    280          CALL lbc_lnk( fraqsr_1lev(:,:), 'T', 1._wp ) 
    281       ENDIF 
     270      ! sea-ice: store the 1st ocean level attenuation coefficient 
     271      DO jj = 2, jpjm1  
     272         DO ji = fs_2, fs_jpim1   ! vector opt. 
     273            IF( qsr(ji,jj) /= 0._wp ) THEN   ;   fraqsr_1lev(ji,jj) = qsr_hc(ji,jj,1) / ( r1_rau0_rcp * qsr(ji,jj) ) 
     274            ELSE                             ;   fraqsr_1lev(ji,jj) = 1._wp 
     275            ENDIF 
     276         END DO 
     277      END DO 
     278      CALL lbc_lnk( fraqsr_1lev(:,:), 'T', 1._wp ) 
    282279      ! 
    283280      IF( iom_use('qsr3d') ) THEN      ! output the shortwave Radiation distribution 
     
    336333      TYPE(FLD_N)        ::   sn_chl   ! informations about the chlorofyl field to be read 
    337334      !! 
    338       NAMELIST/namtra_qsr/  sn_chl, cn_dir, ln_qsr_rgb, ln_qsr_2bd, ln_qsr_bio, ln_qsr_ice, & 
     335      NAMELIST/namtra_qsr/  sn_chl, cn_dir, ln_qsr_rgb, ln_qsr_2bd, ln_qsr_bio, & 
    339336         &                  nn_chldta, rn_abs, rn_si0, rn_si1 
    340337      !!---------------------------------------------------------------------- 
     
    359356         WRITE(numout,*) '      2 band               light penetration       ln_qsr_2bd = ', ln_qsr_2bd 
    360357         WRITE(numout,*) '      bio-model            light penetration       ln_qsr_bio = ', ln_qsr_bio 
    361          WRITE(numout,*) '      light penetration for ice-model (LIM3)       ln_qsr_ice = ', ln_qsr_ice 
    362358         WRITE(numout,*) '      RGB : Chl data (=1) or cst value (=0)        nn_chldta  = ', nn_chldta 
    363359         WRITE(numout,*) '      RGB & 2 bands: fraction of light (rn_si1)    rn_abs     = ', rn_abs 
Note: See TracChangeset for help on using the changeset viewer.