Changeset 8426 for branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/iceforcing.F90

Timestamp:

2017-08-08T17:53:09+02:00 (7 years ago)

Author:

clem

Message:

last routine names to be changed

File:

• branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/iceforcing.F90 (modified) (8 diffs)

branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/iceforcing.F90

@@ -19 +19 @@
    USE sbcblk          ! Surface boundary condition: bulk
    USE sbccpl          ! Surface boundary condition: coupled interface
-   USE icealbedo       ! ice albedo
+   USE icealb          ! ice albedo
    !
    USE iom             ! I/O manager library
@@ -64 +64 @@
       IF( nn_timing == 1 )   CALL timing_start('ice_forcing_tau')
 
+      IF( kt == nit000 .AND. lwp ) THEN
+         WRITE(numout,*)
+         WRITE(numout,*)'ice_forcing_tau'
+         WRITE(numout,*)'~~~~~~~~~~~~~~~'
+      ENDIF
+
       SELECT CASE( ksbc )
          CASE( jp_usr     )   ;    CALL usrdef_sbc_ice_tau( kt )                 ! user defined formulation
@@ -78 +84 @@
             END DO
          END DO
-         CALL lbc_lnk( utau_ice, 'U', -1. )
-         CALL lbc_lnk( vtau_ice, 'V', -1. )
+         CALL lbc_lnk_multi( utau_ice, 'U', -1., vtau_ice, 'V', -1. )
       ENDIF
 
@@ -119 +124 @@
       IF( nn_timing == 1 )   CALL timing_start('ice_forcing_flx')
 
+      IF( kt == nit000 .AND. lwp ) THEN
+         WRITE(numout,*)
+         WRITE(numout,*)'ice_forcing_flx'
+         WRITE(numout,*)'~~~~~~~~~~~~~~~'
+      ENDIF
+
       ! --- cloud-sky and overcast-sky ice albedos --- !
-      CALL ice_albedo( t_su, ht_i, ht_s, a_ip_frac, h_ip, ln_pnd_rad, zalb_cs, zalb_os )
+      CALL ice_alb( t_su, ht_i, ht_s, a_ip_frac, h_ip, ln_pnd_rad, zalb_cs, zalb_os )
 
       ! albedo depends on cloud fraction because of non-linear spectral effects
@@ -140 +151 @@
                                    CALL blk_ice_flx( t_su, alb_ice )
             IF( ln_mixcpl      )   CALL sbc_cpl_ice_flx( picefr=at_i_b, palbi=alb_ice, psst=sst_m, pist=t_su )
-            IF( nn_limflx /= 2 )   CALL ice_lim_flx( t_su, alb_ice, qns_ice, qsr_ice, dqns_ice, evap_ice, devap_ice, nn_limflx )
+            IF( nn_limflx /= 2 )   CALL ice_flx_dist( t_su, alb_ice, qns_ice, qsr_ice, dqns_ice, evap_ice, devap_ice, nn_limflx )
 
          CASE ( jp_purecpl )                                     ! coupled formulation
                                    CALL sbc_cpl_ice_flx( picefr=at_i_b, palbi=alb_ice, psst=sst_m, pist=t_su )
-            IF( nn_limflx == 2 )   CALL ice_lim_flx( t_su, alb_ice, qns_ice, qsr_ice, dqns_ice, evap_ice, devap_ice, nn_limflx )
+            IF( nn_limflx == 2 )   CALL ice_flx_dist( t_su, alb_ice, qns_ice, qsr_ice, dqns_ice, evap_ice, devap_ice, nn_limflx )
       END SELECT
 
@@ -162 +173 @@
 
 
-   SUBROUTINE ice_lim_flx( ptn_ice, palb_ice, pqns_ice, pqsr_ice, pdqn_ice, pevap_ice, pdevap_ice, k_limflx )
-      !!---------------------------------------------------------------------
-      !!                  ***  ROUTINE ice_lim_flx  ***
+   SUBROUTINE ice_flx_dist( ptn_ice, palb_ice, pqns_ice, pqsr_ice, pdqn_ice, pevap_ice, pdevap_ice, k_limflx )
+      !!---------------------------------------------------------------------
+      !!                  ***  ROUTINE ice_flx_dist  ***
       !!
       !! ** Purpose :   update the ice surface boundary condition by averaging and / or
@@ -195 +206 @@
       !!----------------------------------------------------------------------
       !
-      IF( nn_timing == 1 )  CALL timing_start('ice_lim_flx')
+      IF( nn_timing == 1 )  CALL timing_start('ice_flx_dist')
       !
       SELECT CASE( k_limflx )                              !==  averaged on all ice categories  ==!
@@ -231 +242 @@
       END SELECT
       !
-      IF( nn_timing == 1 )  CALL timing_stop('ice_lim_flx')
-      !
-   END SUBROUTINE ice_lim_flx
+      IF( nn_timing == 1 )  CALL timing_stop('ice_flx_dist')
+      !
+   END SUBROUTINE ice_flx_dist