Changeset 4147 for branches/2013/dev_LOCEAN_2013/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_mfs.F90

Timestamp:

2013-11-04T12:51:55+01:00 (11 years ago)

Author:

cetlod

Message:

merge in dev_LOCEAN_2013, the 1st development branch dev_r3853_CNRS9_Confsetting, from its starting point ( r3853 ) on the trunk: see ticket #1169

File:

• branches/2013/dev_LOCEAN_2013/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_mfs.F90 (modified) (2 diffs)

branches/2013/dev_LOCEAN_2013/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_mfs.F90

@@ -104 +104 @@
       INTEGER  :: ifpr     ! dummy loop indice
       INTEGER  :: jj,ji    ! dummy loop arguments
+      INTEGER  ::   ios    ! Local integer output status for namelist read
       REAL(wp) :: act_hour
       !!--------------------------------------------------------------------
@@ -127 +128 @@
       IF( kt == nit000 ) THEN                   !  First call kt=nit000  !
          !                                      ! ====================== !
-      ALLOCATE( sh_now(jpi,jpj), catm(jpi,jpj), alonl(jpi,jpj), alatl(jpi,jpj),     &
+         ALLOCATE( sh_now(jpi,jpj), catm(jpi,jpj), alonl(jpi,jpj), alatl(jpi,jpj),     &
          &        gsst(jpi,jpj),  qbw(jpi,jpj),    ha(jpi,jpj),  elat(jpi,jpj),     &
          &        evap(jpi,jpj), STAT=ierror )
 
          IF( ierror /= 0 )   CALL ctl_warn('sbc_blk_mfs: failed to allocate arrays')
+
+         REWIND( numnam_ref )              ! Namelist namsbc_msf in reference namelist : MFS files
+         READ  ( numnam_ref, namsbc_mfs, IOSTAT = ios, ERR = 901)
+901      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_mfs in reference namelist', lwp )
+
+         REWIND( numnam_cfg )              ! Namelist namsbc_msf in configuration namelist : MFS files
+         READ  ( numnam_cfg, namsbc_mfs, IOSTAT = ios, ERR = 902 )
+902      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_mfs in configuration namelist', lwp )
+         WRITE ( numond, namsbc_mfs )
+         !
+         ! store namelist information in an array
+         slf_i(jp_wndi) = sn_wndi   ;   slf_i(jp_wndj) = sn_wndj
+         slf_i(jp_clc ) = sn_clc    ;   slf_i(jp_msl ) = sn_msl
+         slf_i(jp_tair) = sn_tair   ;   slf_i(jp_rhm)  = sn_rhm
+         slf_i(jp_prec) = sn_prec   ;  
+         !
+         ALLOCATE( sf(jpfld), STAT=ierror )         ! set sf structure
+         IF( ierror > 0 ) THEN
+            CALL ctl_stop( 'sbc_blk_mfs: unable to allocate sf structure' )   ;   RETURN
+         ENDIF
+         DO ifpr= 1, jpfld
+            ALLOCATE( sf(ifpr)%fnow(jpi,jpj,1) )
+            IF( slf_i(ifpr)%ln_tint ) ALLOCATE( sf(ifpr)%fdta(jpi,jpj,1,2) )
+         END DO
+         ! fill sf with slf_i and control print
+         CALL fld_fill( sf, slf_i, cn_dir,'sbc_blk_mfs','bulk formulation for ocean SBC', 'namsbc_mfs' )
             !
-            ! set file information (default values)
-            cn_dir = './'       ! directory in which the model is executed
-            !
-            ! (NB: frequency positive => hours, negative => months)
-            !            !    file     ! frequency !  variable  ! time intep !  clim   ! 'yearly' or ! weights  ! rotation   !
-            !            !    name     !  (hours)  !   name     !   (T/F)    !  (T/F)  !  'monthly'  ! filename ! pairs      !
-            sn_wndi = FLD_N( 'ecmwf'   ,    24     ,  'u10'     ,  .false.   , .false. ,   'daily'   , ''       , ''         )
-            sn_wndj = FLD_N( 'ecmwf'   ,    24     ,  'v10'     ,  .false.   , .false. ,   'daily'   , ''       , ''         )
-            sn_clc  = FLD_N( 'ecmwf'   ,    24     ,  'clc'     ,  .false.   , .false. ,   'daily'   , ''       , ''         )
-            sn_msl  = FLD_N( 'ecmwf'   ,    24     ,  'msl'     ,  .false.   , .false. ,   'daily'   , ''       , ''         )
-            sn_tair = FLD_N( 'ecmwf'   ,    24     ,  't2'      ,  .false.   , .false. ,   'daily'   , ''       , ''         )
-            sn_rhm  = FLD_N( 'ecmwf'   ,    24     ,  'rh'      ,  .false.   , .false. ,   'daily'   , ''       , ''         )
-            sn_prec = FLD_N( 'precip_cmap' ,  -1   ,  'precip'  ,  .true.    ,  .true. ,   'yearly'  , ''       , ''         )
-            !
-            REWIND( numnam )                    ! ... read in namlist namsbc_mfs
-            READ  ( numnam, namsbc_mfs )
-            !
-            ! store namelist information in an array
-            slf_i(jp_wndi) = sn_wndi   ;   slf_i(jp_wndj) = sn_wndj
-            slf_i(jp_clc ) = sn_clc    ;   slf_i(jp_msl ) = sn_msl
-            slf_i(jp_tair) = sn_tair   ;   slf_i(jp_rhm)  = sn_rhm
-            slf_i(jp_prec) = sn_prec   ;  
-            !
-            ALLOCATE( sf(jpfld), STAT=ierror )         ! set sf structure
-            IF( ierror > 0 ) THEN
-               CALL ctl_stop( 'sbc_blk_mfs: unable to allocate sf structure' )   ;   RETURN
-            ENDIF
-            DO ifpr= 1, jpfld
-               ALLOCATE( sf(ifpr)%fnow(jpi,jpj,1) )
-               IF( slf_i(ifpr)%ln_tint ) ALLOCATE( sf(ifpr)%fdta(jpi,jpj,1,2) )
-            END DO
-            ! fill sf with slf_i and control print
-            CALL fld_fill( sf, slf_i, cn_dir,'sbc_blk_mfs','bulk formulation for ocean SBC', 'namsbc_mfs' )
-            !
-         ENDIF
+      ENDIF
 
          CALL fld_read( kt, nn_fsbc, sf )                   ! input fields provided at the current time-step