Skip to content

Instantly share code, notes, and snippets.

View giorgk's full-sized avatar
🏖️
Working on cool projects

Georgios Kourakos giorgk

🏖️
Working on cool projects
View GitHub Profile
@giorgk
giorgk / readCVHMcbc.m
Created August 21, 2024 14:12
Read Cell by Cell Modflow
function CBC = readCVHMcbc(filename)
fid = fopen(filename,'r');
ind1 = 0;
ind2 = 1;
kper = 0;
while ~feof(fid)
temp=[];
KSTP = fread(fid,1,'uint');
KPER = fread(fid,1,'uint');
@giorgk
giorgk / ichnos_test_list.md
Last active November 4, 2022 11:10
Ichnos Tests

Ichnos test list

This is a minimum list of tests that each new version of Ichnos should run succesfully.

The ones with tick have been succesful on the current development version

Current version 0.4.10

  • CLOUD
  • Steady state
@giorgk
giorgk / Faceids_4_Modflow_example.m
Created February 15, 2022 09:04
Faceids for Modflow example
faceVel=nan(15000,22);
faceVel(1,:) = 0;
cnt_fc = 1;
msh_face_ids = zeros(size(msh,1),4);
cnt = 1;
for ii = 1:size(msh_grid,1)
for jj = 1:size(msh_grid,2)
nd = msh_grid(ii,jj).ND;
% Face 1 is the ii-1,jj FRONT
if ii == 1
@giorgk
giorgk / inv_quad.m
Created August 23, 2012 11:03
Inverse transform for quadrilateral elements
[x y]=ginput(4);
x=100*x;
y=100*y;
subplot(1,2,1);plot(x([1:4 1]),y([1:4 1]))
hold on
subplot(1,2,2);plot([-1 -1 1 1 -1],[-1 1 1 -1 -1]);
hold on
subplot(1,2,1);title('Hit enter to exit')
while 1
p1=[];p2=[];