21 |
v1{iFace}=reshape(mygrid.gcm2faces{iFace}*v0,[size(mygrid.XC{iFace}) n3 n4 n5]); |
v1{iFace}=reshape(mygrid.gcm2faces{iFace}*v0,[size(mygrid.XC{iFace}) n3 n4 n5]); |
22 |
end; |
end; |
23 |
|
|
24 |
if mygrid.nFaces==1; gridType='ll'; |
v1=gcmfaces(v1); |
|
elseif mygrid.nFaces==5; gridType='llc'; |
|
|
elseif mygrid.nFaces==6; gridType='cube'; |
|
|
end; |
|
|
v1=gcmfaces(v1,gridType); |
|
25 |
|
|
26 |
elseif ~doGcm2Faces&mygrid.gcm2facesFast; |
elseif ~doGcm2Faces&mygrid.gcm2facesFast; |
27 |
|
|
60 |
|
|
61 |
end; |
end; |
62 |
|
|
63 |
if mygrid.nFaces==1; gridType='ll'; |
v1=gcmfaces(v1); |
|
elseif mygrid.nFaces==5; gridType='llc'; |
|
|
elseif mygrid.nFaces==6; gridType='cube'; |
|
|
end; |
|
|
v1=gcmfaces(v1,gridType); |
|
64 |
|
|
65 |
else; |
else; |
66 |
|
|