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 11395 for NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/DYN – NEMO

Ignore:
Timestamp:
2019-08-02T16:19:00+02:00 (5 years ago)
Author:
mathiot
Message:

ENHANCE-02_ISF_nemo : Initial commit isf simplification (add ISF directory, moved isf routine in and split isf cavity and isf parametrisation, ...) (ticket #2142)

Location:
NEMO/branches/2019/ENHANCE-02_ISF_nemo
Files:
5 edited
1 copied

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/DYN/divhor.F90

    r10425 r11395  
    2222   USE sbc_oce, ONLY : ln_rnf, ln_isf ! surface boundary condition: ocean 
    2323   USE sbcrnf          ! river runoff  
    24    USE sbcisf          ! ice shelf 
     24   USE isfhdiv         ! ice shelf 
    2525   USE iscplhsb        ! ice sheet / ocean coupling 
    2626   USE iscplini        ! ice sheet / ocean coupling 
    … …  
    100100      !  
    101101#endif 
    102       IF( ln_isf )   CALL sbc_isf_div( hdivn )      !==  ice shelf  ==!   (update hdivn field) 
    103       ! 
    104       IF( ln_iscpl .AND. ln_hsb )   CALL iscpl_div( hdivn ) !==  ice sheet  ==!   (update hdivn field) 
     102      IF( ln_isf )   CALL isf_hdiv( hdivn )      !==  ice shelf  ==!   (update hdivn field) 
    105103      ! 
    106104      CALL lbc_lnk( 'divhor', hdivn, 'T', 1. )   !   (no sign change) 

  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/DYN/dynhpg.F90

    r10491 r11395  
    3131   !!---------------------------------------------------------------------- 
    3232   USE oce             ! ocean dynamics and tracers 
     33   USE isf             ! ice shelf 
    3334   USE sbc_oce         ! surface variable (only for the flag with ice shelf) 
    3435   USE dom_oce         ! ocean space and time domain 

  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/DYN/dynnxt.F90

    r10425 r11395  
    2929   USE sbc_oce        ! Surface boundary condition: ocean fields 
    3030   USE sbcrnf         ! river runoffs 
    31    USE sbcisf         ! ice shelf 
     31   USE isfnxt 
    3232   USE phycst         ! physical constants 
    3333   USE dynadv         ! dynamics: vector invariant versus flux form 
    … …  
    241241               ENDIF 
    242242            END IF 
    243  
    244             IF ( ln_isf ) THEN   ! if ice shelf melting 
    245                DO jk = 1, jpkm1 ! Deal with isf separetely, as can be through depth too 
    246                   DO jj = 1, jpj 
    247                      DO ji = 1, jpi 
    248                         IF( misfkt(ji,jj) <=jk .and. jk < misfkb(ji,jj)  ) THEN 
    249                            e3t_b(ji,jj,jk) = e3t_b(ji,jj,jk) - zcoef * ( fwfisf_b(ji,jj) - fwfisf(ji,jj) ) & 
    250                                 &          * ( e3t_n(ji,jj,jk) * r1_hisf_tbl(ji,jj) ) * tmask(ji,jj,jk) 
    251                         ELSEIF ( jk==misfkb(ji,jj) ) THEN 
    252                            e3t_b(ji,jj,jk) = e3t_b(ji,jj,jk) - zcoef * ( fwfisf_b(ji,jj) - fwfisf(ji,jj) ) & 
    253                                 &          * ( e3t_n(ji,jj,jk) * r1_hisf_tbl(ji,jj) ) * ralpha(ji,jj) * tmask(ji,jj,jk) 
    254                         ENDIF 
    255                      END DO 
    256                   END DO 
    257                END DO 
    258             END IF 
     243            ! 
     244            ! ice shelf melting 
     245            IF ( ln_isf ) CALL isf_dynnxt( zcoef ) 
    259246            ! 
    260247            IF( ln_dynadv_vec ) THEN      ! Asselin filter applied on velocity 

  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/DYN/dynspg_ts.F90

    r10742 r11395  
    3333   USE zdf_oce         ! vertical physics: variables 
    3434   USE zdfdrg          ! vertical physics: top/bottom drag coef. 
    35    USE sbcisf          ! ice shelf variable (fwfisf) 
     35   USE isf             ! ice shelf variable (fwfisf) 
    3636   USE sbcapr          ! surface boundary condition: atmospheric pressure 
    3737   USE dynadv    , ONLY: ln_dynadv_vec 
    … …  
    632632      !                                         ! Surface net water flux and rivers 
    633633      IF (ln_bt_fw) THEN 
    634          zssh_frc(:,:) = r1_rau0 * ( emp(:,:) - rnf(:,:) + fwfisf(:,:) ) 
     634         zssh_frc(:,:) = r1_rau0 * ( emp(:,:) - rnf(:,:) + fwfisf_cav(:,:) + fwfisf_par(:,:) ) 
    635635      ELSE 
    636636         zztmp = r1_rau0 * r1_2 
    637637         zssh_frc(:,:) = zztmp * (  emp(:,:) + emp_b(:,:) - rnf(:,:) - rnf_b(:,:)   & 
    638                 &                 + fwfisf(:,:) + fwfisf_b(:,:)                     ) 
     638                &                 + fwfisf_cav(:,:) + fwfisf_cav_b(:,:)             & 
     639                &                 + fwfisf_par(:,:) + fwfisf_par_b(:,:)             ) 
    639640      ENDIF 
    640641      ! 

  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/DYN/sshwzv.F90

    r11293 r11395  
    1717   !!---------------------------------------------------------------------- 
    1818   USE oce            ! ocean dynamics and tracers variables 
     19   USE isf            ! ice shelf 
    1920   USE dom_oce        ! ocean space and time domain variables  
    2021   USE sbc_oce        ! surface boundary condition: ocean 
    … …  
    256257            sshb(:,:) = sshb(:,:) - zcoef * (     emp_b(:,:) - emp   (:,:)   & 
    257258               &                             -    rnf_b(:,:) + rnf   (:,:)   & 
    258                &                             + fwfisf_b(:,:) - fwfisf(:,:)   ) * ssmask(:,:) 
     259               &                             + fwfisf_cav_b(:,:) - fwfisf_cav(:,:)   & 
     260               &                             + fwfisf_par_b(:,:) - fwfisf_par(:,:)   ) * ssmask(:,:) 
    259261         ENDIF 
    260262         sshn(:,:) = ssha(:,:)                              ! now <-- after 
Note: See TracChangeset for help on using the changeset viewer.