Error while running modified LW _Test.tcl script

80 views
Skip to first unread message

Musaib Ahmad

unread,
Nov 30, 2024, 11:58:20 PM11/30/24
to par...@googlegroups.com

Dear ParFlow Users,

I am currently working on a simulation in ParFlow by attempting to run a modified LW_test.tcl script  and have encountered an error that I am unable to resolve. Despite multiple attempts at troubleshooting, I have been unable to pinpoint the root cause of the issue.

To provide context, I have attached the relevant input files used in my simulation as well as the error message I receive.

I am reaching out to this community for guidance. If anyone has encountered a similar error or has suggestions on potential solutions, I would greatly appreciate your help.

Please let me know if further details are needed. I look forward to any insights you can share.

Thank you in advance for your time and assistance.


Thanks and Regards
Musaib Ahmad
Research Scholar
IIT PALAKKAD
LW_Test.tcl
drv_vegp.dat
drv_clmin.dat
drv_vegm.alluv.dat
NLDAS.APCP.000001_to_000006.pfb
slopex.pfb
slopey.pfb
error message.png

Georgios Artavanis

unread,
Dec 2, 2024, 12:50:15 PM12/2/24
to ParFlow
Hello Musaib,

the error you are getting (line 205 at amps_ffopen.c) indicates that there is an issue with the distribution file for some of your inputs. It could be that the the path to the inputs is not correct (e.g. wrong path or typo in the path) or that you haven't distributed your inputs before the run.

Best,
George

smith84

unread,
Dec 3, 2024, 8:04:31 PM12/3/24
to ParFlow
Hello Musaib,

The problem appears to be what George suggested.  The original LW_Test.tcl script is setup to copy files from several locations to the directory ("Outputs") where the run is done.   Did you have files in correct locations for copying?   I modified to copy files from the current working directory where the tclsh command is executed.   Also CLM requires additional files beyond just the NLDAS.ACP file.   Your modified script is only distributing the NLDAS.ACP file.

Here is a slightly edited version of your script that works for me up until CLM fails when trying to read the NLDAS.DSWR.000001_to_000006.pfb.dis file.  You will need NLDAS files for the other values as in the original LW_Test.tcl example.   You need more than just the NLDAS.ACP.file.  Here is example from LW_Test.tcl:

foreach file "NLDAS.DSWR.000001_to_000024 NLDAS.DLWR.000001_to_000024 NLDAS.APCP.000001_to_000024 NLDAS.Temp.000001_to_000024 NLDAS.UGRD.000001_to_000024 NLDAS.VGRD.0000
01_to_000024 NLDAS.Press.000001_to_000024 NLDAS.SPFH.000001_to_000024"
    file copy -force [format "%s/%s.pfb" $path $file] .
    pfdist -nz 6 [format "%s.pfb" $file]
}
Thanks for reporting the problem, I have added some additional information to the output to identify the file that is causing the issue.  I will be pushing up a code change for this shortly.    In the future the <runname>.out.txt file should contain some information on what file is being read and not just exit with little information.   Here is an example of what the new log shows on my modified script:

Node 0: Using process grid (1,1,1)
Node 0: Well Information
Node 0: No Wells.
Node 0: AMPS Error: Can't open the distribution file ./NLDAS.DSWR.000001_to_000006.pfb.dist for reading
Node 0:             Common issue is incorrect name or file location
Node 0:             Common issue is not distributing the file before running
*>

This shows a bit more clearly what file is causing a problem.

Steve
LW_Test_with_edits.tcl

Musaib Ahmad

unread,
Dec 5, 2024, 11:34:45 AM12/5/24
to smith84, ParFlow

Dear Steve and George,

Thank you so much for your prompt and insightful response. Your suggestion worked perfectly! The issue was indeed due to missing input files like DWLR, DSWR, UGRD, etc. After addressing this, the simulation proceeded without errors.

For now, my objective is to simulate precipitation (APCP) only, so I temporarily had ignored the other parameters like DWLR, DSWR, and UGRD. I believe it is possible to simulate precipitation directly in ParFlow without incorporating CLM (please correct me if I’m mistaken). However, I am unsure of the specific TCL script lines required to model variable rainfall intensities for different days or hours directly within ParFlow.

Could you kindly guide me on how to implement variable rainfall intensities in ParFlow? Any example scripts or pointers would be greatly appreciated.

Thank you once again for your assistance and expertise.

Best regards,
Musaib Ahmad
Research Scholar, IIT Palakkad


--
You received this message because you are subscribed to the Google Groups "ParFlow" group.
To unsubscribe from this group and stop receiving emails from it, send an email to parflow+u...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/parflow/17e66441-588a-4995-b138-3793db8f6832n%40googlegroups.com.
Reply all
Reply to author
Forward
0 new messages