Download ISSM at Rocky Linux 8.7

aisundong

The problem seems to come from here. When I calculate using ISSM4.18, it seems that md.smb.monthlyproperties and md.smb.precipitation will automatically match numberofvertices.This can operate normally in 4.18.

But when I switched to ISSM4.23, I had to add a 'mean()' otherwise it wouldn't match, but with this change, I couldn't run PDD correctly.What I mean is that this error seems to come from some changes in the PDD between 4.18 and 4.23? Do you know exactly why?

aisundong

This is my setting for md.smb.

mathieumorlighem

OH! I see you are using ismumgsm=0 and isdelta18o = 0 which we never test... let me see if I can fix it

aisundong

When I remove the 'mean()', which means there are 12 columns as you mentioned, I will receive the following error:

model not consistent: field 'smb.monthlytemperatures' should have only one column as there are md.mesh.numberofvertices lines
model not consistent: field 'smb.precipitation' should have only one column as there are md.mesh.numberofvertices lines

How to solve this problem? How long will it take for you to complete this part of the inspection.
Thank you for your help!

mathieumorlighem

ok, could you change SMBpdd.m and implement the following change (replace 'timeseries',1 by 'size',[md.mesh.numberofvertices 12]):

- md = checkfield(md,'fieldname','smb.monthlytemperatures','timeseries',1,'NaN',1,'Inf',1);
- md = checkfield(md,'fieldname','smb.precipitation','timeseries',1,'NaN',1,'Inf',1);
+ md = checkfield(md,'fieldname','smb.monthlytemperatures','size',[md.mesh.numberofvertices 12],'NaN',1,'Inf',1);
+ md = checkfield(md,'fieldname','smb.precipitation','size',[md.mesh.numberofvertices 12],'NaN',1,'Inf',1);

that should fix your problem. Let me know!
Mathieu

aisundong

mathieumorlighem
The PDD method uses summer temperature (jja) to calculate melting. So, is it related to the order of the 12 month data? For example, when simulating the southern hemisphere, do I need to convert the summer data to the front of the sequence

aisundong

Still reporting the same error.
May I confirm that I need change this file
/lustre/home/2201110305/ISSM/trunk/src/m/classes/SMBpdd.m.
If so, it seens didn't work

mathieumorlighem

ok, to be clear, it is still reporting
model not consistent: field 'smb.monthlytemperatures' should have only one column as there are md.mesh.numberofvertices lines model not consistent: field 'smb.precipitation' should have only one column as there are md.mesh.numberofvertices lines
??

It may be that MATLAB uses the copy /lustre/home/2201110305/ISSM/trunk/bin/SMBpdd.m

aisundong

You're right!
I changed this /lustre/home/2201110305/ISSM/trunk/bin/SMBpdd.m ,and it ran successfully in the end.
Thanks a lot

mathieumorlighem

ok great! Thanks for confirming.
Mathieu

mathieumorlighem

It should be ordered normally: Jan-Feb-mar....
Cheers
Mathieu

aisundong

mathieumorlighem
Will I get the correct result if my seasonal order changes? Where can I check how the specific calculations in this part of the code are implemented

mathieumorlighem

The actual implementation is in src/c/shared/Elements/PddSurfaceMassBalance.cpp and you will see that it uses JJA temperatures (indices 6-7-8). The order needs to be 12 months starting with January.

« Previous Page