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

a few more melt pond corrections

File:
1 edited

Legend:

Unmodified
Added
Removed
  • 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 
Note: See TracChangeset for help on using the changeset viewer.