Changeset 1856


Ignore:
Timestamp:
2010-05-03T12:32:10+02:00 (10 years ago)
Author:
smasson
Message:

fldread_3D: small bugfix and style

Location:
branches/DEV_r1784_3DF/NEMO/OPA_SRC
Files:
9 edited

Legend:

Unmodified
Added
Removed
  • branches/DEV_r1784_3DF/NEMO/OPA_SRC/DTA/dtasal.F90

    r1806 r1856  
    6868       
    6969      !! * Local declarations 
    70       
    71       INTEGER ::   ji, jj, jk, jl, jkk   ! dummy loop indicies 
    72       INTEGER ::   & 
    73            imois, iman, i15, ik           ! temporary integers 
    74       INTEGER            :: ierror 
     70      INTEGER ::   ji, jj, jk, jl, jkk            ! dummy loop indicies 
     71      INTEGER ::   imois, iman, i15 , ik          ! temporary integers 
     72      INTEGER ::   ierror 
    7573#if defined key_tradmp 
    76       INTEGER ::   & 
    77           il0, il1, ii0, ii1, ij0, ij1   ! temporary integers          
    78 #endif 
    79       REAL(wp) ::   zxy, zl 
    80 #if defined key_orca_lev10 
    81       INTEGER   :: ikr, ikw, ikt, jjk 
    82       REAL(wp)  :: zfac 
    83 #endif 
    84       REAL(wp), DIMENSION(jpk) ::   & 
    85           zsaldta            ! auxiliary array for interpolation 
    86       CHARACTER(len=100) :: cn_dir          ! Root directory for location of ssr files 
    87       TYPE(FLD_N)        :: sn_sal 
    88       LOGICAL , SAVE     :: linit_sal = .FALSE. 
     74      INTEGER ::   il0, il1, ii0, ii1, ij0, ij1   ! temporary integers 
     75#endif 
     76      REAL(wp)::   zxy, zl 
     77#if defined key_orca_lev10 
     78      INTEGER ::   ikr, ikw, ikt, jjk  
     79      REAL(wp)::   zfac 
     80#endif 
     81      REAL(wp), DIMENSION(jpk) ::   zsaldta            ! auxiliary array for interpolation 
     82      CHARACTER(len=100)       :: cn_dir          ! Root directory for location of ssr files 
     83      TYPE(FLD_N)              :: sn_sal 
     84      LOGICAL , SAVE           :: linit_sal = .FALSE. 
    8985      !!---------------------------------------------------------------------- 
    9086      NAMELIST/namdta_sal/cn_dir,sn_sal 
     
    114110             CALL ctl_stop( 'dta_sal: unable to allocate sf_sal structure' )   ;   RETURN 
    115111         ENDIF 
    116          ALLOCATE( sf_sal(1)%fnow(jpi,jpj,jpk) ) 
     112#if defined key_orca_lev10 
     113         ALLOCATE( sf_sal(1)%fnow(jpi,jpj,jpkdta  ) ) 
     114         ALLOCATE( sf_sal(1)%fdta(jpi,jpj,jpkdta,2) ) 
     115#else 
     116         ALLOCATE( sf_sal(1)%fnow(jpi,jpj,jpk  ) ) 
    117117         ALLOCATE( sf_sal(1)%fdta(jpi,jpj,jpk,2) ) 
     118#endif 
    118119 
    119120         ! fill sf_sal with sn_sal and control print 
  • branches/DEV_r1784_3DF/NEMO/OPA_SRC/DTA/dtatem.F90

    r1806 r1856  
    7373 
    7474      !! * Local declarations 
    75       INTEGER ::   ji, jj, jk, jl, jkk       ! dummy loop indicies 
    76       INTEGER ::   & 
    77         imois, iman, i15 , ik      ! temporary integers 
    78       INTEGER            :: ierror 
     75      INTEGER ::   ji, jj, jk, jl, jkk            ! dummy loop indicies 
     76      INTEGER ::   imois, iman, i15 , ik          ! temporary integers 
     77      INTEGER ::   ierror 
    7978#if defined key_tradmp 
    80       INTEGER ::   & 
    81          il0, il1, ii0, ii1, ij0, ij1   ! temporary integers 
    82 #endif 
    83       REAL(wp) ::   zxy, zl 
    84 #if defined key_orca_lev10 
    85       !!!REAL(wp), DIMENSION(jpi,jpj,jpkdta,2) :: ztem 
    86       INTEGER   :: ikr, ikw, ikt, jjk  
    87       REAL(wp)  :: zfac 
    88 #endif 
    89       REAL(wp), DIMENSION(jpk) ::   & 
    90          ztemdta            ! auxiliary array for interpolation 
    91       CHARACTER(len=100) :: cn_dir          ! Root directory for location of ssr files 
    92       TYPE(FLD_N)        :: sn_tem 
    93       LOGICAL , SAVE     :: linit_tem = .FALSE. 
     79      INTEGER ::   il0, il1, ii0, ii1, ij0, ij1   ! temporary integers 
     80#endif 
     81      REAL(wp)::   zxy, zl 
     82#if defined key_orca_lev10 
     83      INTEGER ::   ikr, ikw, ikt, jjk  
     84      REAL(wp)::   zfac 
     85#endif 
     86      REAL(wp), DIMENSION(jpk) ::   ztemdta            ! auxiliary array for interpolation 
     87      CHARACTER(len=100)       ::   cn_dir             ! Root directory for location of ssr files 
     88      TYPE(FLD_N)              ::   sn_tem 
     89      LOGICAL , SAVE           ::   linit_tem = .FALSE. 
    9490      !!---------------------------------------------------------------------- 
    9591      NAMELIST/namdta_tem/cn_dir,sn_tem 
     
    107103         sn_tem = FLD_N( 'temperature',  -1.  ,  'votemper',  .false.   , .true.  ,  'yearly'  , ''       , ''         ) 
    108104 
    109          REWIND( numnam )         ! ... read in namlist namdta_tem  
     105         REWIND( numnam )            ! ... read in namlist namdta_tem  
    110106         READ( numnam, namdta_tem )  
    111107 
    112          IF(lwp) THEN              ! control print 
     108         IF(lwp) THEN                ! control print 
    113109            WRITE(numout,*) 
    114110            WRITE(numout,*) 'dta_tem : Temperature Climatology ' 
     
    121117 
    122118#if defined key_orca_lev10 
    123          ALLOCATE( sf_tem(1)%fnow(jpi,jpj,jpkdta ) 
     119         ALLOCATE( sf_tem(1)%fnow(jpi,jpj,jpkdta  ) ) 
    124120         ALLOCATE( sf_tem(1)%fdta(jpi,jpj,jpkdta,2) ) 
    125121#else 
    126          ALLOCATE( sf_tem(1)%fnow(jpi,jpj,jpk ) 
     122         ALLOCATE( sf_tem(1)%fnow(jpi,jpj,jpk  ) ) 
    127123         ALLOCATE( sf_tem(1)%fdta(jpi,jpj,jpk,2) ) 
    128124#endif 
     
    145141          
    146142#if defined key_tradmp 
    147       IF( cp_cfg == "orca"  .AND. jp_cfg == 2 ) THEN 
     143      IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN 
    148144             
    149145         !                                        ! ======================= 
  • branches/DEV_r1784_3DF/NEMO/OPA_SRC/SBC/fldread.F90

    r1824 r1856  
    4848      INTEGER , DIMENSION(2)          ::   nrec_b       ! before record (1: index, 2: second since Jan. 1st 00h of nit000 year) 
    4949      INTEGER , DIMENSION(2)          ::   nrec_a       ! after  record (1: index, 2: second since Jan. 1st 00h of nit000 year) 
    50       REAL(wp) , ALLOCATABLE, DIMENSION(:,:,: ::   fnow       ! input fields interpolated to now time step 
     50      REAL(wp) , ALLOCATABLE, DIMENSION(:,:,:  ) ::   fnow       ! input fields interpolated to now time step 
    5151      REAL(wp) , ALLOCATABLE, DIMENSION(:,:,:,:) ::   fdta       ! 2 consecutive record of input fields 
    5252      CHARACTER(len = 256)            ::   wgtname      ! current name of the NetCDF weight file acting as a key 
     
    146146!CDIR COLLAPSE 
    147147               sd(jf)%fdta(:,:,:,1) = sd(jf)%fdta(:,:,:,2) 
    148                sd(jf)%rotn(1)     = sd(jf)%rotn(2) 
     148               sd(jf)%rotn(1)       = sd(jf)%rotn(2) 
    149149            ENDIF 
    150150 
     
    204204 
    205205            ! read after data 
     206            ipk = SIZE( sd(jf)%fdta, 3 ) 
    206207            IF( LEN(TRIM(sd(jf)%wgtname)) > 0 ) THEN 
    207208               CALL wgt_list( sd(jf), kw ) 
    208                ipk = SIZE(sd(jf)%fdta,3) 
    209                DO jk = 1,ipk 
    210                   CALL fld_interp( sd(jf)%num, sd(jf)%clvar , kw , sd(jf)%fdta(:,:,jk,2) , sd(jf)%nrec_a(1) ) 
    211                ENDDO 
     209               DO jk = 1, ipk 
     210                  CALL fld_interp( sd(jf)%num, sd(jf)%clvar, kw, sd(jf)%fdta(:,:,jk,2), sd(jf)%nrec_a(1) ) 
     211               END DO 
    212212            ELSE 
    213                SELECT CASE( SIZE(sd(jf)%fdta,3) ) 
    214                CASE(1) 
     213               IF( ipk == 1 ) THEN  
    215214                  CALL iom_get( sd(jf)%num, jpdom_data, sd(jf)%clvar, sd(jf)%fdta(:,:,1,2), sd(jf)%nrec_a(1) ) 
    216                CASE(jpk) 
     215               ELSE 
    217216                  CALL iom_get( sd(jf)%num, jpdom_data, sd(jf)%clvar, sd(jf)%fdta(:,:,:,2), sd(jf)%nrec_a(1) ) 
    218                END SELECT 
     217               ENDIF 
    219218            ENDIF 
    220219            sd(jf)%rotn(2) = .FALSE. 
     
    256255                         vtmp(:,:) = 0.0 
    257256                         ! 
    258                          ipk = SIZE( sd(kf)%fdta(:,:,:,nf) ,3 ) 
    259                          DO jk = 1,ipk 
     257                         DO jk = 1, SIZE( sd(kf)%fdta, 3 ) 
    260258                            CALL rot_rep( sd(jf)%fdta(:,:,jk,nf),sd(kf)%fdta(:,:,jk,nf),'T', 'en->i', utmp(:,:) ) 
    261259                            CALL rot_rep( sd(jf)%fdta(:,:,jk,nf),sd(kf)%fdta(:,:,jk,nf),'T', 'en->j', vtmp(:,:) ) 
    262260                            sd(jf)%fdta(:,:,jk,nf) = utmp(:,:) 
    263261                            sd(kf)%fdta(:,:,jk,nf) = vtmp(:,:) 
    264                          ENDDO 
     262                         END DO 
    265263                         ! 
    266264                         sd(jf)%rotn(nf) = .TRUE. 
     
    335333      INTEGER :: inrec          ! number of record existing for this variable 
    336334      INTEGER :: kwgt 
    337       INTEGER :: jk             !vertical loop variable 
    338       INTEGER :: ipk            !number of vertical levels of sdjf%fdta ( 2D: ipk=1 ; 3D: ipk=jpk ) 
     335      INTEGER :: jk             ! vertical loop variable 
     336      INTEGER :: ipk            ! number of vertical levels of sdjf%fdta ( 2D: ipk=1 ; 3D: ipk=jpk ) 
    339337      CHARACTER(LEN=1000) ::   clfmt   ! write format 
    340338      !!--------------------------------------------------------------------- 
     
    401399 
    402400         ! read before data into sdjf%fdta(:,:,2) because we will swap data in the following part of fld_read 
     401         ipk = SIZE( sdjf%fdta, 3 ) 
    403402         IF( LEN(TRIM(sdjf%wgtname)) > 0 ) THEN 
    404403            CALL wgt_list( sdjf, kwgt ) 
    405             ipk = SIZE(sdjf%fdta,3) 
    406             DO jk = 1,ipk 
    407                CALL fld_interp( sdjf%num,sdjf%clvar,kwgt,sdjf%fdta(:,:,jk,2),sdjf%nrec_a(1) ) 
    408             ENDDO 
     404            DO jk = 1, ipk 
     405               CALL fld_interp( sdjf%num, sdjf%clvar, kwgt, sdjf%fdta(:,:,jk,2), sdjf%nrec_b(1) ) 
     406            END DO 
    409407         ELSE 
    410             SELECT CASE ( SIZE(sdjf%fdta,3) ) 
    411             CASE(1) 
    412                 CALL iom_get( sdjf%num, jpdom_data, sdjf%clvar, sdjf%fdta(:,:,1,2), sdjf%nrec_b(1) ) 
    413             CASE(jpk) 
    414                 if(lwp)write(numout,*)'cbr00 ',sdjf%num,SIZE(sdjf%fdta,3)  
    415                 CALL iom_get( sdjf%num, jpdom_data, sdjf%clvar, sdjf%fdta(:,:,:,2), sdjf%nrec_b(1) ) 
    416             END SELECT 
     408            IF( ipk == 1 ) THEN 
     409               CALL iom_get( sdjf%num, jpdom_data, sdjf%clvar, sdjf%fdta(:,:,1,2), sdjf%nrec_b(1) ) 
     410            ELSE 
     411               CALL iom_get( sdjf%num, jpdom_data, sdjf%clvar, sdjf%fdta(:,:,:,2), sdjf%nrec_b(1) ) 
     412            ENDIF 
    417413         ENDIF 
    418414         sdjf%rotn(2) = .FALSE. 
     
    562558      ELSE 
    563559         ! build the new filename if climatological data 
    564          IF( sdjf%cltype == 'monthly' )   WRITE(sdjf%clname, '(a,"m" ,i2.2)' ) TRIM( sdjf%clrootname ), kmonth   ! add month 
     560         IF( sdjf%cltype == 'monthly' )   WRITE(sdjf%clname, '(a,"_m",i2.2)' ) TRIM( sdjf%clrootname ), kmonth   ! add month 
    565561      ENDIF 
    566562      CALL iom_open( sdjf%clname, sdjf%num, ldstop = ldstop, ldiof =  LEN(TRIM(sdjf%wgtname)) > 0 ) 
  • branches/DEV_r1784_3DF/NEMO/OPA_SRC/SBC/sbcblk_clio.F90

    r1806 r1856  
    162162 
    163163         DO ifpr= 1, jpfld 
    164             ALLOCATE( sf(ifpr)%fnow(jpi,jpj,1) ) 
     164            ALLOCATE( sf(ifpr)%fnow(jpi,jpj,1  ) ) 
    165165            ALLOCATE( sf(ifpr)%fdta(jpi,jpj,1,2) ) 
    166166         END DO 
     
    541541               &                / ( ztatm(ji,jj) - 35.86  + MAX( 0.e0, zmt3 ) )  ) 
    542542 
    543             zev = sf(jp_humi)%fnow(ji,jj,1) * zes                      ! vapour pressure   
     543            zev = sf(jp_humi)%fnow(ji,jj,1) * zes                    ! vapour pressure   
    544544            zevsqr(ji,jj) = SQRT( zev * 0.01 )                       ! square-root of vapour pressure 
    545545            zqatm(ji,jj) = 0.622 * zev / ( zpatm - 0.378 * zev )     ! specific humidity  
     
    639639      p_qns(:,:,:) = z_qlw (:,:,:) - z_qsb (:,:,:) - p_qla (:,:,:)      ! Downward Non Solar flux 
    640640!CDIR COLLAPSE 
    641       p_tpr(:,:)   = sf(jp_prec)%fnow(:,:,1) / rday                       ! total precipitation [kg/m2/s] 
     641      p_tpr(:,:)   = sf(jp_prec)%fnow(:,:,1) / rday                     ! total precipitation [kg/m2/s] 
    642642      ! 
    643643!!gm : not necessary as all input data are lbc_lnk... 
     
    741741            zes = 611.0 * EXP(  ABS( ztamr ) * MIN ( zmt1, zmt2 )   &              ! Saturation water vapour 
    742742               &                     / ( sf(jp_tair)%fnow(ji,jj,1) - 35.86  + MAX( 0.e0, zmt3 ) )  ) 
    743             zev(ji,jj) = sf(jp_humi)%fnow(ji,jj,1) * zes * 1.0e-05                   ! vapour pressure   
     743            zev(ji,jj) = sf(jp_humi)%fnow(ji,jj,1) * zes * 1.0e-05                 ! vapour pressure   
    744744         END DO 
    745745      END DO 
     
    814814         DO ji = 1, jpi 
    815815            zlmunoon = ASIN( zps(ji,jj) + zpc(ji,jj) ) / rad                         ! local noon solar altitude 
    816             zcldcor  = MIN(  1.e0, ( 1.e0 - 0.62 * sf(jp_ccov)%fnow(ji,jj,1)   &       ! cloud correction (Reed 1977) 
     816            zcldcor  = MIN(  1.e0, ( 1.e0 - 0.62 * sf(jp_ccov)%fnow(ji,jj,1)   &     ! cloud correction (Reed 1977) 
    817817               &                          + 0.0019 * zlmunoon )                 ) 
    818             pqsr_oce(ji,jj) = zcoef1 * zcldcor * pqsr_oce(ji,jj) * tmask(ji,jj,1)   ! and zcoef1: ellipsity 
     818            pqsr_oce(ji,jj) = zcoef1 * zcldcor * pqsr_oce(ji,jj) * tmask(ji,jj,1)    ! and zcoef1: ellipsity 
    819819         END DO 
    820820      END DO 
     
    871871            zes = 611.0 * EXP(  ABS( ztamr ) * MIN ( zmt1, zmt2 )   &              ! Saturation water vapour 
    872872               &                     / ( sf(jp_tair)%fnow(ji,jj,1) - 35.86  + MAX( 0.e0, zmt3 ) )  ) 
    873             zev(ji,jj) = sf(jp_humi)%fnow(ji,jj,1) * zes * 1.0e-05                   ! vapour pressure   
     873            zev(ji,jj) = sf(jp_humi)%fnow(ji,jj,1) * zes * 1.0e-05                 ! vapour pressure   
    874874         END DO 
    875875      END DO 
  • branches/DEV_r1784_3DF/NEMO/OPA_SRC/SBC/sbcblk_core.F90

    r1806 r1856  
    164164         ENDIF 
    165165         DO ifpr= 1, jfld 
    166             ALLOCATE( sf(ifpr)%fnow(jpi,jpj,1) ) 
     166            ALLOCATE( sf(ifpr)%fnow(jpi,jpj,1  ) ) 
    167167            ALLOCATE( sf(ifpr)%fdta(jpi,jpj,1,2) ) 
    168168         END DO 
     
    262262      ! ocean albedo assumed to be 0.066 
    263263!CDIR COLLAPSE 
    264       qsr (:,:) = ( 1. - 0.066 ) * sf(jp_qsr)%fnow(:,:,1) * tmask(:,:,1)                                 ! Short Wave 
     264      qsr (:,:) = ( 1. - 0.066 ) * sf(jp_qsr)%fnow(:,:,1) * tmask(:,:,1)                                     ! Short Wave 
    265265!CDIR COLLAPSE 
    266266      zqlw(:,:) = (  sf(jp_qlw)%fnow(:,:,1) - Stef * zst(:,:)*zst(:,:)*zst(:,:)*zst(:,:)  ) * tmask(:,:,1)   ! Long  Wave 
     
    463463               ! ... scalar wind at T-point (fld being at T-point) 
    464464               zwndi_t = sf(jp_wndi)%fnow(ji,jj,1) - 0.25 * (  pui(ji,jj+1) + pui(ji+1,jj+1)   & 
    465                   &                                        + pui(ji,jj  ) + pui(ji+1,jj  )  ) 
     465                  &                                          + pui(ji,jj  ) + pui(ji+1,jj  )  ) 
    466466               zwndj_t = sf(jp_wndj)%fnow(ji,jj,1) - 0.25 * (  pvi(ji,jj+1) + pvi(ji+1,jj+1)   & 
    467                   &                                        + pvi(ji,jj  ) + pvi(ji+1,jj  )  ) 
     467                  &                                          + pvi(ji,jj  ) + pvi(ji+1,jj  )  ) 
    468468               z_wnds_t(ji,jj)  = SQRT( zwndi_t * zwndi_t + zwndj_t * zwndj_t ) * tmask(ji,jj,1) 
    469469            END DO 
     
    489489         DO jj = 2, jpjm1 
    490490            DO ji = fs_2, fs_jpim1   ! vect. opt. 
    491                p_taui(ji,jj) = zcoef_wnorm2 * ( z_wnds_t(ji+1,jj) + z_wnds_t(ji,jj) )                          & 
     491               p_taui(ji,jj) = zcoef_wnorm2 * ( z_wnds_t(ji+1,jj  ) + z_wnds_t(ji,jj) )                          & 
    492492                  &          * ( 0.5 * (sf(jp_wndi)%fnow(ji+1,jj,1) + sf(jp_wndi)%fnow(ji,jj,1) ) - pui(ji,jj) ) 
    493                p_tauj(ji,jj) = zcoef_wnorm2 * ( z_wnds_t(ji,jj+1) + z_wnds_t(ji,jj) )                          & 
     493               p_tauj(ji,jj) = zcoef_wnorm2 * ( z_wnds_t(ji,jj+1  ) + z_wnds_t(ji,jj) )                          & 
    494494                  &          * ( 0.5 * (sf(jp_wndj)%fnow(ji,jj+1,1) + sf(jp_wndj)%fnow(ji,jj,1) ) - pvi(ji,jj) ) 
    495495            END DO 
     
    517517               p_qsr(ji,jj,jl) = ( 1. - palb(ji,jj,jl) ) * sf(jp_qsr)%fnow(ji,jj,1) * tmask(ji,jj,1) 
    518518               ! Long  Wave (lw) 
    519                z_qlw(ji,jj,jl) = 0.95 * (  sf(jp_qlw)%fnow(ji,jj,1)       &                          
    520                   &                   - Stef * pst(ji,jj,jl) * zst3  ) * tmask(ji,jj,1) 
     519               z_qlw(ji,jj,jl) = 0.95 * (  sf(jp_qlw)%fnow(ji,jj,1) - Stef * pst(ji,jj,jl) * zst3  ) * tmask(ji,jj,1) 
    521520               ! lw sensitivity 
    522521               z_dqlw(ji,jj,jl) = zcoef_dqlw * zst3                                                
  • branches/DEV_r1784_3DF/NEMO/OPA_SRC/SBC/sbcflx.F90

    r1806 r1856  
    126126         ENDIF 
    127127         DO ji= 1, jpfld 
    128             ALLOCATE( sf(ji)%fnow(jpi,jpj,1) ) 
     128            ALLOCATE( sf(ji)%fnow(jpi,jpj,1  ) ) 
    129129            ALLOCATE( sf(ji)%fdta(jpi,jpj,1,2) ) 
    130130         END DO 
  • branches/DEV_r1784_3DF/NEMO/OPA_SRC/SBC/sbcice_if.F90

    r1806 r1856  
    8181            CALL ctl_stop( 'sbc_ice_if: unable to allocate sf_ice structure' )   ;   RETURN 
    8282         ENDIF 
    83          ALLOCATE( sf_ice(1)%fnow(jpi,jpj,1) ) 
     83         ALLOCATE( sf_ice(1)%fnow(jpi,jpj,1  ) ) 
    8484         ALLOCATE( sf_ice(1)%fdta(jpi,jpj,1,2) ) 
    8585 
     
    107107               ! 
    108108               zt_fzp  = fr_i(ji,jj)                        ! freezing point temperature 
    109                zfr_obs = sf_ice(1)%fnow(ji,jj,1)              ! observed ice cover 
     109               zfr_obs = sf_ice(1)%fnow(ji,jj,1)            ! observed ice cover 
    110110               !                                            ! ocean ice fraction (0/1) from the freezing point temperature 
    111111               IF( sst_m(ji,jj) <= zt_fzp ) THEN   ;   fr_i(ji,jj) = 1.e0 
  • branches/DEV_r1784_3DF/NEMO/OPA_SRC/SBC/sbcrnf.F90

    r1806 r1856  
    7575               CALL ctl_stop( 'sbc_rnf: unable to allocate sf_rnf structure' )   ;   RETURN 
    7676            ENDIF 
    77             ALLOCATE( sf_rnf(1)%fnow(jpi,jpj,1) ) 
     77            ALLOCATE( sf_rnf(1)%fnow(jpi,jpj,1  ) ) 
    7878            ALLOCATE( sf_rnf(1)%fdta(jpi,jpj,1,2) ) 
    7979         ENDIF 
  • branches/DEV_r1784_3DF/NEMO/OPA_SRC/SBC/sbcssr.F90

    r1806 r1856  
    115115               CALL ctl_stop( 'sbc_ssr: unable to allocate sf_sst structure' )   ;   RETURN 
    116116            ENDIF 
    117             ALLOCATE( sf_sst(1)%fnow(jpi,jpj,1) ) 
     117            ALLOCATE( sf_sst(1)%fnow(jpi,jpj,1  ) ) 
    118118            ALLOCATE( sf_sst(1)%fdta(jpi,jpj,1,2) ) 
    119119            ! 
     
    128128               CALL ctl_stop( 'sbc_ssr: unable to allocate sf_sss structure' )   ;   RETURN 
    129129            ENDIF 
    130             ALLOCATE( sf_sss(1)%fnow(jpi,jpj,1) ) 
     130            ALLOCATE( sf_sss(1)%fnow(jpi,jpj,1  ) ) 
    131131            ALLOCATE( sf_sss(1)%fdta(jpi,jpj,1,2) ) 
    132132            ! 
Note: See TracChangeset for help on using the changeset viewer.