Dear Ferreters,
I am facing a problem in converting the depth axis from sigma coordinates to the "meter" axis. The formula that is provided in the model output is shown below.
double lev(lev) ;
lev:long_name = "ocean sigma coordinate" ;
lev:positive = "up" ;
lev:axis = "Z" ;
lev:bounds = "lev_bnds" ;
lev:formula = "z(n,k,j,i) = eta(n,j,i) + sigma(k)*(depth(j,i)+eta(n,j,i))" ; !!!!! formula !!!!!
the structure of the file is as follows