Changeset 13686 for NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/MY_SRC/icestp.F90

Timestamp:

2020-10-28T10:56:07+01:00 (3 years ago)

Author:

laurent

Message:

Improvements for "oce+ice capable" STATION_ASF!

File:

• NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/MY_SRC/icestp.F90 (modified) (5 diffs)

NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/MY_SRC/icestp.F90

@@ -74 +74 @@
    !!----------------------------------------------------------------------
    !! NEMO/ICE 4.0 , NEMO Consortium (2018)
-   !! $Id: icestp.F90 13641 2020-10-19 18:16:58Z clem $
+   !! $Id: icestp.F90 13655 2020-10-21 14:15:13Z laurent $
    !! Software governed by the CeCILL license (see ./LICENSE)
    !!----------------------------------------------------------------------
@@ -127 +127 @@
          !    utau_ice, vtau_ice = surface ice stress [N/m2]
          !------------------------------------------------!
+                                        CALL ice_sbc_tau( kt, ksbc, utau_ice, vtau_ice )
+         !
          !------------------------------------------------------!
          ! --- Thermodynamical coupling with the atmosphere --- !
@@ -150 +152 @@
          !!=============================================================
 
-         !IF(lwp) WRITE(numout,*) ''
-         !IF(lwp) WRITE(numout,*) 'LOLO:ice_stp()@icestp.F90 => calling "ice_update_flx()" at kt =', kt
-         !                               CALL ice_update_flx( kt )     ! -- Update ocean surface mass, heat and salt fluxes
-         !IF(lwp) WRITE(numout,*) 'LOLO:ice_stp()@icestp.F90 => Done with "ice_update_flx()" at kt =', kt
-         !IF(lwp) WRITE(numout,*) ''
-
          ! --- case we bypass ice thermodynamics --- !
          !IF( .NOT. ln_icethd ) THEN   ! we suppose ice is impermeable => ocean is isolated from atmosphere
@@ -196 +192 @@
       ! --- Ocean time step --- !
       !-------------------------!
-      !IF( ln_icedyn )                   CALL ice_update_tau( kt, uu(:,:,1,Kbb), vv(:,:,1,Kbb) )   ! -- update surface ocean stresses
-      !!gm   remark, the ocean-ice stress is not saved in ice diag call above .....  find a solution!!!
       !
       IF( ln_timing )   CALL timing_stop('ice_stp')
@@ -238 +232 @@
       !
       CALL ice_sbc_init                ! set ice-ocean and ice-atm. coupling parameters
+      !
+      fr_i  (:,:)   = at_i(:,:)        ! initialisation of sea-ice fraction
       !
       IF( ln_rstart )   CALL iom_close( numrir )  ! close input ice restart file