Changeset 8063


Ignore:
Timestamp:
2017-05-24T15:05:23+02:00 (3 years ago)
Author:
vancop
Message:

a few more melt pond corrections

Location:
branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90

    r8061 r8063  
    665665            !------------------------------------------             
    666666            !  Place part of the melt pond volume into the ocean.  
    667             IF ( ( nn_pnd_scheme > 0 ) .AND. ( ( nn_pnd_cpl .EQ. 1 ) .OR. ( nn_pnd_cpl .EQ. 3 ) ) ) 
     667            IF ( ( nn_pnd_scheme > 0 ) .AND. ( ( nn_pnd_cpl == 1 ) .OR. ( nn_pnd_cpl == 3 ) ) ) 
    668668               wfx_pnd(ji,jj) = wfx_pnd(ji,jj) + ( rhofw * vprdg(ij) * ( 1._wp - rn_fpondrdg )   &  
    669669               &                                 + rhofw * vprft(ij) * ( 1._wp - rn_fpondrft ) ) * r1_rdtice  ! fresh water source for ocean 
  • branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limmp.F90

    r8061 r8063  
    193193       REAL(wp)                            :: z1_rhofw          ! inverse freshwater density 
    194194       REAL(wp)                            :: z1_zpnd_aspect    ! inverse pond aspect ratio 
     195       REAL(wp)                            :: zvold             ! dummy pond volume 
    195196 
    196197       INTEGER                             :: ji, jj, jl, ij    ! loop indices 
     
    254255               h_ip(ji,jj,jl)      = 0._wp 
    255256 
    256                IF ( ( nn_pnd_cpl .EQ. 1 ) .OR. ( nn_pnd_cpl .EQ. 3 ) ) & !--- Give freshwater to the ocean 
     257               IF ( ( nn_pnd_cpl == 1 ) .OR. ( nn_pnd_cpl == 3 ) ) & !--- Give freshwater to the ocean 
    257258                  wfx_pnd(ji,jj)   = wfx_pnd(ji,jj) + v_ip(ji,jj,jl)  
    258259 
     
    272273               !--- Dump meltwater due to refreezing ( of course this is wrong 
    273274               !--- but this parameterization is too simple ) 
    274                IF ( ( nn_pnd_cpl .EQ. 1 ) .OR. ( nn_pnd_cpl .EQ. 3 ) ) THEN 
     275               IF ( ( nn_pnd_cpl == 1 ) .OR. ( nn_pnd_cpl == 3 ) ) THEN 
    275276 
    276277                  wfx_pnd(ji,jj)   = wfx_pnd(ji,jj) + rhofw * ( v_ip(ji,jj,jl) - zvpold ) * r1_rdtice 
  • branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limvar.F90

    r8061 r8063  
    9292 
    9393      ! MV MP 2016 
    94       IF ( nn_pnd_scheme > 0 ) 
     94      IF ( nn_pnd_scheme > 0 ) THEN 
    9595         at_ip(:,:) = SUM( a_ip, dim=3 ) 
    9696         vt_ip(:,:) = SUM( v_ip, dim=3 ) 
  • branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limwri.F90

    r8061 r8063  
    171171      CALL iom_put( "vfxice"     , wfx_ice * ztmp       )        ! total ice growth/melt  
    172172 
     173      IF ( nn_pnd_scheme >= 0 ) & 
     174         CALL iom_put( "vfxpnd"  , wfx_pnd * ztmp       )        ! melt pond water flux 
     175 
    173176      IF ( iom_use( "vfxthin" ) ) THEN   ! ice production for open water + thin ice (<20cm) => comparable to observations   
    174177         WHERE( htm_i(:,:) < 0.2 .AND. htm_i(:,:) > 0. ) ; z2d = wfx_bog 
     
    367370      CALL histwrite( kid, "vfxbom", kt, wfx_bom        , jpi*jpj, (/1/) ) 
    368371      CALL histwrite( kid, "vfxsum", kt, wfx_sum        , jpi*jpj, (/1/) ) 
     372      IF ( nn_pnd_scheme >= 0 ) & 
     373         CALL histwrite( kid, "vfxpnd", kt, wfx_pnd     , jpi*jpj, (/1/) ) 
    369374 
    370375      CALL histwrite( kid, "sithicat", kt, ht_i        , jpi*jpj*jpl, (/1/) )     
Note: See TracChangeset for help on using the changeset viewer.