Ignore:
Timestamp:
2018-12-19T22:54:16+01:00 (2 years ago)
Author:
smasson
Message:

trunk: merge back dev_r10164_HPC09_ESIWACE_PREP_MERGE@10424 into the trunk

Location:
NEMO/trunk/src/OCE/SBC
Files:
17 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/SBC/cpl_oasis3.F90

    r10404 r10425  
    456456         !--- Fill the overlap areas and extra hallows (mpp) 
    457457         !--- check periodicity conditions (all cases) 
    458          IF( .not. llfisrt )   CALL lbc_lnk( pdata(:,:,jc), srcv(kid)%clgrid, srcv(kid)%nsgn )    
     458         IF( .not. llfisrt ) THEN 
     459            CALL lbc_lnk( 'cpl_oasis3', pdata(:,:,jc), srcv(kid)%clgrid, srcv(kid)%nsgn )    
     460         ENDIF 
    459461  
    460462      ENDDO 
  • NEMO/trunk/src/OCE/SBC/fldread.F90

    r10068 r10425  
    669669               IF( sdjf%ln_tint ) THEN 
    670670                  CALL iom_get( sdjf%num, sdjf%clvar, sdjf%fdta(2,2,1,2), sdjf%nrec_a(1) ) 
    671                   CALL lbc_lnk( sdjf%fdta(:,:,1,2),'Z',1. ) 
     671                  CALL lbc_lnk( 'fldread', sdjf%fdta(:,:,1,2),'Z',1. ) 
    672672               ELSE 
    673673                  CALL iom_get( sdjf%num, sdjf%clvar, sdjf%fnow(2,2,1  ), sdjf%nrec_a(1) ) 
    674                   CALL lbc_lnk( sdjf%fnow(:,:,1  ),'Z',1. ) 
     674                  CALL lbc_lnk( 'fldread', sdjf%fnow(:,:,1  ),'Z',1. ) 
    675675               ENDIF 
    676676            ELSE 
     
    683683               IF( sdjf%ln_tint ) THEN 
    684684                  CALL iom_get( sdjf%num, jpdom_unknown, sdjf%clvar, sdjf%fdta(2,2,:,2), sdjf%nrec_a(1) ) 
    685                   CALL lbc_lnk( sdjf%fdta(:,:,:,2),'Z',1. ) 
     685                  CALL lbc_lnk( 'fldread', sdjf%fdta(:,:,:,2),'Z',1. ) 
    686686               ELSE 
    687687                  CALL iom_get( sdjf%num, jpdom_unknown, sdjf%clvar, sdjf%fnow(2,2,:  ), sdjf%nrec_a(1) ) 
    688                   CALL lbc_lnk( sdjf%fnow(:,:,:  ),'Z',1. ) 
     688                  CALL lbc_lnk( 'fldread', sdjf%fnow(:,:,:  ),'Z',1. ) 
    689689               ENDIF 
    690690            ELSE 
  • NEMO/trunk/src/OCE/SBC/geo2ocean.F90

    r10370 r10425  
    152152         &      gsinv(jpi,jpj), gcosv(jpi,jpj),   &   
    153153         &      gsinf(jpi,jpj), gcosf(jpi,jpj), STAT=ierr ) 
    154       IF(lk_mpp)   CALL mpp_sum( ierr ) 
     154      CALL mpp_sum( 'geo2ocean', ierr ) 
    155155      IF( ierr /= 0 )   CALL ctl_stop( 'angle: unable to allocate arrays' ) 
    156156      ! 
     
    276276      ! =========================== ! 
    277277      !           ! lateral boundary cond.: T-, U-, V-, F-pts, sgn 
    278       CALL lbc_lnk_multi( gcost, 'T', -1., gsint, 'T', -1., gcosu, 'U', -1., gsinu, 'U', -1., &  
     278      CALL lbc_lnk_multi( 'geo2ocean', gcost, 'T', -1., gsint, 'T', -1., gcosu, 'U', -1., gsinu, 'U', -1., &  
    279279                      &   gcosv, 'V', -1., gsinv, 'V', -1., gcosf, 'F', -1., gsinf, 'F', -1.  ) 
    280280      ! 
     
    304304         ALLOCATE( gsinlon(jpi,jpj,4) , gcoslon(jpi,jpj,4) ,   & 
    305305            &      gsinlat(jpi,jpj,4) , gcoslat(jpi,jpj,4) , STAT=ierr ) 
    306          IF( lk_mpp    )   CALL mpp_sum( ierr ) 
     306         CALL mpp_sum( 'geo2ocean', ierr ) 
    307307         IF( ierr /= 0 )   CALL ctl_stop('geo2oce: unable to allocate arrays' ) 
    308308      ENDIF 
     
    381381         ALLOCATE( gsinlon(jpi,jpj,4) , gcoslon(jpi,jpj,4) ,   & 
    382382            &      gsinlat(jpi,jpj,4) , gcoslat(jpi,jpj,4) , STAT=ierr ) 
    383          IF( lk_mpp    )   CALL mpp_sum( ierr ) 
     383         CALL mpp_sum( 'geo2ocean', ierr ) 
    384384         IF( ierr /= 0 )   CALL ctl_stop('oce2geo: unable to allocate arrays' ) 
    385385      ENDIF 
  • NEMO/trunk/src/OCE/SBC/sbc_ice.F90

    r10068 r10425  
    150150 
    151151      sbc_ice_alloc = MAXVAL( ierr ) 
    152       IF( lk_mpp            )   CALL mpp_sum ( sbc_ice_alloc ) 
     152      CALL mpp_sum ( 'sbc_ice', sbc_ice_alloc ) 
    153153      IF( sbc_ice_alloc > 0 )   CALL ctl_warn('sbc_ice_alloc: allocation of arrays failed') 
    154154   END FUNCTION sbc_ice_alloc 
     
    197197      ALLOCATE( snwice_mass(jpi,jpj) , snwice_mass_b(jpi,jpj), snwice_fmass(jpi,jpj) , STAT=ierr(1) ) 
    198198      sbc_ice_alloc = MAXVAL( ierr ) 
    199       IF( lk_mpp            )   CALL mpp_sum ( sbc_ice_alloc ) 
     199      CALL mpp_sum ( 'sbc_ice', sbc_ice_alloc ) 
    200200      IF( sbc_ice_alloc > 0 )   CALL ctl_warn('sbc_ice_alloc: allocation of arrays failed') 
    201201   END FUNCTION sbc_ice_alloc 
  • NEMO/trunk/src/OCE/SBC/sbc_oce.F90

    r10068 r10425  
    186186         ! 
    187187      sbc_oce_alloc = MAXVAL( ierr ) 
    188       IF( lk_mpp            )   CALL mpp_sum ( sbc_oce_alloc ) 
     188      CALL mpp_sum ( 'sbc_oce', sbc_oce_alloc ) 
    189189      IF( sbc_oce_alloc > 0 )   CALL ctl_warn('sbc_oce_alloc: allocation of arrays failed') 
    190190      ! 
     
    216216         END DO 
    217217      END DO 
    218       CALL lbc_lnk( wndm(:,:) , 'T', 1. ) 
     218      CALL lbc_lnk( 'sbc_oce', wndm(:,:) , 'T', 1. ) 
    219219      ! 
    220220   END SUBROUTINE sbc_tau2wnd 
  • NEMO/trunk/src/OCE/SBC/sbcapr.F90

    r10068 r10425  
    9494      ! 
    9595      IF( ln_ref_apr ) THEN                        !* Compute whole inner domain mean masked ocean surface 
    96          tarea = glob_sum( e1e2t(:,:) ) 
     96         tarea = glob_sum( 'sbcapr', e1e2t(:,:) ) 
    9797         IF(lwp) WRITE(numout,*) '         Variable ref. Patm computed over a ocean surface of ', tarea*1e-6, 'km2' 
    9898      ELSE 
     
    141141         ! 
    142142         !                                                  !* update the reference atmospheric pressure (if necessary) 
    143          IF( ln_ref_apr )   rn_pref = glob_sum( sf_apr(1)%fnow(:,:,1) * e1e2t(:,:) ) / tarea 
     143         IF( ln_ref_apr )   rn_pref = glob_sum( 'sbcapr', sf_apr(1)%fnow(:,:,1) * e1e2t(:,:) ) / tarea 
    144144         ! 
    145145         !                                                  !* Patm related forcing at kt 
  • NEMO/trunk/src/OCE/SBC/sbcblk.F90

    r10190 r10425  
    147147         &      cdn_oce(jpi,jpj), chn_oce(jpi,jpj), cen_oce(jpi,jpj), STAT=sbc_blk_alloc ) 
    148148      ! 
    149       IF( lk_mpp             )   CALL mpp_sum ( sbc_blk_alloc ) 
    150       IF( sbc_blk_alloc /= 0 )   CALL ctl_warn('sbc_blk_alloc: failed to allocate arrays') 
     149      CALL mpp_sum ( 'sbcblk', sbc_blk_alloc ) 
     150      IF( sbc_blk_alloc /= 0 )   CALL ctl_stop( 'STOP', 'sbc_blk_alloc: failed to allocate arrays' ) 
    151151   END FUNCTION sbc_blk_alloc 
    152152 
     
    236236      !Activated wave module but neither drag nor stokes drift activated 
    237237         IF ( .NOT.(ln_cdgw .OR. ln_sdw .OR. ln_tauwoc .OR. ln_stcor ) )   THEN 
    238             CALL ctl_warn( 'Ask for wave coupling but ln_cdgw=F, ln_sdw=F, ln_tauwoc=F, ln_stcor=F') 
     238            CALL ctl_stop( 'STOP',  'Ask for wave coupling but ln_cdgw=F, ln_sdw=F, ln_tauwoc=F, ln_stcor=F' ) 
    239239      !drag coefficient read from wave model definable only with mfs bulk formulae and core  
    240240         ELSEIF (ln_cdgw .AND. .NOT. ln_NCAR )       THEN        
     
    407407         END DO 
    408408      END DO 
    409       CALL lbc_lnk_multi( zwnd_i, 'T', -1., zwnd_j, 'T', -1. ) 
     409      CALL lbc_lnk_multi( 'sbcblk', zwnd_i, 'T', -1., zwnd_j, 'T', -1. ) 
    410410      ! ... scalar wind ( = | U10m - U_oce | ) at T-point (masked) 
    411411      wndm(:,:) = SQRT(  zwnd_i(:,:) * zwnd_i(:,:)   & 
     
    485485         END DO 
    486486      END DO 
    487       CALL lbc_lnk_multi( utau, 'U', -1., vtau, 'V', -1. ) 
     487      CALL lbc_lnk_multi( 'sbcblk', utau, 'U', -1., vtau, 'V', -1. ) 
    488488 
    489489      !  Turbulent fluxes over ocean 
     
    727727         END DO 
    728728      END DO 
    729       CALL lbc_lnk( wndm_ice, 'T',  1. ) 
     729      CALL lbc_lnk( 'sbcblk', wndm_ice, 'T',  1. ) 
    730730      ! 
    731731      ! Make ice-atm. drag dependent on ice concentration 
     
    761761         END DO 
    762762      END DO 
    763       CALL lbc_lnk_multi( utau_ice, 'U', -1., vtau_ice, 'V', -1. ) 
     763      CALL lbc_lnk_multi( 'sbcblk', utau_ice, 'U', -1., vtau_ice, 'V', -1. ) 
    764764      ! 
    765765      ! 
     
    11861186         END DO 
    11871187      END DO 
    1188       CALL lbc_lnk_multi( Cd, 'T',  1., Ch, 'T', 1. ) 
     1188      CALL lbc_lnk_multi( 'sbcblk', Cd, 'T',  1., Ch, 'T', 1. ) 
    11891189      ! 
    11901190   END SUBROUTINE Cdn10_Lupkes2015 
  • NEMO/trunk/src/OCE/SBC/sbccpl.F90

    r10404 r10425  
    224224 
    225225      sbc_cpl_alloc = MAXVAL( ierr ) 
    226       IF( lk_mpp            )   CALL mpp_sum ( sbc_cpl_alloc ) 
     226      CALL mpp_sum ( 'sbccpl', sbc_cpl_alloc ) 
    227227      IF( sbc_cpl_alloc > 0 )   CALL ctl_warn('sbc_cpl_alloc: allocation of arrays failed') 
    228228      ! 
     
    11621162                  END DO 
    11631163               END DO 
    1164                CALL lbc_lnk_multi( frcv(jpr_otx1)%z3(:,:,1), 'U',  -1., frcv(jpr_oty1)%z3(:,:,1), 'V',  -1. ) 
     1164               CALL lbc_lnk_multi( 'sbccpl', frcv(jpr_otx1)%z3(:,:,1), 'U',  -1., frcv(jpr_oty1)%z3(:,:,1), 'V',  -1. ) 
    11651165            ENDIF 
    11661166            llnewtx = .TRUE. 
     
    11891189               END DO 
    11901190            END DO 
    1191             CALL lbc_lnk( frcv(jpr_taum)%z3(:,:,1), 'T', 1. ) 
     1191            CALL lbc_lnk( 'sbccpl', frcv(jpr_taum)%z3(:,:,1), 'T', 1. ) 
    11921192            llnewtau = .TRUE. 
    11931193         ELSE 
     
    15701570         END SELECT 
    15711571         IF( srcv(jpr_itx1)%clgrid /= 'U' ) THEN  
    1572             CALL lbc_lnk_multi( p_taui, 'U',  -1., p_tauj, 'V',  -1. ) 
     1572            CALL lbc_lnk_multi( 'sbccpl', p_taui, 'U',  -1., p_tauj, 'V',  -1. ) 
    15731573         ENDIF 
    15741574          
     
    23352335                  END DO 
    23362336               END DO 
    2337                CALL lbc_lnk_multi( zitx1, 'T', -1., zity1, 'T', -1. ) 
     2337               CALL lbc_lnk_multi( 'sbccpl', zitx1, 'T', -1., zity1, 'T', -1. ) 
    23382338            CASE( 'mixed oce-ice'        )      ! Ocean and Ice on C-grid ==> T 
    23392339               DO jj = 2, jpjm1 
     
    23462346               END DO 
    23472347            END SELECT 
    2348             CALL lbc_lnk_multi( zotx1, ssnd(jps_ocx1)%clgrid, -1.,  zoty1, ssnd(jps_ocy1)%clgrid, -1. ) 
     2348            CALL lbc_lnk_multi( 'sbccpl', zotx1, ssnd(jps_ocx1)%clgrid, -1.,  zoty1, ssnd(jps_ocy1)%clgrid, -1. ) 
    23492349            ! 
    23502350         ENDIF 
     
    24182418                END DO 
    24192419             END DO 
    2420              CALL lbc_lnk_multi( zitx1, 'T', -1.,  zity1, 'T', -1. )  
     2420             CALL lbc_lnk_multi( 'sbccpl', zitx1, 'T', -1.,  zity1, 'T', -1. )  
    24212421          CASE( 'mixed oce-ice'        )      ! Ocean and Ice on C-grid ==> T   
    24222422             DO jj = 2, jpjm1  
     
    24292429             END DO 
    24302430          END SELECT 
    2431          CALL lbc_lnk_multi( zotx1, ssnd(jps_ocxw)%clgrid, -1., zoty1, ssnd(jps_ocyw)%clgrid, -1. )  
     2431         CALL lbc_lnk_multi( 'sbccpl', zotx1, ssnd(jps_ocxw)%clgrid, -1., zoty1, ssnd(jps_ocyw)%clgrid, -1. )  
    24322432         !  
    24332433         !  
  • NEMO/trunk/src/OCE/SBC/sbcdcy.F90

    r10068 r10425  
    4444            &      rtmd(jpi,jpj) , rdawn(jpi,jpj) , rdusk(jpi,jpj) , rscal(jpi,jpj) , STAT=sbc_dcy_alloc ) 
    4545            ! 
    46          IF( lk_mpp             )   CALL mpp_sum ( sbc_dcy_alloc ) 
    47          IF( sbc_dcy_alloc /= 0 )   CALL ctl_warn('sbc_dcy_alloc: failed to allocate arrays') 
     46         CALL mpp_sum ( 'sbcdcy', sbc_dcy_alloc ) 
     47         IF( sbc_dcy_alloc /= 0 )   CALL ctl_stop( 'STOP', 'sbc_dcy_alloc: failed to allocate arrays' ) 
    4848      END FUNCTION sbc_dcy_alloc 
    4949 
  • NEMO/trunk/src/OCE/SBC/sbcflx.F90

    r10068 r10425  
    157157         END DO 
    158158         taum(:,:) = taum(:,:) * tmask(:,:,1) ; wndm(:,:) = wndm(:,:) * tmask(:,:,1) 
    159          CALL lbc_lnk( taum(:,:), 'T', 1. )   ;   CALL lbc_lnk( wndm(:,:), 'T', 1. ) 
     159         CALL lbc_lnk( 'sbcflx', taum(:,:), 'T', 1. )   ;   CALL lbc_lnk( 'sbcflx', wndm(:,:), 'T', 1. ) 
    160160 
    161161         IF( nitend-nit000 <= 100 .AND. lwp ) THEN                ! control print (if less than 100 time-step asked) 
  • NEMO/trunk/src/OCE/SBC/sbcfwb.F90

    r10068 r10425  
    7171      REAL(wp), ALLOCATABLE, DIMENSION(:,:) ::   ztmsk_neg, ztmsk_pos, z_wgt ! 2D workspaces 
    7272      REAL(wp), ALLOCATABLE, DIMENSION(:,:) ::   ztmsk_tospread, zerp_cor    !   -      - 
     73      REAL(wp)   ,DIMENSION(1) ::   z_fwfprv   
     74      COMPLEX(wp),DIMENSION(1) ::   y_fwfnow   
    7375      !!---------------------------------------------------------------------- 
    7476      ! 
     
    8688         IF( kn_fwb == 3 .AND. ln_isfcav    )   CALL ctl_stop( 'sbc_fwb: nn_fwb = 3 with ln_isfcav = .TRUE. not working, we stop ' ) 
    8789         ! 
    88          area = glob_sum( e1e2t(:,:) * tmask(:,:,1))           ! interior global domain surface 
     90         area = glob_sum( 'sbcfwb', e1e2t(:,:) * tmask(:,:,1))           ! interior global domain surface 
    8991         ! isf cavities are excluded because it can feedback to the melting with generation of inhibition of plumes 
    9092         ! and in case of no melt, it can generate HSSW. 
     
    102104         ! 
    103105         IF( MOD( kt-1, kn_fsbc ) == 0 ) THEN 
    104             z_fwf = glob_sum( e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf(:,:) - snwice_fmass(:,:) ) ) / area   ! sum over the global domain 
    105             zcoef = z_fwf * rcp 
    106             emp(:,:) = emp(:,:) - z_fwf              * tmask(:,:,1) 
     106            y_fwfnow(1) = local_sum( e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf(:,:) - snwice_fmass(:,:) ) ) 
     107            CALL mpp_delay_sum( 'sbcfwb', 'fwb', y_fwfnow(:), z_fwfprv(:), kt == nitend - nn_fsbc + 1 ) 
     108            z_fwfprv(1) = z_fwfprv(1) / area 
     109            zcoef = z_fwfprv(1) * rcp 
     110            emp(:,:) = emp(:,:) - z_fwfprv(1)        * tmask(:,:,1) 
    107111            qns(:,:) = qns(:,:) + zcoef * sst_m(:,:) * tmask(:,:,1) ! account for change to the heat budget due to fw correction 
    108112         ENDIF 
     
    127131            a_fwb_b = a_fwb                           ! mean sea level taking into account the ice+snow 
    128132                                                      ! sum over the global domain 
    129             a_fwb   = glob_sum( e1e2t(:,:) * ( sshn(:,:) + snwice_mass(:,:) * r1_rau0 ) ) 
     133            a_fwb   = glob_sum( 'sbcfwb', e1e2t(:,:) * ( sshn(:,:) + snwice_mass(:,:) * r1_rau0 ) ) 
    130134            a_fwb   = a_fwb * 1.e+3 / ( area * rday * 365. )     ! convert in Kg/m3/s = mm/s 
    131135!!gm        !                                                      !!bug 365d year  
     
    154158            WHERE( erp < 0._wp )   ztmsk_pos = 0._wp 
    155159            ztmsk_neg(:,:) = tmask_i(:,:) - ztmsk_pos(:,:) 
    156             ! 
    157             zsurf_neg = glob_sum( e1e2t(:,:)*ztmsk_neg(:,:) )  ! Area filled by <0 and >0 erp  
    158             zsurf_pos = glob_sum( e1e2t(:,:)*ztmsk_pos(:,:) ) 
    159160            !                                                  ! fwf global mean (excluding ocean to ice/snow exchanges)  
    160             z_fwf     = glob_sum( e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf(:,:) - snwice_fmass(:,:) ) ) / area 
     161            z_fwf     = glob_sum( 'sbcfwb', e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf(:,:) - snwice_fmass(:,:) ) ) / area 
    161162            !             
    162163            IF( z_fwf < 0._wp ) THEN         ! spread out over >0 erp area to increase evaporation 
    163                 zsurf_tospread      = zsurf_pos 
    164                 ztmsk_tospread(:,:) = ztmsk_pos(:,:) 
     164               zsurf_pos = glob_sum( 'sbcfwb', e1e2t(:,:)*ztmsk_pos(:,:) ) 
     165               zsurf_tospread      = zsurf_pos 
     166               ztmsk_tospread(:,:) = ztmsk_pos(:,:) 
    165167            ELSE                             ! spread out over <0 erp area to increase precipitation 
    166                 zsurf_tospread      = zsurf_neg 
    167                 ztmsk_tospread(:,:) = ztmsk_neg(:,:) 
     168               zsurf_neg = glob_sum( 'sbcfwb', e1e2t(:,:)*ztmsk_neg(:,:) )  ! Area filled by <0 and >0 erp  
     169               zsurf_tospread      = zsurf_neg 
     170               ztmsk_tospread(:,:) = ztmsk_neg(:,:) 
    168171            ENDIF 
    169172            ! 
    170             zsum_fwf   = glob_sum( e1e2t(:,:) * z_fwf )         ! fwf global mean over <0 or >0 erp area 
     173            zsum_fwf   = glob_sum( 'sbcfwb', e1e2t(:,:) * z_fwf )         ! fwf global mean over <0 or >0 erp area 
    171174!!gm :  zsum_fwf   = z_fwf * area   ???  it is right?  I think so.... 
    172175            z_fwf_nsrf =  zsum_fwf / ( zsurf_tospread + rsmall ) 
    173176            !                                                  ! weight to respect erp field 2D structure  
    174             zsum_erp   = glob_sum( ztmsk_tospread(:,:) * erp(:,:) * e1e2t(:,:) ) 
     177            zsum_erp   = glob_sum( 'sbcfwb', ztmsk_tospread(:,:) * erp(:,:) * e1e2t(:,:) ) 
    175178            z_wgt(:,:) = ztmsk_tospread(:,:) * erp(:,:) / ( zsum_erp + rsmall ) 
    176179            !                                                  ! final correction term to apply 
     
    178181            ! 
    179182!!gm   ===>>>>  lbc_lnk should be useless as all the computation is done over the whole domain ! 
    180             CALL lbc_lnk( zerp_cor, 'T', 1. ) 
     183            CALL lbc_lnk( 'sbcfwb', zerp_cor, 'T', 1. ) 
    181184            ! 
    182185            emp(:,:) = emp(:,:) + zerp_cor(:,:) 
  • NEMO/trunk/src/OCE/SBC/sbcice_cice.F90

    r10068 r10425  
    100100      !!---------------------------------------------------------------------- 
    101101      ALLOCATE( png(jpi,jpj,jpnij), STAT=sbc_ice_cice_alloc ) 
    102       IF( lk_mpp                 )   CALL mpp_sum ( sbc_ice_cice_alloc ) 
     102      CALL mpp_sum ( 'sbcice_cice', sbc_ice_cice_alloc ) 
    103103      IF( sbc_ice_cice_alloc > 0 )   CALL ctl_warn('sbc_ice_cice_alloc: allocation of arrays failed.') 
    104104   END FUNCTION sbc_ice_cice_alloc 
     
    217217      ENDDO 
    218218 
    219       CALL lbc_lnk_multi( fr_iu , 'U', 1.,  fr_iv , 'V', 1. ) 
     219      CALL lbc_lnk_multi( 'sbcice_cice', fr_iu , 'U', 1.,  fr_iv , 'V', 1. ) 
    220220 
    221221      ! set the snow+ice mass 
     
    513513         ENDDO 
    514514      ENDDO 
    515       CALL lbc_lnk( ss_iou , 'U', -1. ) 
     515      CALL lbc_lnk( 'sbcice_cice', ss_iou , 'U', -1. ) 
    516516 
    517517! y comp of ocean-ice stress  
     
    525525         ENDDO 
    526526      ENDDO 
    527       CALL lbc_lnk( ss_iov , 'V', -1. ) 
     527      CALL lbc_lnk( 'sbcice_cice', ss_iov , 'V', -1. ) 
    528528 
    529529! x and y comps of surface stress 
     
    578578      fmmflx(:,:) = ztmp1(:,:) !!Joakim edit 
    579579       
    580       CALL lbc_lnk_multi( emp , 'T', 1., sfx , 'T', 1. ) 
     580      CALL lbc_lnk_multi( 'sbcice_cice', emp , 'T', 1., sfx , 'T', 1. ) 
    581581 
    582582! Solar penetrative radiation and non solar surface heat flux 
     
    604604#endif 
    605605      qsr(:,:)=qsr(:,:)+ztmp1(:,:) 
    606       CALL lbc_lnk( qsr , 'T', 1. ) 
     606      CALL lbc_lnk( 'sbcice_cice', qsr , 'T', 1. ) 
    607607 
    608608      DO jj=1,jpj 
     
    619619      qns(:,:)=qns(:,:)+nfrzmlt(:,:)+ztmp1(:,:) 
    620620 
    621       CALL lbc_lnk( qns , 'T', 1. ) 
     621      CALL lbc_lnk( 'sbcice_cice', qns , 'T', 1. ) 
    622622 
    623623! Prepare for the following CICE time-step 
     
    639639      ENDDO 
    640640 
    641       CALL lbc_lnk_multi( fr_iu , 'U', 1., fr_iv , 'V', 1. ) 
     641      CALL lbc_lnk_multi( 'sbcice_cice', fr_iu , 'U', 1., fr_iv , 'V', 1. ) 
    642642 
    643643      ! set the snow+ice mass 
     
    863863!     A. Ensure all haloes are filled in NEMO field (pn) 
    864864 
    865       CALL lbc_lnk( pn , cd_type, psgn ) 
     865      CALL lbc_lnk( 'sbcice_cice', pn , cd_type, psgn ) 
    866866 
    867867#if defined key_nemocice_decomp 
     
    10401040!     D. Ensure all haloes are filled in pn 
    10411041 
    1042       CALL lbc_lnk( pn , cd_type, psgn ) 
     1042      CALL lbc_lnk( 'sbcice_cice', pn , cd_type, psgn ) 
    10431043 
    10441044   END SUBROUTINE cice2nemo 
  • NEMO/trunk/src/OCE/SBC/sbcisf.F90

    r10068 r10425  
    157157 
    158158         ! lbclnk 
    159          CALL lbc_lnk_multi( risf_tsc(:,:,jp_tem), 'T', 1., risf_tsc(:,:,jp_sal), 'T', 1., fwfisf,'T', 1., qisf, 'T', 1.) 
     159         CALL lbc_lnk_multi( 'sbcisf', risf_tsc(:,:,jp_tem), 'T', 1., risf_tsc(:,:,jp_sal), 'T', 1., fwfisf,'T', 1., qisf, 'T', 1.) 
    160160         ! output 
    161161         IF( iom_use('iceshelf_cea') )   CALL iom_put( 'iceshelf_cea', -fwfisf(:,:)                      )   ! isf mass flux 
     
    244244               &    STAT= sbc_isf_alloc ) 
    245245         ! 
    246          IF( lk_mpp             )   CALL mpp_sum ( sbc_isf_alloc ) 
    247          IF( sbc_isf_alloc /= 0 )   CALL ctl_warn('sbc_isf_alloc: failed to allocate arrays.') 
     246         CALL mpp_sum ( 'sbcisf', sbc_isf_alloc ) 
     247         IF( sbc_isf_alloc /= 0 )   CALL ctl_stop( 'STOP', 'sbc_isf_alloc: failed to allocate arrays.' ) 
    248248         ! 
    249249      ENDIF 
     
    721721            END DO 
    722722         END DO 
    723          CALL lbc_lnk_multi( pgt, 'T', 1., pgs, 'T', 1.) 
     723         CALL lbc_lnk_multi( 'sbcisf', pgt, 'T', 1., pgs, 'T', 1.) 
    724724      END SELECT 
    725725      ! 
     
    778778            END DO 
    779779         END DO 
    780          CALL lbc_lnk(pvarout,'T',-1.) 
     780         CALL lbc_lnk('sbcisf', pvarout,'T',-1.) 
    781781       
    782782      CASE ( 'V' ) ! compute V in the top boundary layer at T- point  
     
    810810            END DO 
    811811         END DO 
    812          CALL lbc_lnk(pvarout,'T',-1.) 
     812         CALL lbc_lnk('sbcisf', pvarout,'T',-1.) 
    813813 
    814814      CASE ( 'T' ) ! compute T in the top boundary layer at T- point  
  • NEMO/trunk/src/OCE/SBC/sbcmod.F90

    r10190 r10425  
    115115      IF(lwm) WRITE( numond, namsbc ) 
    116116      ! 
     117#if defined key_mpp_mpi 
     118      ncom_fsbc = nn_fsbc    ! make nn_fsbc available for lib_mpp 
     119#endif 
    117120      !                             !* overwrite namelist parameter using CPP key information 
    118121#if defined key_agrif 
     
    440443         ! icebergs may advect into haloes during the icb step and alter emp. 
    441444         ! A lbc_lnk is necessary here to ensure restartability (#2113) 
    442          IF( .NOT. ln_passive_mode ) CALL lbc_lnk( emp, 'T', 1. ) ! ensure restartability with icebergs 
     445         IF( .NOT. ln_passive_mode ) CALL lbc_lnk( 'sbcmod', emp, 'T', 1. ) ! ensure restartability with icebergs 
    443446      ENDIF 
    444447 
     
    457460!!$!RBbug do not understand why see ticket 667 
    458461!!$!clem: it looks like it is necessary for the north fold (in certain circumstances). Don't know why. 
    459 !!$      CALL lbc_lnk( emp, 'T', 1. ) 
     462!!$      CALL lbc_lnk( 'sbcmod', emp, 'T', 1. ) 
    460463      ! 
    461464      IF( kt == nit000 ) THEN                          !   set the forcing field at nit000 - 1    ! 
  • NEMO/trunk/src/OCE/SBC/sbcrnf.F90

    r10068 r10425  
    8383         &      rnf_tsc_b(jpi,jpj,jpts) , rnf_tsc (jpi,jpj,jpts) , STAT=sbc_rnf_alloc ) 
    8484         ! 
    85       IF( lk_mpp            )   CALL mpp_sum ( sbc_rnf_alloc ) 
     85      CALL mpp_sum ( 'sbcrnf', sbc_rnf_alloc ) 
    8686      IF( sbc_rnf_alloc > 0 )   CALL ctl_warn('sbc_rnf_alloc: allocation of arrays failed') 
    8787   END FUNCTION sbc_rnf_alloc 
     
    423423         IF( nn_rnf_depth_file == 1 ) THEN      !  save  output nb levels for runoff 
    424424            IF(lwp) WRITE(numout,*) '   ==>>>   create runoff depht file' 
    425             CALL iom_open  ( TRIM( sn_dep_rnf%clname ), inum, ldwrt = .TRUE., kiolib = jprstlib ) 
     425            CALL iom_open  ( TRIM( sn_dep_rnf%clname ), inum, ldwrt = .TRUE. ) 
    426426            CALL iom_rstput( 0, 0, inum, 'rodepth', h_rnf ) 
    427427            CALL iom_close ( inum ) 
  • NEMO/trunk/src/OCE/SBC/sbcssm.F90

    r10068 r10425  
    1818   USE sbcapr         ! surface boundary condition: atmospheric pressure 
    1919   USE eosbn2         ! equation of state and related derivatives 
     20   USE traqsr, ONLY: ln_traqsr 
    2021   ! 
    2122   USE in_out_manager ! I/O manager 
     
    253254      ENDIF 
    254255      ! 
     256      IF( .NOT. ln_traqsr )   fraqsr_1lev(:,:) = 1._wp   ! default definition: qsr 100% in the fisrt level  
     257      ! 
    255258      IF( lwxios.AND.nn_fsbc > 1 ) THEN 
    256259         CALL iom_set_rstw_var_active('nn_fsbc') 
  • NEMO/trunk/src/OCE/SBC/sbcwave.F90

    r10068 r10425  
    215215      ENDIF 
    216216 
    217       CALL lbc_lnk_multi( usd, 'U', -1., vsd, 'V', -1. ) 
     217      CALL lbc_lnk_multi( 'sbcwave', usd, 'U', -1., vsd, 'V', -1. ) 
    218218 
    219219      ! 
     
    240240#endif 
    241241      ! 
    242       CALL lbc_lnk( ze3divh, 'T', 1. ) 
     242      CALL lbc_lnk( 'sbcwave', ze3divh, 'T', 1. ) 
    243243      ! 
    244244      IF( ln_linssh ) THEN   ;   ik = 1   ! none zero velocity through the sea surface 
     
    301301            END DO 
    302302         END DO 
    303          CALL lbc_lnk_multi( utau(:,:), 'U', -1. , vtau(:,:), 'V', -1. , taum(:,:) , 'T', -1. ) 
     303         CALL lbc_lnk_multi( 'sbcwave', utau(:,:), 'U', -1. , vtau(:,:), 'V', -1. , taum(:,:) , 'T', -1. ) 
    304304      ENDIF 
    305305      ! 
Note: See TracChangeset for help on using the changeset viewer.