Opened 11 years ago
Closed 11 years ago
#364 closed upgrade (Done)
Updated balance operator input parameters
Reported by: | arango | Owned by: | arango |
---|---|---|---|
Priority: | major | Milestone: | Adjoint Based Algorithms |
Component: | Adjoint | Version: | 3.3 |
Keywords: | Cc: |
Description
Updated the balance operator to specify its input parameters from the s4dvar.in script:
! Number of iterations in the biconjugate gradient algorithm used to solve ! the elliptic equation for sea surface height in the error covariance ! balance operator, [1:Ngrids]. Nbico == 200 ! Parameters used to compute balanced salinity in terms of temperature using ! empirical T-S relationships in the error covariance balance operator, ! [1:Ngrids]. dTdz_min == 0.001d0 ! minimun dT/dz (Celsius/m) ml_depth == 100.0d0 ! mixed-layer depth (m; positive)
Notice that we need as many Nbico iterations are required to decrease the error value of the reference free-surface to 1E-8 or smaller. In many applications Nbico=200 will do the job.
Warning: be aware that there are 4 arrays that are allocated using Nbico and its value maybe constrained by available memory:
- FOURDVAR(ng) % p_r2d (LBi:UBi,LBj:UBj,Nbico(ng))
- FOURDVAR(ng) % r_r2d (LBi:UBi,LBj:UBj,Nbico(ng))
- FOURDVAR(nd) % bp_r2d(LBi:UBi,LBj:UBj,Nbico(ng))
- FOURDVAR(ng) % br_r2d(LBi:UBi,LBj:UBj,Nbico(ng))
All the values are needed in the backward iterations of the adjoint.
I also corrected couple of bugs. Many thanks to Andy more fore reporting and correcting these problems.
Change History (1)
comment:1 Changed 11 years ago by arango
- Resolution set to Done
- Status changed from new to closed
Note: See
TracTickets for help on using
tickets.