Changeset 13082


Ignore:
Timestamp:
2020-06-10T09:38:08+02:00 (3 months ago)
Author:
smueller
Message:

Removal of variable fwfold from module sbcfwb (ticket #2377)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/ticket2377/src/OCE/SBC/sbcfwb.F90

    r12608 r13082  
    3737   REAL(wp) ::   a_fwb     ! annual domain averaged freshwater budget from the 
    3838                           ! previous year 
    39    REAL(wp) ::   fwfold    ! fwfold to be suppressed 
    4039   REAL(wp) ::   area      ! global mean ocean surface (interior domain) 
    4140 
     
    112111         ENDIF 
    113112         ! 
    114       CASE ( 2 )                             !==  fwf budget adjusted from the previous year  ==! 
     113      CASE ( 2 )                             !==  fw adjustment based on fw budget at the end of the previous year  ==! 
    115114         ! 
    116115         IF( kt == nit000 ) THEN                      ! initialisation 
    117             !                                         ! Read the corrective factor on precipitations (fwfold) 
     116            !                                         ! read the fw adjustment (a_fwb) 
    118117            IF ( ln_rstart .AND. iom_varid( numror, 'a_fwb',   ldstop = .FALSE. ) > 0 ) THEN 
    119118               IF(lwp) WRITE(numout,*) 'sbc_fwb : reading FW-budget adjustment from restart file' 
     
    124123               CLOSE( inum ) 
    125124            END IF 
    126             fwfold = a_fwb                            ! current year freshwater budget correction 
    127             !                                         ! estimate from the previous year budget 
     125            ! 
    128126            IF(lwp)WRITE(numout,*) 
    129             IF(lwp)WRITE(numout,*)'sbc_fwb : year = ',iyear  , ' freshwater budget correction = ', fwfold 
     127            IF(lwp)WRITE(numout,*)'sbc_fwb : year = ',iyear  , ' freshwater budget correction = ', a_fwb 
    130128            IF(lwp)WRITE(numout,*)'          year = ',iyear-1, ' freshwater budget read       = ', a_fwb 
    131129            ! 
     
    134132            END IF 
    135133         ENDIF    
    136          !                                         ! Update fwfold if new year start 
     134         !                                         ! Update a_fwb if new year start 
    137135         ikty = 365 * 86400 / rn_Dt                  !!bug  use of 365 days leap year or 360d year !!!!!!! 
    138136         IF( MOD( kt, ikty ) == 0 ) THEN 
     
    142140            a_fwb   = a_fwb * 1.e+3 / ( area * rday * 365. )     ! convert in Kg/m3/s = mm/s 
    143141!!gm        !                                                      !!bug 365d year  
    144             fwfold =  a_fwb                           ! current year freshwater budget correction 
    145             !                                         ! estimate from the previous year budget 
    146142         ENDIF 
    147143         !  
    148144         IF( MOD( kt-1, kn_fsbc ) == 0 ) THEN         ! correct the freshwater fluxes 
    149             zcoef = fwfold * rcp 
    150             emp(:,:) = emp(:,:) + fwfold             * tmask(:,:,1) 
     145            zcoef = a_fwb * rcp 
     146            emp(:,:) = emp(:,:) + a_fwb              * tmask(:,:,1) 
    151147            qns(:,:) = qns(:,:) - zcoef * sst_m(:,:) * tmask(:,:,1) ! account for change to the heat budget due to fw correction 
    152148         ENDIF 
     
    161157         END IF 
    162158         ! 
    163          IF( kt == nitend .AND. lwm ) THEN            ! save fwfold value in a file (only one required) 
     159         IF( kt == nitend .AND. lwm ) THEN            ! save a_fwb value in a file (only one required) 
    164160            CALL ctl_opn( inum, 'EMPave.dat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE., narea ) 
    165161            WRITE( inum, "(24X,I8,2ES24.16)" ) nyear, a_fwb 
Note: See TracChangeset for help on using the changeset viewer.