Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

remove 3d variable capability in getsfcensmeanp - not needed and broken #9

Open
wants to merge 6 commits into
base: develop
Choose a base branch
from
Open

remove 3d variable capability in getsfcensmeanp - not needed and broken #9

Changes from all commits
Commits
Show all changes
6 commits
Select commit Hold shift + click to select a range
50bafbf
remove 3d variable capability - not needed and broken
jswhit2 Jul 18, 2022
6103a6a
Merge remote-tracking branch 'upstream/develop' into develop
jswhit2 Aug 8, 2022
0aa4ddd
Merge remote-tracking branch 'upstream/develop' into develop
jswhit2 Mar 11, 2023
1c29b1d
Merge remote-tracking branch 'upstream/develop' into develop
jswhit2 Mar 15, 2024
3bf702a
Merge remote-tracking branch 'upstream/develop' into develop
jswhit2 May 15, 2024
80a329a
Merge remote-tracking branch 'upstream/develop' into develop
May 29, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
17 changes: 3 additions & 14 deletions src/EnKF/gfs/src/getsfcensmeanp.fd/getsfcensmeanp.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -42,9 +42,8 @@ program getsfcensmeanp
logical:: nemsio, sfcio, ncio

type(Dataset) :: dset,dseto
type(Dimension) :: londim,latdim,levdim
type(Dimension) :: londim,latdim
real(4), allocatable, dimension(:,:) :: values_2d, values_2d_avg
real(4), allocatable, dimension(:,:,:):: values_3d, values_3d_avg

character*500 filenamein,filenameout,datapath,fileprefix
character*3 charnanal
Expand Down Expand Up @@ -309,12 +308,10 @@ program getsfcensmeanp
elseif (ncio) then
londim = get_dim(dset,'grid_xt'); lonb = londim%len
latdim = get_dim(dset,'grid_yt'); latb = latdim%len
levdim = get_dim(dset,'pfull'); levs = levdim%len
allocate(values_2d_avg(lonb,latb))
allocate(values_3d_avg(lonb,latb,levs))
if (mype == 0) then
dseto = create_dataset(filenameout, dset, copy_vardata=.true.)
print *,'opened netcdf file ',trim(filenameout)
print *,'opened netcdf file ',trim(filenameout),lonb,latb,levs
endif
do nvar=1,dset%nvars
! Following fields are not averaged
Expand All @@ -339,22 +336,14 @@ program getsfcensmeanp
print *,'writing ens mean ',trim(dset%variables(nvar)%name)
call write_vardata(dseto,trim(dset%variables(nvar)%name),values_2d_avg)
endif
elseif (dset%variables(nvar)%ndims == 4) then
call read_vardata(dset,trim(dset%variables(nvar)%name),values_3d)
call mpi_allreduce(values_3d,values_3d_avg,lonb*latb*levs,mpi_real4,mpi_sum,new_comm,iret)
values_3d_avg = values_3d_avg * rnanals
if (mype == 0) then
print *,'writing ens mean ',trim(dset%variables(nvar)%name)
call write_vardata(dseto,trim(dset%variables(nvar)%name),values_3d_avg)
endif
else
write(6,*)'***ERROR*** invalid ndims= ',dset%variables(nvar)%ndims,&
' for variable ',trim(dset%variables(nvar)%name)
call MPI_Abort(MPI_COMM_WORLD,98,iret)
endif
enddo
if (mype == 0) call close_dataset(dseto)
deallocate(values_2d,values_2d_avg,values_3d,values_3d_avg)
deallocate(values_2d,values_2d_avg)
endif
! Jump here if more mpi processors than files to process
else
Expand Down