JRA3Q_Reanalysis/DownLoad_JRA3Q_anl125_Pressure_Level.sh at main · QianyeSu/JRA3Q_Reanalysis · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
#!/bin/bash

opts="-N"
cert_opt="--no-check-certificate"

# Loop through the years 1948 to 1978
for year in {1948..2023}
do
    # Loop through each month
    for month in {01..12}
    do
        # Construct the month_year string used in the base URL
        month_year="${year}${month}"
        base_url="https://osdf-director.osg-htc.org/ncar/gdex/d640000/anl_p125/${month_year}"

        # Define all variable codes and their respective GRIB variable names
        variables=(
            "0_0_0.tmp-pres-an-ll125" "0_1_0.spfh-pres-an-ll125" \
            "0_1_1.rh-pres-an-ll125" "0_2_2.ugrd-pres-an-ll125" \
            "0_2_3.vgrd-pres-an-ll125" "0_2_8.vvel-pres-an-ll125" \
            "0_3_5.hgt-pres-an-ll125" "0_0_7.depr-pres-an-ll125" \
            "0_2_12.relv-pres-an-ll125" "0_2_13.reld-pres-an-ll125" \
            "0_2_4.strm-pres-an-ll125" "0_2_5.vpot-pres-an-ll125"
        )

        # Determine the end day of the month to handle February and leap years
        if [ "$month" == "02" ]; then
            if ((year % 4 == 0 && (year % 100 != 0 || year % 400 == 0))); then
                end_day="29" # Leap year February
            else
                end_day="28" # Non-leap year February
            fi
        elif [[ "$month" == "04" || "$month" == "06" || "$month" == "09" || "$month" == "11" ]]; then
            end_day="30"
        else
            end_day="31"
        fi

        # Loop through each variable code
        for var in "${variables[@]}"
        do
            start_time="${year}${month}0100"
            end_time="${year}${month}${end_day}18"
            file_name="jra3q.anl_p125.${var}.${start_time}_${end_time}.nc"

            # Check if the file already exists and is not empty
            if [ -f "$file_name" ] && [ -s "$file_name" ]; then
                echo "File $file_name already exists and is not empty. Skipping download."
            else
                echo "Downloading $file_name..."
                # Download the file
                wget $cert_opt $opts $base_url/$file_name
            fi
        done
    done
done