make: f2py: 没有那个文件或目录

[YoungEnzo]

I encountered a problem while installing Pyroms, please help me! Thank you! I have already installed all the necessary libraries, but I encountered a problem when installing Scrip after downloading Pyroms.

the neccessary libraries:

from mpl_toolkits.basemap import Basemap
import lpsolve55
mpl_toolkits

I found an error while making sudo:
(base) zeng@zeng-virtual-machine:~/zeng/pyroms/pyroms/external/scrip/source$ sudo make
f2py --fcompiler=gnu95 --f90exec=gfortran --f90flags='-g -fdefault-real-8 -ffixed-form -O2 -fPIC'
-L. -L/usr/lib -lnetcdf -lnetcdff -ljpeg -lmfhdf -ldf -lhdf5 -lhdf5_fortran -lhdf5_hl -lhdf5hl_fortran -I. -I/usr/include -m scrip -c kinds_mod.o constants.o iounits.o netcdf.o grids.o remap_vars.o remap_distwgt.o remap_conserv.o remap_bilinear.o remap_bicubic.o remap_read.o remap.o timers.o remap_write.o ./pyscrip.f90
make: f2py: 没有那个文件或目录 (There is no such file or directory)
make: *** [makefile:120：scrip_ext] 错误 127

To my knowledge, f2py is included in numpy, and I have also provided variables:
export f2py=/home/zeng/anaconda3/envs/py2/lib/python2.7/site-packages/numpy/f2py/

but it didn't work..

I am not very familiar with Linux, so this error has made me crazy. Is anyone willing to help me? Thank you very much!!

[kshedstrom]

As I've said before, I have totally given up on scrip. Even if it compiles, it doesn't load.

Instead, use xesmf and the Arctic_HYCOM_GLBy example.

[YoungEnzo]

Wow, Kate, thanks for your quick reply and your contrubutions in roms or pyroms which really help me a lot in roms learning. You have absolutely lighten the world! haha..

I used the HYCOM_GLBY example to create my input file before, but recently my hard drive broke, so I had to reinstall Pyroms. I didn't encounter this error during my previous installation, so I'm worried it's a version issue. I used to think that you had to install Scrip before you could install pyroms. But after I import pyroms, the terminal did not report an error. So it seems that installing pyroms is not related to installing Scrip? If I successfully imprt pyroms, it means that I have successfully installed pyroms and can use it normally, right?

[kshedstrom]

Did you download a fresh pyroms? I told it to stop loading scrip when loading pyroms. Yes, you can use pyroms, but not the scrip remapping part of it.

[YoungEnzo]

Yes, I download the new pyroms by clone git.... And I need remapping. because I remember that when I create the input files (bc, ic, and clm), I fixed the remap*.py scripts to create input files. How can I download the scrip...Or can I copy the whole pyroms directory to my new hard drive from the broke one. I am afraid the copied directory cannot build the connection with my new hard drive. Dear Kate, could you please give me some suggestions? Thank you!

[kshedstrom]

You need a fresh conda environment with xesmf in it. At least sometimes, it has worked best to ask for that package first. It will allow you to do remapping, along the lines of the HYCOM_GLBy example. That's what I use now. If you can get xesmf to work, you don't need scrip ever again.

[YoungEnzo]

Thank you, Kate, but I wanna know, do I need modify the remap*.py? Because One neccessary library has changed from SCRIP to XESMF. Or, I can just run my previous remap*.py scripts and then they will create the input files directly, right?

[kshedstrom]

Well, no, you do have to change your scripts.

[YoungEnzo]

Dear Kate, is there some documentations which can guide me how to modify the code? Because I am also the new in python...Sorry to bother you with these simple questions.

[kshedstrom]

What can I say? I asked the original pyroms author for more documentation and he said no. He has since moved on, no longer works with ROMS.
I no longer work with ROMS.
I am also an idiot in Python. I'm currently hacking some Python code in a jupyter notebook, adding print statements. Do you know how to run a jupyter notebook? Or do a web search on the Python debugger? Try things and watch it fail, learning every time.

[YoungEnzo]

Thank you ,Kate, I didn't ues the JUPYTER before, but ater I searched it in the Internet, I know it.s kind of python IDE, seems like it can help you with solving error and coding? I will try to use it, thank you again. But before this, I'd like to downgrade the NUMPY. Because in my new OS, is UBUNTU22.04 and all the libraries I use are the latest, which may lead a crash to pyroms.... There is no reason why I can install SCRIP successfully befor, but it didn't work now.. So I guess creat the same environment is also kind of solution.. If everything is usefulless.. I can noly use the XESMF and modify the code which is a bit difficult for me.

[YoungEnzo]

Dear Kate, I am trying to use the xESMF and example Arctic_HYCOM_GLBy now. I used the example Palau_HYCOM beore, and I got the merged .nc file: HYCOM_GLBy0.08_2022_188.nc by running the scripts in Palau_HYCOM. I wanna ask that can I use the .nc file directly and then run scripts in Arctic_HYCOM_GLBy for remapping and creating bdry/ic/clm files?
And I also found that, Kate you said that we should not run the make_remap_weights_file.py if we use the Arctic_HYCOM_GLBy example, right?

[kshedstrom]

Yes, you download HYCOM on the HYCOM grid, then remap to your grid using the Arctic_HYCOM_GLBy example. This example should not have the make_remap_weights.py - if it does, I forgot to delete it.

[YoungEnzo]

Dear Kate, I'm sorry to bother you again, I encountered a problem..I want to use ROMS do the simulations from 2022-7-1 to 2022-7-7, but at the first step of the simulation, it shows the simulation time is 2022-7-7. After simulating July 7th, the mode was terminated at 2022-7-8 00:00:00. After I checked the input files, I found the initial .nc file only has data for July 7th.

So I try to re-create the input files with Pyroms, when I run the python get_hycom_GLBy0.08_salt/ssh/temp/u/v_2022.py, the terminal display：

(base) enzo@enzo:~/pyroms/inputfiles/yangjiang$ python get_hycom_GLBy0.08_temp_2022.py 
Processing file for water_temp, day 01, month 07, year 2022
Got water_temp from server...
Processing file for water_temp, day 02, month 07, year 2022
Got water_temp from server...
Processing file for water_temp, day 03, month 07, year 2022
Got water_temp from server...
Processing file for water_temp, day 04, month 07, year 2022
Got water_temp from server...
Processing file for water_temp, day 05, month 07, year 2022
Got water_temp from server...
Processing file for water_temp, day 06, month 07, year 2022
Got water_temp from server...
Processing file for water_temp, day 07, month 07, year 2022
Got water_temp from server...
Write with file data/HYCOM_GLBy0.08_temp_2022_188.nc
Done with file data/HYCOM_GLBy0.08_temp_2022_188.nc

But when I checked the HYCOM_GLBy0.08_temp_2022_188.nc, I found that the .nc file has only one ocean_time which is 44747.0 (which means 2022-7-7), and I think it should contain seven ocean_time values, from 44741.0 to 44747.0.
Ocean time 1 of 1=44747.0 days since 1900-01-01 00:00:00
I don't know what's wrong with it, and my code is shown below, dear Kate, could you please point out the reason? Thank you so much!

import matplotlib
matplotlib.use('Agg')

import numpy as np
import netCDF4
from datetime import datetime
import pyroms
import pyroms_toolbox
import sys
import pandas as pd


def create_HYCOM_file(name, time, lon, lat, z, var):

    print('Write with file %s' %name)

    #create netCDF file
    nc = netCDF4.Dataset(name, 'w', format='NETCDF3_64BIT')
    nc.Author = sys._getframe().f_code.co_name
    nc.Created = datetime.now().strftime("%Y-%m-%d %H:%M:%S")
    nc.title = 'HYCOM + NCODA Global 1/12 Analysis (GLBy0.08) 3 hourly'

    #create dimensions
    Mp, Lp = lon.shape
    N = len(z)
    nc.createDimension('lon', Lp)
    nc.createDimension('lat', Mp)
    nc.createDimension('z', N)
    nc.createDimension('ocean_time', None)

    #create variables        
    nc.createVariable('lon', 'f', ('lat', 'lon'))
    nc.variables['lon'].long_name = 'longitude'
    nc.variables['lon'].units = 'degrees_east'
    nc.variables['lon'][:] = lon

    nc.createVariable('lat', 'f', ('lat', 'lon'))
    nc.variables['lat'].long_name = 'latitude'
    nc.variables['lat'].units = 'degrees_north'
    nc.variables['lat'][:] = lat

    nc.createVariable('z', 'f', ('z'))
    nc.variables['z'].long_name = 'depth'
    nc.variables['z'].units = 'meter'
    nc.variables['z'][:] = z

    nc.createVariable('ocean_time', 'f', ('ocean_time'))
    nc.variables['ocean_time'].units = 'days since 1900-01-01 00:00:00'
    nc.variables['ocean_time'].calendar = 'LEAP'
    nc.variables['ocean_time'][0] = time

    nc.createVariable(outvarname, 'f', ('ocean_time', 'z', 'lat', 'lon'), fill_value=spval)
    nc.variables[outvarname].long_name = long_name
    nc.variables[outvarname].units = units
    nc.variables[outvarname].coordinates = 'lon lat'
    nc.variables[outvarname][0] = var

    nc.close()

    print('Done with file %s' %name)

# get HYCOM Northeast Pacific data from 2007 to 2011

year = 2022
retry='False'

invarname = 'water_temp'
outvarname = 'temp'

year_tag = '%04d' %year
#month_tag = '01'
#day_tag = '01'
#date_tag = year_tag + month_tag + day_tag

#read grid and variable attributes from the first file
url='https://tds.hycom.org/thredds/dodsC/datasets/GLBy0.08/expt_93.0/data/hindcasts/2022/hycom_glby_930_2022070112_t000_ts3z.nc'
dataset = netCDF4.Dataset(url)

lon = dataset.variables['lon'][1381:1423]
repetitions_lon = len(dataset.variables['lat'][2506:2558])
lon = np.tile(lon,(repetitions_lon,1))
lon = np.asarray(lon)

lat = dataset.variables['lat'][2506:2558]
repetitions_lat = len(dataset.variables['lon'][1381:1423])
lat = np.transpose([lat]*repetitions_lat)

z = dataset.variables['depth'][:]
#spval = dataset.variables[invarname]._FillValue
units = dataset.variables[invarname].units
long_name = dataset.variables[invarname].long_name
dataset.close()


retry_day = []
####
if year%4 == 0:
# daysinyear = 366
	daysinmonth = ([0, 31, 29, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31])
else:
# daysinyear = 365
	daysinmonth = ([0, 31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31])

# loop over daily files
for month in range(7,7+1):
    month_tag = '%02d' %(month)
    for day in range(1,7+1):
        day_tag = '%02d' %(day)
        date_tag = year_tag + month_tag + day_tag
        print('Processing file for %s, day %02d, month %02d, year %04d' %(invarname,day,month,year))
        url='http://tds.hycom.org/thredds/dodsC/datasets/GLBy0.08/expt_93.0/data/hindcasts/'+ year_tag +'/hycom_glby_930_'+ date_tag +'12_t000_ts3z.nc' #Changed URL according to previous url
 #get data from server
        try:
            dataset = netCDF4.Dataset(url)
            var = dataset.variables[invarname][0,:,2506:2558,1381:1423]
    #Range are set obtain from MATLAB
            spval = var.get_fill_value()
            dataset.close()
            print('Got %s from server...' %invarname)
        except:
            print('No file on the server... We skip this day.')
            retry_day.append(day)
            continue

  #create netCDF file
date_dash = year_tag + '-' + month_tag + '-' + day_tag
period = pd.Period(date_dash, freq='H') #Using panda library to convert date format
day = period.dayofyear #Convert to Julian day
outfile = 'data/HYCOM_GLBy0.08_%s_%04d_%03d.nc' %(outvarname,year,day)
jday = pyroms_toolbox.date2jday(datetime(year, 1, 1)) + day - 1
create_HYCOM_file(outfile, jday, lon, lat, z, var)

[kshedstrom]

Where did that script come from? I can't debug random scripts. The script I used to download HYCOM is in that Arctic_HYCOM_GLBy directory.

The initial file should have just one time - the initial time. What needs to have multiple times is the boundary file.

The HYCOM scripts download a lot of individual files. I would have day 188 for each field, then day 189 for each field, etc. I guess I didn't include the bash script to merge all the files for each day into one file for each day - but there are still separate daily files.

#!/bin/bash

#echo -n "Enter year to process: "; read year; echo $year
year=2023

leap=`echo $(($year % 4))`

if [[ $leap == 0 ]] ; then
    nday=366
else
    nday=365
fi
nday=37

#for day in $(eval echo {1..$nday}); do
for day in $(eval echo {11..$nday}); do
    day=`echo $day | awk '{printf "%03d", $1}'`
    echo $day
    ncks -O HYCOM_GLBy0.08_ssh_${year}_${day}.nc HYCOM_GLBy0.08_${year}_${day}.nc
    ncks -A HYCOM_GLBy0.08_temp_${year}_${day}.nc HYCOM_GLBy0.08_${year}_${day}.nc
    ncks -A HYCOM_GLBy0.08_salt_${year}_${day}.nc HYCOM_GLBy0.08_${year}_${day}.nc
    ncks -A HYCOM_GLBy0.08_u_${year}_${day}.nc HYCOM_GLBy0.08_${year}_${day}.nc
    ncks -A HYCOM_GLBy0.08_v_${year}_${day}.nc HYCOM_GLBy0.08_${year}_${day}.nc
    rm -f HYCOM_GLBy0.08_ssh_${year}_${day}.nc HYCOM_GLBy0.08_temp_${year}_${day}.nc HYCOM_GLBy0.08_salt_${year}_${day}.nc HYCOM_GLBy0.08_u_${year}_${day}.nc HYCOM_GLBy0.08_v_${year}_${day}.nc
done

[YoungEnzo]

Dear Kate,
①I modified the scripts followed other's recording vedios, it seems that these scripts can only download one day data. So I run get_hycom_GLBy0.08_salt/ssh/temp/u/v_2022.py for seven times to get one week data (2022-7-1~7). After running merge_HYCOM_GLBy0.08_daily.ksh, there are 7 files, HYCOM_GLBy0.08_2022_182-188.nc, then I run make_bdry_file.py, but I got seven boundary files(182-188). So, could I add all these 7 boundary files into roms.in which is similar as adding many forcings into roms.in ? Or, should I merge HYCOM_GLBy0.08_2022_182-188.nc into only one file and then run make_bdry_file.py?
② As for ic file, in my configuration, the simulation time from 2022-7-1 to 7-7, and my understanding is that I only need data of 1st July to create the initial file. Is it right?
③ At the same time, I used Arctic_HYCOM_GLBy directory to make input files. My simulation domain is (20.20 , 22.28 N ; 110.40 , 113.76 E). I modified the codes of get_hycom_GLBy0.08_salt_2019_rob.py:

import matplotlib
matplotlib.use('Agg')

import numpy as np
import netCDF4
from netCDF4 import num2date, date2num
from datetime import datetime, timedelta
import pyroms
import pyroms_toolbox
import os, sys, time

def create_HYCOM_file(name, time, lon, lat, z, vard, spval):

    print('Write with file %s' %name)

    #create netCDF file
    nc = netCDF4.Dataset(name, 'w', format='NETCDF3_64BIT')
    nc.Author = sys._getframe().f_code.co_name
    nc.Created = datetime.now().strftime("%Y-%m-%d %H:%M:%S")
    nc.title = 'HYCOM + NCODA Global 1/12 Analysis (GLBy0.08) expt_93.0'

    #create dimensions
    #Mp, Lp = lon.shape
    Lp = len(lon)
    Mp = len(lat)
    N = len(z)
    nc.createDimension('lon', Lp)
    nc.createDimension('lat', Mp)
    nc.createDimension('z', N)
    nc.createDimension('ocean_time', None)

    #create variables
    #nc.createVariable('lon', 'f', ('lat', 'lon'))
    nc.createVariable('lon', 'f', ('lon'))
    nc.variables['lon'].long_name = 'longitude'
    nc.variables['lon'].units = 'degrees_east'
    nc.variables['lon'][:] = lon

    #nc.createVariable('lat', 'f', ('lat', 'lon'))
    nc.createVariable('lat', 'f', ('lat'))
    nc.variables['lat'].long_name = 'latitude'
    nc.variables['lat'].units = 'degrees_north'
    nc.variables['lat'][:] = lat

    nc.createVariable('z', 'f', ('z'))
    nc.variables['z'].long_name = 'depth'
    nc.variables['z'].units = 'meter'
    nc.variables['z'][:] = z

    nc.createVariable('ocean_time', 'f', ('ocean_time'))
    nc.variables['ocean_time'].units = 'days since 1900-01-01 00:00:00'
    nc.variables['ocean_time'].calendar = 'gregorian'
    nc.variables['ocean_time'][0] = time

    nc.createVariable(outvarname, 'f', ('ocean_time', 'z', 'lat', 'lon'), fill_value=spval)
    nc.variables[outvarname].long_name = long_name
    nc.variables[outvarname].units = units
    nc.variables[outvarname].coordinates = 'lon lat'
    print("the vard shape is ", vard.shape)
    print("target shape for the variable is ", nc.variables[outvarname].shape)
    nc.variables[outvarname][0] = vard

    nc.close()

    print('Done with file %s' % name)

# get HYCOM Northeast Pacific data from 2019

year = 2022

invarname = 'salinity'
outvarname = 'salt'

#read grid and variable attributes from thredds server
url='https://tds.hycom.org/thredds/dodsC/datasets/GLBy0.08/expt_93.0/data/hindcasts/2022/hycom_glby_930_2022070112_t000_ts3z.nc'
dataset = netCDF4.Dataset(url)
#lon = dataset.variables['Longitude'][2100:,550:4040]
lon = dataset.variables['lon'][:]
lon_coord = np.where((110.40<=lon) & (lon<=113.76))
lonStart = lon_coord[0][0]
lonEnd = lon_coord[0][-1]+1
lon = dataset.variables['lon'][lonStart:lonEnd]
#lat = dataset.variables['Latitude'][2100:,550:4040]
# Lat 45.0+N
lat = dataset.variables['lat'][:]
lat_coord = np.where((20.20<=lat) & (lat<=22.28))
latStart = lat_coord[0][0]
latEnd = lat_coord[0][-1]+1
lat = dataset.variables['lat'][latStart:latEnd]
z = dataset.variables['depth'][:]
mt = dataset.variables['time'][:]
mtStart = mt[0]
mtEnd = mt[-1]
timeUnits = dataset.variables['time'].units
# Convert times and units to dates we can use
mtDates = num2date(mt,timeUnits)
units = dataset.variables[invarname].units
long_name = dataset.variables[invarname].long_name
dataset.close()

# loop over daily files
if year%4 == 0:
    daysinyear = 366
else:
    daysinyear = 365

# Use python dates to specify start and stop time
startTime = datetime(year, 7, 1, 0, 0)
iTime = startTime
endTime = datetime(year, 7, 7, 0, 0)
#endTime = datetime(year, 12, 31, 0, 0)
missingDates = [
]
# Convert start and end Time(s) to indicies from thredds server
#startNum = date2num(startTime,timeUnits)
#endNum = date2num(endTime,timeUnits)
# we assume numeric times within the dataset on the server are
# consecutive and that none are missing.
#if startNum < mtStart:
#    print("Requested start: %s (index=%d)" % (startTime.strftime("%Y-%m-%d %H:%M:%S"),startNum))
#    print("Dataset start: %s (index=%d)" % (mtDates[0].strftime("%Y-%m-%d %H:%M:%S"),mtStart))
#    print("Dataset reported time units: %s" % (timeUnits))
#    sys.exit("Start date is outside available time index of dataset.")
#if endNum > mtEnd:
#    print("Requested end: %s (index=%d)" % (endTime.strftime("%Y-%m-%d %H:%M:%S"),endNum))
#    print("Dataset end: %s (index=%d)" % (mtDates[-1].strftime("%Y-%m-%d %H:%M:%S"),mtEnd))
#    print("Dataset reported time units: %s" % (timeUnits))
#    sys.exit("End date is outside available time index of dataset.")

#sys.exit()
#import pdb; pdb.set_trace()

# Setup the correct indicies to request from the server
# 0 for thredds is the first index
#istartNum = int(startNum) - int(mtStart)
#iendNum = int(endNum) - int(mtStart)

# Thredds server limits amount of data that can be recalled at one time.
# For 3D fields, 13 chunks seems to work.

layerStart = [0, 3, 6, 9,12,15,18,21,24,27,30,33,36]
layerEnd =   [3, 6, 9,12,15,18,21,24,27,30,33,36,40]
layerChunks = len(layerStart)

#for day in range(1,daysinyear+1):
#for datasetIndex in range(istartNum,iendNum+1):
while iTime <= endTime:

    #nowNum = datasetIndex + int(mtStart)

    # create netCDF output filename
    #day = nowNum - int(date2num(datetime(year, 1, 1, 0, 0),timeUnits)) + 1
    day = iTime.timetuple().tm_yday
    outfile = 'data/HYCOM_GLBy0.08_%s_%04d_%03d.nc' %(outvarname,year,day)

    # sometimes years have specific days missing
    if iTime in missingDates:
        iTime = iTime + timedelta(days=1)
        continue

    # skip files that exist
    if os.path.isfile(outfile):
        iTime = iTime + timedelta(days=1)
        continue

    print('Processing file for %s, %s' % (invarname, iTime.strftime("%Y-%m-%d")))

    #get data from server
    url='https://tds.hycom.org/thredds/dodsC/datasets/GLBy0.08/expt_93.0/data/hindcasts/%04d/hycom_glby_930_%s12_t000_ts3z.nc' % (year, iTime.strftime("%Y%m%d"))
    writeFile = True
    for layerChunk in range(0,layerChunks):

        retries = 0
        success = False

        while success == False:
            try:
                dataset = netCDF4.Dataset(url)
                var = dataset.variables[invarname][:,layerStart[layerChunk]:layerEnd[layerChunk],latStart:latEnd,:]
                #print(var.shape)
                # strip off time dimension
                var = np.squeeze(var)

                # Server sometimes delivers blocks back with all zeros.
                # If server is returning zeros, try to give it a [90*number_of_failures second] break.
                vmin = var.min()
                vmax = var.max()
                if vmin == 0.0 and vmax == 0.0:
                    print("Field arrived as all zeros, retrying...")
                    success = False
                    retries = retries + 1
                    time.sleep(90*retries)
                    continue
                if layerChunk == 0:
                  vard = var
                  spval = var.get_fill_value()
                else:
                  #print(vard.shape,var.shape)
                  vard = np.concatenate((vard,var),axis=0)
                #print(vard.shape)
                dataset.close()
                print('Got %s(chunk %d of %d) from server...' % (invarname,layerChunk+1,layerChunks))
                success = True
            except:
                # On any server failure, wait 120 seconds.
                err = sys.exc_info()[1]
                #import pdb; pdb.set_trace()
                # If it is a missing data file, automatically skip it
                if 'strerror' in dir(err) and err.strerror == 'NetCDF: file not found':
                    #iTime = iTime + timedelta(days=1)
                    writeFile = False
                    success = True
                else:
                    print("Unexpected error:",err)
                    print('Server failed to deliver data, retrying...')
                    retries = retries + 1
                    time.sleep(120)
                    continue

    if writeFile:
        #create netCDF file
        jday = pyroms_toolbox.date2jday(datetime(year, 1, 1)) + day - 1
        create_HYCOM_file(outfile, jday, lon, lat, z, vard, spval)
    else:
        print("File not found on server %s" % (url))

    iTime = iTime + timedelta(days=1)

After downloading data for 1st July , It reported an error. When I re-run it, it can download the data for 2nd July. But there is no data in the both two output file, dear Kate, could you please tell what the matter? Many thanks for your kind and warm help.

(base) enzo@enzo:~/pyroms/inputfiles/try$ python get_hycom_GLBy0.08_salt_2019_rob.py 
Processing file for salinity, 2022-07-01
Got salinity(chunk 1 of 13) from server...
Got salinity(chunk 2 of 13) from server...
Got salinity(chunk 3 of 13) from server...
Got salinity(chunk 4 of 13) from server...
Got salinity(chunk 5 of 13) from server...
Got salinity(chunk 6 of 13) from server...
Got salinity(chunk 7 of 13) from server...
Got salinity(chunk 8 of 13) from server...
Got salinity(chunk 9 of 13) from server...
Got salinity(chunk 10 of 13) from server...
Got salinity(chunk 11 of 13) from server...
Got salinity(chunk 12 of 13) from server...
Got salinity(chunk 13 of 13) from server...
Write with file data/HYCOM_GLBy0.08_salt_2022_182.nc
Traceback (most recent call last):
  File "src/netCDF4/_netCDF4.pyx", line 5490, in netCDF4._netCDF4.Variable.__setitem__
  File "/home/enzo/anaconda3/lib/python3.10/site-packages/numpy/ma/core.py", line 3429, in shape
    super(MaskedArray, type(self)).shape.__set__(self, shape)
ValueError: cannot reshape array of size 9360000 into shape (1,40,52,42)

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
  File "/home/enzo/pyroms/inputfiles/try/get_hycom_GLBy0.08_salt_2019_rob.py", line 224, in <module>
    create_HYCOM_file(outfile, jday, lon, lat, z, vard, spval)
  File "/home/enzo/pyroms/inputfiles/try/get_hycom_GLBy0.08_salt_2019_rob.py", line 59, in create_HYCOM_file
    nc.variables[outvarname][0] = vard
  File "src/netCDF4/_netCDF4.pyx", line 5492, in netCDF4._netCDF4.Variable.__setitem__
  File "<__array_function__ internals>", line 180, in broadcast_to
  File "/home/enzo/anaconda3/lib/python3.10/site-packages/numpy/lib/stride_tricks.py", line 412, in broadcast_to
    return _broadcast_to(array, shape, subok=subok, readonly=True)
  File "/home/enzo/anaconda3/lib/python3.10/site-packages/numpy/lib/stride_tricks.py", line 348, in _broadcast_to
    it = np.nditer(
ValueError: operands could not be broadcast together with remapped shapes [original->remapped]: (40,52,4500)  and requested shape (1,40,52,42)
(base) enzo@enzo:~/pyroms/inputfiles/try$ python get_hycom_GLBy0.08_salt_2019_rob.py 
Processing file for salinity, 2022-07-02
Got salinity(chunk 1 of 13) from server...
Got salinity(chunk 2 of 13) from server...
Got salinity(chunk 3 of 13) from server...
Got salinity(chunk 4 of 13) from server...
Got salinity(chunk 5 of 13) from server...
Got salinity(chunk 6 of 13) from server...
Got salinity(chunk 7 of 13) from server...
Got salinity(chunk 8 of 13) from server...
Got salinity(chunk 9 of 13) from server...
Got salinity(chunk 10 of 13) from server...
Got salinity(chunk 11 of 13) from server...
Got salinity(chunk 12 of 13) from server...
Got salinity(chunk 13 of 13) from server...
Write with file data/HYCOM_GLBy0.08_salt_2022_183.nc
the vard shape is  (40, 52, 4500)
target shape for the variable is  (1, 40, 52, 42)
Traceback (most recent call last):
  File "src/netCDF4/_netCDF4.pyx", line 5490, in netCDF4._netCDF4.Variable.__setitem__
  File "/home/enzo/anaconda3/lib/python3.10/site-packages/numpy/ma/core.py", line 3429, in shape
    super(MaskedArray, type(self)).shape.__set__(self, shape)
ValueError: cannot reshape array of size 9360000 into shape (1,40,52,42)

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
  File "/home/enzo/pyroms/inputfiles/try/get_hycom_GLBy0.08_salt_2019_rob.py", line 226, in <module>
    create_HYCOM_file(outfile, jday, lon, lat, z, vard, spval)
  File "/home/enzo/pyroms/inputfiles/try/get_hycom_GLBy0.08_salt_2019_rob.py", line 61, in create_HYCOM_file
    nc.variables[outvarname][0] = vard
  File "src/netCDF4/_netCDF4.pyx", line 5492, in netCDF4._netCDF4.Variable.__setitem__
  File "<__array_function__ internals>", line 180, in broadcast_to
  File "/home/enzo/anaconda3/lib/python3.10/site-packages/numpy/lib/stride_tricks.py", line 412, in broadcast_to
    return _broadcast_to(array, shape, subok=subok, readonly=True)
  File "/home/enzo/anaconda3/lib/python3.10/site-packages/numpy/lib/stride_tricks.py", line 348, in _broadcast_to
    it = np.nditer(
ValueError: operands could not be broadcast together with remapped shapes [original->remapped]: (40,52,4500)  and requested shape (1,40,52,42)

I am very sorry to trouble you with reading so much content, and I just want to express my sincere gratitude once again.