Changeset 11657


Ignore:
Timestamp:
2019-10-07T11:50:50+02:00 (12 months ago)
Author:
dancopsey
Message:

Couple the correct sea ice meltpond variables

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0_fix_cpl_oce_only/src/OCE/SBC/sbccpl.F90

    r11479 r11657  
    156156   INTEGER, PARAMETER ::   jps_wlev   = 32   ! water level  
    157157   INTEGER, PARAMETER ::   jps_fice1  = 33   ! first-order ice concentration (for semi-implicit coupling of atmos-ice fluxes) 
    158    INTEGER, PARAMETER ::   jps_a_p    = 34   ! meltpond area 
     158   INTEGER, PARAMETER ::   jps_a_p    = 34   ! meltpond area fraction 
    159159   INTEGER, PARAMETER ::   jps_ht_p   = 35   ! meltpond thickness 
    160160   INTEGER, PARAMETER ::   jps_kice   = 36   ! sea ice effective conductivity 
     
    23082308      !                                                      !      Ice melt ponds       !  
    23092309      !                                                      ! ------------------------- ! 
    2310       ! needed by Met Office 
     2310      ! needed by Met Office - 1) fraction of ponded ice; 2) local/actual pond depth 
    23112311      IF( ssnd(jps_a_p)%laction .OR. ssnd(jps_ht_p)%laction ) THEN  
    23122312         SELECT CASE( sn_snd_mpnd%cldes)   
     
    23142314            SELECT CASE( sn_snd_mpnd%clcat )   
    23152315            CASE( 'yes' )   
    2316                ztmp3(:,:,1:jpl) =  a_ip(:,:,1:jpl) 
    2317                ztmp4(:,:,1:jpl) =  v_ip(:,:,1:jpl)   
     2316               ztmp3(:,:,1:jpl) =  a_ip_frac(:,:,1:jpl) 
     2317               ztmp4(:,:,1:jpl) =  h_ip(:,:,1:jpl) 
    23182318            CASE( 'no' )   
    23192319               ztmp3(:,:,:) = 0.0   
    23202320               ztmp4(:,:,:) = 0.0   
    23212321               DO jl=1,jpl   
    2322                  ztmp3(:,:,1) = ztmp3(:,:,1) + a_ip(:,:,jpl)   
    2323                  ztmp4(:,:,1) = ztmp4(:,:,1) + v_ip(:,:,jpl)  
     2322                 ztmp3(:,:,1) = ztmp3(:,:,1) + a_ip_frac(:,:,jpl) 
     2323                 ztmp4(:,:,1) = ztmp4(:,:,1) + h_ip(:,:,jpl) 
    23242324               ENDDO   
    23252325            CASE default   ;   CALL ctl_stop( 'sbc_cpl_snd: wrong definition of sn_snd_mpnd%clcat' )   
Note: See TracChangeset for help on using the changeset viewer.