Changeset 1242 for trunk/NEMO/OPA_SRC/SBC/sbcblk_core.F90

Timestamp:

2009-01-06T11:21:38+01:00 (15 years ago)

Author:

rblod

Message:

Fix runtime issues with AGRIF on NEC and add the ability to run without sea-ice on the fine grid hierarchy

File:

• trunk/NEMO/OPA_SRC/SBC/sbcblk_core.F90 (modified) (3 diffs)

trunk/NEMO/OPA_SRC/SBC/sbcblk_core.F90

@@ -166 +166 @@
 
 
-      CALL fld_read( kt, nn_fsbc, sf )                ! input fields provided at the current time-step
+      CALL fld_read( kt, nn_fsbc, sf )                   ! input fields provided at the current time-step
 
 #if defined key_lim3
@@ -173 +173 @@
 
       IF( MOD( kt - 1, nn_fsbc ) == 0 ) THEN
-          CALL blk_oce_core( sst_m, ssu_m, ssv_m )    ! compute the surface ocean fluxes using CLIO bulk formulea
+          CALL blk_oce_core( sf, sst_m, ssu_m, ssv_m )   ! compute the surface ocean fluxes using CLIO bulk formulea
       ENDIF
-      !                                               ! using CORE bulk formulea
+      !                                                  ! using CORE bulk formulea
    END SUBROUTINE sbc_blk_core
    
    
-   SUBROUTINE blk_oce_core( pst, pu, pv )
+   SUBROUTINE blk_oce_core( sf, pst, pu, pv )
       !!---------------------------------------------------------------------
       !!                     ***  ROUTINE blk_core  ***
@@ -196 +196 @@
       !!              - tprecip : Total precipitation                   (Kg/m2/s)
       !!              - sprecip : Solid precipitation                   (Kg/m2/s)
+      !!
+      !!  ** Nota  :   sf has to be a dummy argument for AGRIF on NEC
       !!---------------------------------------------------------------------
-      REAL(wp), INTENT(in), DIMENSION(jpi,jpj) ::   pst   ! surface temperature                      [Celcius]
-      REAL(wp), INTENT(in), DIMENSION(jpi,jpj) ::   pu    ! surface current at U-point (i-component) [m/s]
-      REAL(wp), INTENT(in), DIMENSION(jpi,jpj) ::   pv    ! surface current at V-point (j-component) [m/s]
+      TYPE(fld), INTENT(in), DIMENSION(:)       ::   sf    ! input data
+      REAL(wp),  INTENT(in), DIMENSION(jpi,jpj) ::   pst   ! surface temperature                      [Celcius]
+      REAL(wp),  INTENT(in), DIMENSION(jpi,jpj) ::   pu    ! surface current at U-point (i-component) [m/s]
+      REAL(wp),  INTENT(in), DIMENSION(jpi,jpj) ::   pv    ! surface current at V-point (j-component) [m/s]
 
       INTEGER  ::   ji, jj     ! dummy loop indices